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SOME CRITICISMS OF “ON THE ORIGIN OF THE LUNAR SURFACE 
FEATURES” BY G. P. KUIPER* 
By Haroup C. Urry 
INSTITUTE FOR NUCLEAR STUDIES, UNIVERSITY OF CHICAGO 
Communicated May 12, 1955 


In a recent issue of these PRocEEDINGS, Professor Kuiper! has presented some 
views in regard to the structure and history of the moon which are in sharp dis- 
agreement with views presented by the writer? * and others*: ° on these matters. 
It is the purpose of this note to answer Kuiper’s argument. 

Kuiper has made additional observations on the moon with the 82-inch reflector 
at MeDonald Observatory; but, so far as I can deduce from his discussion, he has 
not observed anything markedly different from what has been previously observed 
and recorded in the extensive literature on the moon. It may be that he has ob- 
served some new details, but the information bearing on his theories has been avail- 
able before. His discussion, which takes issue with my ideas, really derives from 
his view that the moon was completely melted in its early history because of radio- 
active heating. In the fall of 1953 I remarked to Professor Kuiper in a discussion 
of the increased age of the solar system as presented by Patterson® that it seemed 
to me that the constants, i.e., the age of the solar system and amounts of the radio- 
active elements, should be such that the moon would not have melted. From this 
remark it should have been obvious that, at that time, I had already made Dr. 
Kuiper’s calculations in regard to the melting of the moon. In a recent paper,’ 
I discussed the abundances of the radioactive isotopes relative to this question of a 
melted moon, as well as the heat balance of Mars and the earth. From this paper 
it is evident that serious difficulties in regard to the explanation of the properties 
of the earth, the moon, and Mars are encountered, unless the abundances of K, U, 
and Th are much less in the planets than they are in the chondritic meteorites. 
That such should be the case is not surprising in view of the quite obvious differ- 
entiation of the meteorites in other ways. 

Kuiper maintains that the moon’s irregular shape is not due to a frozen tidal 
wave but was produced by the irregular fall of planetesimals onto its surface and 
that this occurred while the moon was completely molten. Certainly, if the moon 
were molten, its surface would be an equipotential surface appropriate to its dis- 
tance from the earth at that time. In this case the irregular shape must be due toa 
frozen tidal wave, in contradiction to Kuiper’s other arguments. 

He argues that the moon’s elongated axis, now pointed toward and away from 
the earth, was once oriented along its orbital path. For this he uses dynamical 
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arguments which assume that the moon had a rigid structure as it moved from a 
position nearer the earth; but he also assumes that at the time when it was oriented 
with its newly acquired elongated axis in the direction of its path it was also com- 
pletely melted and hence not rigid at all. Sir George Darwin has discussed the 
question of plastic deformation of the earth and has shown convincingly that dynam- 
ical problems of this kind are very difficult to handle rigorously. Kuiper’s very 
brief discussion of this subject is at least internally inconsistent. 

Kuiper postulates that the moon accumulated at low temperatures, an assump- 
tion which had previously been made explicitly. He then points out that radio- 
active heating would raise its temperature. If the amounts of radioactive sub- 
stances in the moon are the same as those in the meteorites as reported by Patter- 
son,® melting would occur in about 5 X 10° years or more’ (the timing of events 
here is difficult to understand and will be discussed later). It seems to me that ex- 
pansion would occur during this heating process both because of thermal expansion 
and because of expansion on melting, and hence cracks would occur in the moon’s 
surface, and the liquid silicates would flow up through such cracks and flood the 
surface. All these effects are ignored in Kuiper’s presentation. In fact, it is im- 
possible to decide whether in his opinion the liquid lunar silicates are less or more 
dense than the solids. It seems most improbable that any surface features of the 
moon acquired previous to complete melting would remain after this melting proc- 
ess, as he assumes. 

Next he postulates that solidification of the moon took place with the radioactive 
elements concentrating almost entirely at the surface. Now the solid silicates are 
definitely assumed to be more dense than their liquids, and hence solidification be- 
gins to occur in the moon’s interior; this is compared with a similar process of 
solidification of the earth. In previous publications?» * I showed that the moon 
would lose heat from its interior very slowly by conduction, and a revision of these 
calculations with more recent data confirms the previous conclusions. — It is difficult 
to present briefly the serious objections to Kuiper’s assumptions. Even in 4.5 X 
10° years the center of the moon would lose very little heat through conduction and 
hence would remain near its melting point because of original heat alone, if at any 
time, for any reason, it had been completely melted. Even a small fraction of the 
postulated concentrations of radioactive elements, only 15 per cent of those re- 
ported for the chondritic meteorites, would maintain the interior of the moon out 
to 0.8 of its radius at the melting point of silicates or, indeed, cause extensive melt- 
ing. It is difficult to believe that all the radioactive elements were removed by a 
solidification process. This certainly has not occurred in the case of the earth. If 
complete separation of the radioactive elements has occurred in the case of the 
earth, the mean concentration of potassium, uranium, and thorium in the earth as a 
whole must be very low, and very little heating due to radioactive substances 
would occur in matter of this concentration. It would then be necessary to assume 
much higher concentrations in the moon and the parent-bodies of the meteorites if 
Kuiper’s postulated heating of these objects occurs. This necessitates some con- 
sideration of possible processes for making large astronomical objects of markedly 
different chemical composition, which Kuiper has not done. Even using much 
lower abundances for the radioactive elements in the earth,’ it appears that only 
one-half to one-fourth of these elements are in the earth’s crust. No consistent 
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explanation for the rigidity of the moon as indicated by its irregular shape has been 
offered by Kuiper, for the moon could hardly retain an irregular shape if its in- 
terior were melted. It seems to me that the moon accumulated at low temperatures 
and that the concentrations of radioactive elements are so low that the moon has 
never melted. I advanced this view some years ago, and there appears to be no 
other way to account for the facts. 

Kuiper believes that a ‘‘sediment ring,”’ consisting of solid objects of the kind that 
bombarded the moon early in its history, previously existed about the earth. He 
believes that this ring accumulated by molecular condensation of a “silicate gas,”’ 
i.e., molecular weight of 50, at some 100° K. during some 10?—10* years. In an- 
other publication’ he accounted for the planetesimals, previously postulated by the 
writer to explain the chemical properties of the earth and meteorites, by a similar 
mechanism of accumulation during some 10’ years. At the concentrations postu- 
lated, ‘silicate gas’’ or, in fact, water vapor, which he also discusses, would be highly 
supersaturated. Radioactive substances were present, and in their presence (or 
even in their absence) nucleation would occur in seconds, and no growth of large 
objects would occur. Kuiper says: ‘The absence of such bodies would be more 
puzzling than their presence.”” That such objects existed is certain, since objects 
did collide with the moon; but if physical-chemical theory and observation are 
valid, they could not have accumulated as he assumes. Also it should be noted 
that no such sediment rings exist now in the neighborhood of any planets except 
Saturn, and the existence of this ring seems to be due to the breakup of a satellite 
which formed or tried to form within the Roche limit of the planet. 

The maria are markedly concentrated in the center of the moon’s disk. Kuiper 
has turned the moon through 90° as it passed through his postulated “sediment 
ring” and then has the planetesimals hitting the advanced and trailing regions of 
the moon. But in most cases the moon would surely perturb the orbits of these 
planetesimals before collision occurred, and they would be scattered out of the ring 
into approximately elliptical orbits. Subsequently, collisions from all directions 
would occur, and hence this mechanism gives no explanation for the distribution of 
the maria. 

Kuiper postulates that the planetesimal which produced the Hyginus crater initi- 
ated a crack some hundreds of kilometers long in a rigid crust of the moon, a crust 
which was floating in some way on melted silicates. Subsequently, according to 
Kuiper, these liquid silicates rose in the crack to form a level floor. Kuiper’s solid 
lunar crust must have had most unusual properties, lower densities than the liquids 
and great tensile strength, in order for this floating block to exist under great tension 
and form a crack some 5-10 km. wide when ‘“‘triggered”’ by the fall of a planetesimal. 
Silicates have low tensile strength even in mechanically perfect samples, and the 
surface of the moon subjected to enormous forces by an intense bombardment would 
have the strength of sand. Gilbert® accounted for this floor of the Hyginus Cleft by 
postulating partial filling of the crack with liquid silicates produced during the col- 
lisions that formed the lunar maria. It seems to me that Gilbert’s explanation is at 
least in accord with physical possibilities, while Kuiper’s is not. 

There are other doubtfui points in Kuiper’s paper. I recently discussed the 
question of the origin of tektites.’ Kuiper suggests that they originated from the 
moon, a very old suggestion, made again recently by Nininger.'’ It is difficult to 














126 ASTRONOMY: H. C. UREY Proc, N. ALS. 


understand their chemical composition and terrestrial distribution unless they are of 
terrestrial origin. 

Melted lunar silicates might be more dense than solids with which they are in 
equilibrium, or less dense, but they should not sometimes rise in cracks only to a 
certain distance below the surface and at other times rise to the tops of high moun- 
tains, as Kuiper assumes. Also it seems most improbable that nearly all the lunar 
surface was subjected to intense bombardment and yet substantial areas were 
missed completely, so that an earlier accreted surface is preserved through the time 
of assumed melting and bombardment, as Kuiper assumes for the region north of 
Mare Imbrium. Gilbert’ maintained that vestiges of old craters could be seen in 
these areas and that they were covered with debris during the Imbrian collision, 
which seems to be a more probable explanation. 

Kuiper postulates that the Alpine Valley is a crack produced during the Imbrian 
collision and not a groove produced by an iron-nickel object plowing through the 
lunar surface, as the writer has suggested. It may be that this is the case. As he 
points out, it has a slight bend in it, which is evident in good photographs; but it is 
remarkably straight for a fracture, and, if it was produced in this way, one would 
expect other fractures in its neighborhood. If it was produced by a missile, one 
wonders where the displaced material is. One also wonders whether such a large 
crack, some 130 km. long and 10 km. wide, could be made in material of low tensile 
strength. Isn’t the problem similar to Hubbert’s'™ problem of lifting the state of 
Texas? The lava in either case could have flowed into the valley from Mare Im- 
brium and need not have risen from the interior, since the large mountainous masses 
probably did not completely close the entrance to the mare. The Alpine Valley 
is 10 km. wide and, if it is a crack, must have been equally deep when it was formed. 
In this case slumping should have occurred, causing irregularities in the two sides, 
and any corresponding features on the two sides as noted by Kuiper must be acciden- 
tal. 

Kuiper believes that the central peaks of the craters cannot be due to a rebound 
of the lunar material after a collision. He states: “While one can visualize a re- 
bound in a liquid or plastic, there seems no reason to suppose that a solid can act 
in thismanner.” I find that many physical scientists, particularly those conversant 
with ballistic problems, have observed that solids flow under high pressure. Gilbert 
in 1893 specifically pointed this out. On the other hand, Kuiper speaks of 
“squeezed-up”’ lavas, an idea which this writer is wholly unable to entertain as a 
serious physical postulate. Lavas are liquid, not putty, and could hardly form 
great, simple “functional” structures. And what sort of mechanism squeezed 
them up in the center of the craters? Kuiper’s suggestions are similar to Shaler’s, '* 
which I studied and rejected five years ago. 

No clear picture of the order of events or the times at which they occurred 
emerges from Kuiper’s paper. Probably such a picture is impossible unless some defi- 
nite conclusion has been reached in regard to the amounts of the radioactive elements 
in the substance of the moon, and Kuiper does not state what he believes these 
amounts to be. In fact without specifying such data, one cannot estimate the 
time required for melting or, indeed, whether melting would occur at all. Kuiper 
does assume that objects suitable for producing the meteorites melted 4.5 X 10° 
years ago, owing to radioactive heating. Hence substantial objects, some hundreds 
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of kilometers in radius perhaps, must have accumulated some 5 or 5.5 x 10° 
years ago, acquired a high temperature during 0.5-1.0 X 10% years, as estimated by 
Urey,’ been broken up by some collision processes, and the debris accumulated in the 
chondritic meteorites, which have remained essentially unchanged for the last 4.5 X 
10° years. But when did the earth and moon accumulate? Kuiper believes that 
the moon did not melt until some 5 x 10° years after the earth and moon were 
formed. 

It would be a thankless task to review adequately this paper in all details. Fea- 
tures of the moon that have been observed and recorded during the last century 
are again described as though they were newly observed and then are followed by « 
discussion of the completely melted moon theory. On page 1104 there is a long 
paragraph listing facts well known for many years, and then at the bottom of the 
page is the statement: “It is at first somewhat surprising that the mountain ranges 
around Mare Imbrium should have formed during the period of maximum melting, 
at which time the crust must have been fairly thin and possibly was not able to sup- 
port the extra load.”’ It is certainly surprising to this reader. Then comes a 
discussion of the Alpine Valley that still leaves me at a complete loss as to how these 
massive mountains were supported. 

It seems to me that Kuiper’s theory is not internally consistent. Perhaps a 
theory of the moon which assumes that it was completely molten in its early history 
can be devised, even though he has not done so. Such a theory would probably 
best assume that the irregular shape of the moon is due to a frozen tidal wave and 
that the thermal conductivity of lunar material is considerably higher than current 
estimates for terrestrial materials. Then it should be assumed that the irregular 
shape and the surface features were formed very late in its history when it had 
acquired a very rigid structure. 

In a recent publication’ I made an attempt to explore possible radioactive heating 
effects during the early history of the solar system. Very complicated and puzzling 
histories must be assumed. Though the process of the moon’s origin was undoubt- 
edly more intricate than anyone has the courage to imagine, it is my belief today, 
as it has been for some five years, that very probably the moon was accumulated at 
low temperatures from a primitive dust cloud of solar composition with the iron in 
oxidized states and that the concentrations of radioactive substances within the 
moon are sufficiently low that melting has never occurred. In fact, I believe that 
present temperatures are so low that the interior of the moon has a high strength 
and that such low temperatures require the moon to have been formed at low tem- 
peratures and never to have been melted at any time. 

* Presented at a conference on the abundances of the elements held at Williams Bay, Wisconsin, 
September, 1953. 
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EFFECT OF ADENINE MONONUCLEOTIDES ON NEURAL TUBE 
FORMATION OF FROG EMBRYO 


By ELISABETH AMBELLAN 
ZOOLOGY DEPARTMENT, COLUMBIA UNIVERSITY, NEW YORK, NEW YORK 
Communicated by L. C. Dunn, May 10, 1955 


This is the first of a series of studies now in progress on the biochemical nature and 
energy relationships of developmental processes of Rana pipiens treated with 
solutions of materials thought to be essential growth substances. This treatment 
has the advantage of identifying effective agents without artifacts of graft or 
transplant experiments that may confuse inducing materials with products of 
cytolysis, and avoids the necessity of extrapolating from what “‘breis’’ or explants 
may do under experimental conditions to what the intact embryo does under rela- 
tively normal conditions. 

Frog embryos require little energy during cleavage, blastulation, and early 
gastrulation,’ after which energy demands are markedly increased. This period 
is marked biochemically by the synthesis of purines,’ ribonucleic acids,* and new 
specific proteins,‘ and these are presumably connected with ensuing morphological 
processes. The measurable respiratory rate of the gastrula increases, and there 
is a high local concentration of adenosine triphosphate (ATP) in areas of synthetic 
activity.° 

One method for synthesis of the nucleic acids and new proteins that precedes 
visible differentiation requires, besides the necessary raw materials and enzymes, 
high energy sources to build larger molecules. The first experiments, therefore, 
were to test whether or not the growing embryo could utilize ATP as an exogenous 
source of energy in furtherance of normal synthetic activities. 

Methods.—Ten to twenty dejellied eggs per Stender dish in about 15 ml. of 
freshly made solutions were kept in an 18° C. constant-temperature bath. Solution 
volumes and number of eggs were constant within any one experiment. Dishes 
were taken out of the water bath for observation one at a time (in later experiments) 
and were never out of the bath for more than a few minutes. Selection of an 
apparently homogeneous starting population by eliminating any abnormal eggs 
and any that were ahead of or behind the average stage of development for the 
clutch resulted in 90-100 per cent of eggs within dishes maintaining the same rate 
of development throughout most of the experiments. 

Estimations of growth rate were first made by arbitrary substaging of Shumway 
stages (as pictured in Rugh’s manual), but in the last three experiments direct time 
measurements were made. Samples of experimentals were fixed in Bouin’s solu- 
tion, and the time taken for controls grown in 10 per cent Ringer’s solution to reach 
the same degree of neural tube closure as the fixed experimentals was considered 
the measure of advanced growth rate of experimentals over controls. 
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Schwartz (chromatographically pure) disodium salt of adenosine triphosphate 
(Na,-ATP-3H.O) and Schwartz yeast adenylic acid (AMP-3) which have a pH of 
about 3.6 in 10 per cent Ringer’s solution were brought to pH 5.6 with “ Ringer’s 
hydroxides,” a mixture of Na, K, and Ca hydroxides with the same cation concen- 
trations as Ringer’s solution. Sigma adenosine diphosphate (ADP) having neutral 
pH in 10 per cent Ringer’s solution was brought to pH 5.6 with 0.1 N HCl. In- 
organic pyrophosphate (Na4sP.O;-10H»O), H;PO,, Schwartz adenine and Schwartz 
adenosine were used in various experiments. 

Results.—EKighteen separate experiments involving treatment of R. pipiens em- 
bryos with ATP in 10 per cent Ringer’s solution at pH 5.6 for short time periods 
beginning during late blastula stages resulted in precocious formation of neural 
tubes that developed 2-3 days after removal of embryos from treatment solutions. 
Duration of treatments ranged from 12 to 24 hours, and, when the treatments 
were applied at stages 7 to late 9, the effects were consistent and specific for this 
one morphological process. When controls were at late neural fold stage (Shum- 
way, late 14) with open anterior folds 0.6 mm. wide, the neural ridges in experi- 
mentals were completely fused and rotation was observed (Shumway, late 15). 
The average advance in growth rate at this period of development was estimated 
by substaging methods at about 4 hours for all eighteen experiments; in the last 
three experiments in which this time was actually recorded, one showed 3*/,, an- 
other 4, and another 4!/2 hour’s advance. This advanced development continued 
through the next neural tube stage (Shumway, 16), gradually diminished, and by 
tail bud stage (Shumway, 17) the treated embryos usually exhibited typical abnor- 
malities of open blastopores, microcephaly, bloated bellies, and stunted bodies. 

Growth effects, toxicity, and abnormalities were generally found in proportion 
to concentration of materials, length of treatment and developmental stage at which 
given. Treatment of blastulae at pH 5.6 with 0.008 M ATP was toxic in three 
experiments and was not used again; 0.004 M was toxic to some clutches, or there 
were many casualties within experimental dishes; 0.001 ./ sometimes showed less 
advance, and 0.0005 M showed little or no advance over controls. In general, 
0.002 M gave the maximum advance in neural tube formation without early toxicity 
or later abnormalities and was used as the additional control in experiments testing 
effects of other materials where negative results were anticipated. 

Other experiments, not included in the eighteen covered by this report, showed 
that continuous ATP treatment in similar concentrations never produced addi- 
tional growth increase beyond that obtained by short-term treatment for the same 
clutch. The continuous treatment was more toxic, and the few survivors showed 
typical but exaggerated abnormalities in proportion to concentration. In another 
experiment, 0.0001 M ATP was effective with 8 hours of treatment applied right 
after first cleavage, although 0.002 and 0.004 M were toxic when applied this early. 
One experiment with R. sylvatica embryos showed that ATP, 0.002 M, pH 5.6, had 
the same effect in neural tube formation as with the R. pipiens embryos. 

The exact limit of effective pH ranges of ATP was not determined, but those 
used in most of the experiments ranged from pH 5.5 to 5.8. At pH 3.6 both con- 
trols and experimentals died. In one experiment ATP was effective at 0.001, 
0.002, and 0.004 M at pH 5.6, but at pH 7.1 the two higher concentrations were 
toxic, and the 0.001 M had only slight positive effects, followed by extreme ab- 
normalities. 
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Experiments showed that additional Na ion concentration (equivalent to the 
Nay-ATP added), inorganic phosphate, inorganic pyrophosphate, adenine, and 
adenosine in concentrations 0.001—0.008 M at pH 5.6 had no effect on neurulation. 

Three separate experiments using adenosine diphosphate (ADP) at pH 5.6 re- 
sulted in the same rate of neural tube formation as with ATP treatment, and the 
embryos could tolerate twice the concentrations without cytolysis or abnormalities. 
Three experiments using yeast adenylic acid (AMP-3) in 0.001—0.004 M concentra- 
tions resulted in the same rate of increase in neurulation, but after effects at tail bud 
stage varied between clutches. Results of one experiment with a single clutch of 
embryos treated for 20 hours beginning during late blastula and then returned to 
10 per cent Ringer’s solution (minus bicarbonate and phosphate) are given in 
Table 1. Stages in the second column were simultaneous observations. ‘‘ Hours 
advance”’ was time recorded at 18° C. for controls to reach the same stage of neural 
tube closure as a typical experimental fixed in Bouin’s solution as a sample. 


TABLE 1 


Errect oF SHort-TeERM TREATMENT WITH ADENINE MONONUCLEOTIDE SOLUTIONS ON FROG 
DEVELOPMENT 
Hours Advance 
Shumway Stage over Controls Mm. Length 
at One at Neurulation Observations at after 
Material in Solution Observation (Approx.) Tail Bud Stage Hatching 
10 per cent Ringer’s (controls )* Late 14 Normal 5.4 av. 
ATP 0.004 M Late 15 4!'/, Stubby bodies, 
microcephaly t¢ 
ATP 0.002 M Late 15 4!/, Stubby 
AMP-3 0.004 7 Late 15 4!/, Stubbyt 
AMP-3 0.002 M Late 15 41/, Stubbyt 
ADP 0.008 M Late 15 4!'/, Normal 
ADP 0.004 M Late 15 41/, Normal 
ADP 0.002 M Late 15 41/, Normal 
ATP 0.004 M + AMP-3 0.004 Late 15 4'/s Normal 
ATP 0.002 M + AMP-3 0.002 Late 15 41/, Normal 
ATP 0.001 M + AMP-3 0.001 Late 15 41/, Normal 


9 
4.9 


rnowocok & ee 


—) 


* Three replicates of controls dishes had no apparent variability in development between dishes or in individual 
eggs within dishes. 

t Only 3 eggs were left in this dish at this stage; others died during previous treatment. Ten eggs in each of all 
other control and experimental dishes lived until hatching, except those fixed in Bouin’s solution during experiment. 

t In a later experiment these tail bud stage abnormalities following advanced neurulation did not appear. 


Discussion.—The effects of the adenine mononucleotides, ATP, ADP, and 
AMP-3, were found to be generally proportional to concentration, were effective 
long after embryos were removed from treatment solutions, and were highly specific 
for one morphological process. It is assumed, therefore, that there was penetration 
of these materials that affected a specific metabolism of the embryos. While 
the effectiveness of materials applied was measured in terms of closure of the neural 
tube, such organized morphogenetic movements must result, in turn, from pre- 
vious synthetic cellular activities and must be an integral part of the constructive 
pattern of development that constitutes “growth” during this period. Since these 
mononucleotides enhance a normal morphogenetic process, it is suggested that they 
can be utilized in the normal biochemical processes causing this directed mor- 
phology. The advanced neural tube formation resulting from treatment with all 
three mononucleotides indicates that it is the adenylic acid common to all of them 
that is the effective agent—adenine, adenosine, inorganic phosphate, and pyro- 
phosphate having been experimentally eliminated from consideration when they 
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had no effect on development under similar conditions. It may be that the ef- 
fective materials found here are used as raw materials for synthesis of nucleic acids 
and new proteins in a manner similar to the role of the mononucleotides suggested 
by Brachet® in neural plate formation. Essential enzymes for this synthesis, or the 
ability to induce them when presented with the proper substrate, are assumed to 
be present in the embryo, since they are needed for development under normal 
conditions. If other mononucleotides or combinations of them are found to have 
the same effect, this will lend support to Brachet’s hypothesis.* In any case, this 
hypothesis can now be investigated further, since the specific metabolic pathway 
taken from the time the adenine mononucleotides are first introduced to their final 
effect in neural tube formation can be followed by isotope, bio-assay, or cyto- 
chemical methods. 

If the adenylic acid is the effective part of these molecules, and not the high- 
energy pyrophosphate bond in ATP, as was first thought, these experiments indi- 
cate that both energy mechanisms and nucleic acids play a role in differentiation. 
In the experiment reported in Table 1, 0.004 7 ATP was highly toxic alone but was 
not toxic when mixed with equimolar AMP-3. It is possible that the equilibrium 
found in experiments, in vitro, 


2ADP = ATP + AMP 


is here serving to maintain the energy equilibrium required for normal growth in the 
entire organism. 

The observed pH dependence of the ATP reactions could be explained as being 
due to higher penetration of less ionized molecules (expected with ATP at the 
more acid pH), but this is not likely here, since toxicity of ATP at pH 7.1 indicates 
penetrance at this pH also. Experiments with other materials used in similar or 
higher concentrations that had no effect on development showed that additional 
salts in the medium were not by themselves toxic. It may be that there is active 
transport by, or activation of, enzymes that function at this lower pH, but there 
is no evidence on this point. This pH dependence may explain Brachet’s results’ 
finding that ATP, AMP, and other nucleic acid derivatives were toxic for gastrula 
cells, “although the pH of the solution was carefully adjusted.”’ It is possible that 
his solutions were adjusted to a more alkaline pH. It is also possible that 
there is a difference in behavior between the intact organisms used here and the 
explants used by Brachet. 

Summary.—In twenty-four experiments FR. pipiens blastulae treated with solu- 
tions of ATP, ADP, or AMP-3 (with specific conditions of pH and concentrations) 
for short periods of time resulted in precocious closure of neural tubes about 4 hours 
(18° C.) ahead of controls. In one experiment with R. sylvatica, using ATP, the 
results were identical. The effects of solution treatments appeared 2-3 days after 
embryos were returned to normal medium. Adenine, adenosine, inorganic phos- 
phate, and pyrophosphate in similar experimental conditions had no effect on 
neurulation, nor had the effect of pH alone. A tentative hypothesis for the method 
of utilization of mononucleotides is presented. 

* A recent experiment using separately cytidylic, uridylic, and guanylic acids at pH 5.6, 0.004M 
for 35 hours of treatment indicates that these other mononucleotides are also effective in acceler- 
ating neural tube formation. 
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Solution treatments that enhance a normal growth process offer a simple system 
for study of biochemical factors of morphogenesis in an intact organism under 
comparatively normal conditions. The mechanism by which the adenine mono- 
nucleotides, or other materials that may be found to favor a specific growth process, 
affect metabolism after penetration can be followed further by more refined tech- 
niques. 


I wish to express thanks to Professor Lester G. Barth and Dr. Lucena J. Barth for 
numerous suggestions, helpful criticisms, manuscript revisions—in short, for in- 
valuable hours of time given to this work. 


1L. G. and L. J. Barth, Energetics of Development (New York: Columbia University Press, 
1954). 

2S. Graff and L. Barth, Cold Spring Harbor Symposia Quant. Biol., 6, 103-108, 1938. 

3M. Steinert, Biochim. et Biophys. Acta, 10, 427-431, 1953, quoted by Brachet, Arch. Biol., 
15, 44, 1954. 

4 Barth and Barth, op. cit., p. 70. 

5 J. Brachet and H. Chantrenne, Acta Biol. Belg., 4, 455, 1942, quoted by Brachet, Chemical 
Embryology (New York: Interscience Publishers, 1950), p. 403. 

6 J. Brachet, Cold Spring Harbor Symposia Quant. Biol., 13, 18-25, 1947. 

7 Ibid. 


A RAPID METHOD FOR VIABLE CELL TITRATION AND CLONE 
PRODUCTION WITH Hela CELLS IN TISSUE CULTURE: THE USE 
OF X-IRRADIATED CELLS TO SUPPLY CONDITIONING FACTORS* 


By THreopore T. Puck AND Puiuip I. Marcus 


DEPARTMENT OF BIOPHYSICS, FLORENCE R. SABIN LABORATORIES, UNIVERSITY OF COLORADO 
MEDICAL CENTER, DENVER 


Communicated by O. H. Robertson, May 5, 1955 


Studies of many aspects of the genetics and metabolism of animal cells are 
seriously handicapped by lack of a simple, effective technique for large-scale colony 
production from single cells. Thus it is not easily possible to develop a new popu- 
lation from the single-cell survivors isolated after the action of high-energy irradi- 
ation, virus invasion, or any of a variety of toxic drugs or nutritional] stresses. This 
difficulty also makes almost impossible the preparation of mutant cell strains, a 
technique which has so enormously advanced studies of cell and virus genetics and 
biochemical transformations in microbial systems.' In addition, it is not possible 
to determine what fraction of any animal cell population is able to initiate multi- 
plication, and therefore it becomes difficult to evaluate the mechanism of action of 
agents which change growth rates. 

Earle and his associates? have studied quantitatively the minimum number of 
L-strain fibroblasts required in a given volume of medium in order to insure growth. 
Their careful investigations have emphasized the necessity for ‘‘conditioning”’ by 
the living cells of the nutrient media currently employed in tissue culture before 
multiplication can be supported. They were able to produce clones from a variety 
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of cells by a fairly exacting technique involving sealing of individual cells in capillary 
tubes, where diffusible metabolic products essential for multiplication remain in 
association with the cells. However, at most, only a few per cent of the single 
cells isolated by this method could be carried to the point of self-sustaining colony 
development. In attempting to develop a rapid method for growing clones from 
every viable cell, as is routinely accomplished with bacteria, we were impressed 
by the demonstration of Earle and his colleagues of the need to conserve diffusible 
cellular products, and so directed our efforts toward furnishing isolated cells with 
maximal amounts of these substances. 


METHODS AND MATERIALS 


The HeLa strain of epithelial cell obtained by Gey et al.’ from a human cervical carcinoma 
was used throughout. These cells, initially supplied by the George Washington Carver 
Foundation, are maintained in this laboratory in continuous subeulture. Growth medium 
of the following composition was employed: synthetic solution of amino acids, vitamins, 
and other growth factors,‘ 40 per cent; horse serum, 10 per cent; human serum, 20 per cent; 
Hank’s solution, 30 per cent. Penicillin (100 units/cc) and streptomycin (100 ugm/ce) were 
routinely added. Cultures were incubated at 36.5° C. in an atmosphere of water-saturated 
air containing 3-5 per cent COs. 

Suspensions of HeLa cells were prepared by trypsinization, according to standard tech- 
niques,® of the stock cultures which are grown on glass. Cell counts were made with a 
hemacytometer, and suspensions were diluted in nutrient medium to provide the requisite 
cell densities. These standardized suspensions were used in two types of experiments. In 
one procedure a series of approximately globular droplets, with volumes varying from 0.02 
to 0.001 ce. and with progressively varying numbers of cells, was deposited on glass micro- 
scope slides, which were then placed in a Petri dish, empty except for cotton pledgets satu- 
rated with the growth medium, for humidity maintenance. These dishes were incubated, 
and the droplets subjected to daily microscopic count of cell numbers, from which growth 
rates were determined. In the other procedure, 0.05 ce. of standard cell suspensions was 
pipetted onto a piece of microscope slide (1.0 2.5 em.) and spread with the tip of the 
pipette over the entire surface, to insure separation of the single cells. Microscopic exami- 
nation of such slides confirmed that the cells so deposited lodge at widely separated points and 
in numbers to be expected within the uncertainty of sampling error. The slides were incu- 
bated for 5 hours, during which the cells attached to the glass. They were then placed in 
Petri dishes containing nutrient medium and nonmultiplying feeder cells, as described below, 
and reincubated. Fresh nutrient was exchanged for the old medium after approximately 
3-4 days. 


EXPERIMENTAL RESULTS 


1. Hstablishment of the Minimum Number of Cells Necessary To Initiate Sus- 
tained Growth in Nutrient Medium Alone.—tThe first step consisted in determination 
of the minimum colony size which would permit uniformly successful subculture and 
regrowth to any desired cell number. It was found that colonies with as little as 
25-50 cells could regularly be trypsinized and transferred to drops on new micro- 
scope slides or test tubes, where they would eventually yield confluent growth. 
By successive transfer from a slide colony to a test tube and then to a 1,000-ce. 
bottle, cell populations as great as 10’ could be developed. To provide a margin 
of safety, the figure of 100 cells was selected as a minimum colony size from which 
routine subculture leading to any desired cell number could reliably be expected. 
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Next, experiments were undertaken to ascertain conditions under which iso- 
lated single cells would grow up to colonies of 100 or more. Varying numbers of 
cells, in droplets of nutrient medium of controlled size, were deposited on gless 
slides and incubated as described, in the thought that sufficiently small droplets 
might limit dilution of any essential, diffusible molecules as effectively as a capillary. 
However, the smallest drop size that could conveniently be manipulated contained 
a volume of 0.001 cc., which was still very much greater than that of the capil- 
laries used by Earle. Such droplets regularly promoted growth past the stage of 
100 cells only if they contained an initial inoculum greater than 11 cells (Table 1). 


TABLE 1 
GrowTH oF HeLa CeLts iN Discrete Drop.ets or 0.001-0.0025 Ce. 


No. of No. of Cells in Per Cent of Tests Which Resulted 
Trials Original Inoculum in Cell Proliferation Past the 100 Mark* 
56 1-5 18 
12 6-10 66 
32 11-50 100 


* The criterion taken for unlimited growth was the development of a colony of 100 or more, since experience 
demonstrated that this number of cells could reliably be grown up to 107? by successive subculture in progressively 
larger containers. 


The use of conditioned medium, as recommended by Earle, did not produce 
appreciable improvement. The medium was harvested after having supported a 
growing cell culture for 18-24 hours. Regardless of whether such conditioned 
medium was used full strength or diluted 1:1 with fresh medium, it produced no 
more success than fresh medium alone in promoting growth of small numbers of 
cells. It was concluded that, if the bulk of the cell population is really capable of 
reproducing and if the sole reason for failure of a single cell to initiate sustained mul- 
tiplication is lack of a sufficient concentration of cell-manufactured, diffusible 
metabolites, some of these may have a half-life at 37° C. considerably less than 24 
hours. Experiments were therefore designed to take account of this possibility. 

2. Growth of Clones from Single Cells Using X-trradiated ‘Feeder’ Cells To 
Supply Conditioning Factors—Since the identity of the presumed short-lived, dif- 
fusible, conditioning factors was unknown, an attempt was made to provide a con- 
tinuing supply of these by a large number of metabolizing cells which are them- 
selves incapable of multiplication (feeder cells). If such growth-inhibited cells 
were arranged in close juxtaposition to a single multiplying cell, they could supply 
the latter with the necessary growth factors which should permit multiplication 
into a colony. Since high-energy irradiation operates to stop cell multiplication 
long before general metabolism is appreciably affected, this seemed an ideal agent 
for preparation of feeder cells to supply the needs of the immature clone. 

Details of one experimental arrangement which gives successful results are as fol- 
lows: 3 X 10° HeLa cells in 7 ml. of growth medium were pipetted into a 50-mm. 
Petri dish and incubated overnight to form the feeder layer. The following morn- 
ing the medium was removed, replaced with 1 ml. of balanced saline, and the dishes 
exposed to a beam of hard X-rays (210-40 kv.; 15-20 ma.; 0.5-mm. Cu and 1.0-mm. 
Al filters). Doses varying between 5,000 and 40,000 roentgens appeared equally 
effective. After exposure, the original medium, centrifuged to remove any con- 
tained cells, was replaced and a sterile, U-shaped piece of Plexiglas laid over the 
feeder-cell layer. Two or three segments of glass microscope slide, on which vary- 
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ing aliquots of a standardized HeLa cell suspension had previously been spread, 
were then placed on the Plexiglas holder. The arrangement employed is shown in 
Figure 1.6 
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Fic. 1.—Schematic diagram of the arrangement used to grow clones of HeLa cells. Both the 
clones (which actually achieve a diameter of about | mm. in 8-10 days) and the cells are enor- 
mously exaggerated in size. A—U-shaped plastic holder supporting the microscope slides; B— 
segments of microscope slide on which the clones develop; C—clone; D—Petri dish; E—level of 
liquid nutrient medium; /—layer of feeder cells. 
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The plates were incubated for 8-10 days. Microscopic examination during this 
period revealed that the feeder layer undergoes little apparent change during the 
first 48-72 hours. Thereafter, this bottom cell layer degenerates, rounded forms 
appearing and debris accumulating in the supernate. After 4 days, only 10-20 




















436 BIOCHEMISTRY; PUCK AND MARCUS Proc, N. A. 5S. 


per cent of the area of the Petri dish remains covered with recognizable cells, 
whereas, without irradiation, growth would be confluent. 

During the static period of the feeder layer, the test cells on the microscope slide 
multiply steadily, reaching a clone population averaging about 800-2,000 and a 
diameter of about 1.0 mm. in 8-10 days. At this point the microscope slide can be 
removed and the colonies fixed, stained, and counted directly. A typical set of 
such colonies is illustrated in Figure 2, and Table 2 presents a comparison between * 


, 





Fig, 2.—A—Typical clones produced from single HeLa cells growing on glass slides above a 
feeder layer (X1.7). The top slide had 5 times as great an inoculum as the bottom slide. B 
Photomicrograph of a typical clone (X50). 


TABLE 2 


COMPARISON BETWEEN THEORETICAL AND ACTUAL NUMBERS OF CLONES DEVELOP- 
ING ON Microscope Siipes INOCULATED WITH SEPARATED HeLa CELLS AND 
GROWN AS DESCRIBED 

Average No. of Single Cells 


Spread on Each Slide No. of Clones Found Mean and Standard Deviation 
{ 7,3,9,5,4,2 | * 
5 (91287,4 f 5.14 2.8 
(10,13, 10,20 | 
| 7,9,8,10,14 | 
10 { 7,5, 20 5 96+ 3.9 
| 8,9,7,11,6 | 
| 9, 11, 5, 6, 6, 10) 
{35, 44, 48 1 yee 
50 44, 43) 45, 44, 30/ i chaeceliad 
100 77+, 68+ 72.5+* 


* This figure needs correction upward, because the colonies on these slides were so crowded as 
to make it certain that an appreciable number had been lost by coalescence. 
the colony counts obtained and the average number of single cells contained in the 
aliquot deposited on the microscope slide. | Correspondence between these num- 
bers is as good as that obtained in bacterial assays. More than fifty such titrations 
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have been carried out, with no values falling outside the limits listed in the experi- 
ments of Table 2. 


Experiments in which an attempt was made to secure similar results without irradiation 
of the feeder-cell layer were unsuccessful because cells from the bottom layer detach and 
are carried up to the test layer, where they may reattach and initiate colony formation. 
Thus, when an empty slide is placed on top of an unirradiated cell layer, as many as 150 
colonies may develop which have originated from the bottom layer. Experiments which 
will be described later have demonstrated that the doses of X-irradiation employed inactivate 
the multiplication of 99.9 per cent of the feeder cells, so that, even if an occasional cell from 
the bottom finds its way to the test slide, it does not form a colony. 


Periodic counts of the cells in the developing clones permit estimation of the 
doubling time, which was found to be 20 hours in young clones. 

Preliminary results essentially similar were obtained with feeder cells produced 
by ultraviolet irradiation. However, unless quartz covers are available, the neces- 
sity for uncovering the Petri dishes during irradiation increases the risk of con- 
tamination. 

DISCUSSION 

The success of this titration permits the conclusion that almost all the HeLa cells 
obtained by trypsinization of a growing monolayer are capable of reproduction in 
an appropriate medium. The effectiveness of the feeder-cell principle also lends 
support to the idea that the sole reason for failure to obtain growth from single HeLa 
cells in the presence of large quantities of liquid is the loss of a short-lived, diffusible 
factor. New designs of simple growth arrangements for single cells, some of which 
eliminate the need for feeder cells, are being tested. 

Since clones of these cells can easily be produced which exceed the critical size 
necessary for transfer to tubes and eventually to large bottles, it is possible 
to produce the large clonal populations needed for many types of genetic studies. 
Experiments with this technique on genetic and physiological actions of several 
agents are in progress. The applicability of this procedure to other types of animal 
cells is also under study. 

SUMMARY 


A simple method for production of clones from almost all the cells of a HeLa 
population is described. The technique can be used to titrate the number of cells 
capable of reproduction in a suspension and should find application in genetic and 
metabolic studies involving the action on animal cells of viruses, drugs, nutritional 
stresses, and mutagenic agents. 


* Contribution No. 31. This study has been aided by a grant from the National Foundation 
for Infantile Paralysis and from the Damon Runyon Fund. 

1 J. Lederberg, Microbial Genetics (Madison: University of Wisconsin Press, 1951). 

2 W. R. Earle, J. C. Bryant, and E. L. Schilling, Ann. N.Y. Acad. Sci., 58, 1000, 1954. 

3G. O. Gey, W. D. Coffman, and M. T. Kubicek, Cancer Research, 12, 264, 1952. 

4 The composition of this solution, containing fewer constituents than Medium 199 (J. F. Mor- 
gan, H. J. Mortan, and R. C. Parker, Proc. Soc. Exptl. Biol. Med., 73, 1, 1950), was kindly sup- 
plied by Dr. Charity Weymouth, through Dr. Dulbecco, and will be published. 

5 Peyton Rous and F.S. Jones, J. Exptl. Med., 23, 549, 1916; R. Dulbecco and M. Vogt, J. Expil. 
Med., 99, 167, 1954; J.S. Youngner, Proc. Soc. Exptl. Biol. Med., 85, 202, 1954. 

6 In our earliest experiments test cells and feeder cells were placed in the same layer. We wish 
to thank Dr. Leo Szilard, who suggested the more advantageous geometrical arrangement in 
which the test cells were placed on top of the layer of feeder cells. 
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One of the geometrical isomers of vitamin A and retinene, arbitrarily called 
“neo-b,” is the precursor of the visual pigments rhodopsin and iodopsin.' 2 Vita- 
min A possesses five conjugated double bonds, to which retinene—its aldehyde— 


adds a conjugated carbonyl group (Fig. 1). 
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AS NEAINEA ENB 
H H H all-trans 
vitamin A 


c 
Ne 
S C No 


H 
retinene 
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Fic. 1.—Structural formulas of the unhindered geo- 
metrical isomers of vitamin A and retinene. The cis 
linkage occurs without serious steric hindrance only at 
double bonds 9 and 13, and only the forms shown are 
relatively coplanar. 


It was believed until recently 
that, in this class of com- 
pounds, only the double bonds 
from methylated carbon atoms 
—here numbered 9 and 13— 
could assume the cis configu- 
ration.* * At the other double 
bonds, the cis configuration 
encounters serious steric hin- 
drance, necessitating a large 
departure from coplanarity—a 
large twist in the molecule. 
It was thought that the conse- 
quent interference with reso- 
nance should make hindered cis 
forms highly unstable. Only 
four unhindered geometrical 
isomers of vitamin A or reti- 
nene were therefore expected: 
all-trans, 9-cis, 13-cis, and 9,13- 
dicis. 

Five crystalline geometrical 
isomers of retinene, however, 
have been described: all-trans, 
first prepared by Ball e¢ al. 
neoretinenes a and b, prepared 
in the Harvard laboratory ; and 
isoretinenes a and b, prepared 
in the Organic Research Labo- 
ratory of Distillation Products 
Industries in Rochester, New 
York (hereafter referred to as 
“EyPy”).** 
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The DPI workers assigned tentative structures to three of the cis isomers, which 
they have since confirmed.’ These are the unhindered forms, neo-a (13-cis); 
iso-a (9-c7s), first synthesized by Graham, Van Dorp, and Arens; and iso-b (9,13- 
dicis). Table 1 indicates the various numbering systems which have been attached 
to these compounds in the literature. 


TABLE | 


NOMENCLATURE OF GEOMETRIC ISOMERS OF VITAMIN A AND RETINENE* 


———____——— —NamME———— ——_—_—_— 
NUMBERING SYSTEM Neo-b Iso-a Neo-c Neo-a Iso-b 
Present 7-cis 9-cis 1 1-cis 13-cis 9,13-dicis 
Geneva system 8-cis 6-cis 4-cis 2-cis 2,6-dicis 
Hubbard and Wald 2-cis 3-cis 4-cis 5-18 3,4-dicis 


* The present paper follows Karrer’s convention for numbering the carotenoids, adopted by the 
Organic and Biological Nomenclature Commission of the International Union of Pure and Applied 
Chemistry (1947), following a recommendation of the Committee on Biological Nomenclature of 
the National Research Council (Chem. Eng. News, 24, 1236, 1946). It has the advantage of keep- 
ing the same numbers attached to corresponding carbon atoms in isoprenic polyenes of various 
sizes. Other workers have used the ‘‘Geneva’™’ system.’ 22. This proceeds from the end group 
and would assign different numbers to corresponding C eset in chains of various lengths. The 
arbitrary numbering system formerly used by Hubbard and Wald (J. Gen. Physiol., 36, 269, 
1952-1953) is also listed. 


The fifth isomer, neo-b, must therefore possess a hindered cis linkage. The 
existence of stable polyenes containing hindered cis configurations has recently 
been placed beyond doubt by the synthesis of a number of such molecules.?~'? 
Neoretinene b displays in its infrared spectrum evidence of a cis —CH==CH— 
linkage, necessarily hindered in this type of molecule.'* It is the first natural 
compound known to possess a hindered c7zs structure. 

Such a hindered cis linkage can occur only at double bond 7 or 11. The hin- 
drances in these positions are shown in Figure 2. Even ‘‘unhindered”’ cis linkages 


Cs 





A 


Fig. 2.—Types of steric hindrance encountered in cis isomers of 
vitamin A and retinene. Diagram A shows the relatively small 
hindrance between H atoms in such “‘unhindered”’ cis linkages as at 
double bonds 9 and 13. B shows the much larger hindrance between 
—H and —CH; groups involved in the 1l-cis linkage. C shows 
the still larger hindrance between methyl! piety By in the 
7-cis linkage. These diagrams are drawn as though the structures 
were coplanar. In reality, the superpositions shown are avoided by 
twisting of the molecule. 


involve some overlapping of —-H atoms as shown in Figure 2, A. The small twist 
in the molecule that this necessitates introduces a partial block to resonance. 
This is believed to be responsible for the fall in extinction and displacement of 
spectrum toward shorter wave lengths that ordinarily accompany an “unhindered”’ 
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cis linkage. Linkages of the 11-cis type (Fig. 2, B) are much more strongly hin- 
dered, through the overlapping of an —H with a—CH; group. Molecules contain- 
ing such linkages have been prepared synthetically.'? The 7-c7s linkage (Fig. 2, 
C), in which —-CH; groups overlap with one another, involves a degree of hindrance 
not heretofore encountered in polyenes. 

It has been suggested that neo-bcontains an unhindered in addition to the hindered 
cis linkage—indeed, that it is the 11,13-dicis isomer.'* The purpose of this paper 
is to show (a) that neo-b is a monocis isomer and hence must be either 7-cis or 11- 
cis; (b) through synthesis of the 11-czs isomer, that it is not 11-czs; and (c) that, 
for this reason and on the evidence of its own properties, neo-b is the 7-cis isomer. 
This structure raises an array of special physiological problems. 

I. GENERAL CONSIDERATIONS 

The absorption spectra of the unhindered isomers of retinene appear as solid 
lines in Figure 3. Zechmeister' has shown in plant carotenoids that a single czs 
linkage shifts the main absorption maximum in hexane 4-6 my toward shorter 
wave lengths and depresses its extinction. Two such linkages are additive in their 
effects." 
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Fic. 3.—Absorption spectra of crystalline geometrical isomers of retinene in ethyl alcohol 
solution. The four unhindered isomers are shown with solid lines, the hindered cis isomer neo- 
retinene b (7-cis) with a broken line. The spectra exhibit three general regions of absorption: 
the main band, appearing at 383 my in all-trans retinene and displaced 6-7 my toward shorter 
wave lengths for each cis linkage; the cis peak region, at about 285 mz; and subsidiary bands, 
at 255-260 mu. 


The unhindered retinene isomers follow this simple rule. The spectrum of the 
all-trans isomer lies at the longest wave length and possesses the highest specific 
extinction. The spectra of the 9-cis and 13-cis isomers lie 5—6 my lower in wave 
length (6-7 my in ethanol). The specific extinctions vary accordingly.® 
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The corresponding stereoisomers of vitamin A display similar relationships, 
with the exception of neo-a. In this isomer Ajax lies 3 mu higher in wave length 
than in all-trans vitamin A and appears, therefore, to violate the ordinary rules for 
cts isomers. ” 

As Figure 1 shows, however, neo-a vitamin A differs from the all-trans configu- 
ration only in the position of the —H and —CH,.OH groups about carbon atom 14. 
The length and shape of the conjugated system are not thereby altered in any way. 
There is no particular reason for such a terminal cis linkage to induce the kinds of 
effect ordinarily associated with interior cis linkages. Indeed, a shift of Amax to 
longer wave lengths, as in neo-a vitamin A, may be characteristic of terminal cis link- 
ages.!? The corresponding retinene, in which the 13-czs linkage is interior to the 
aldehyde group, displays the ordinary cis relationships. 

As shown in Figure 1, a czs linkage bends the polyene chain. ‘This creates a new 
electric vector at right angles to the main axis of the conjugated system. To the 
principal mode of electronic oscillation along the main axis, responsible for the main 
absorption band, there is added a smaller perpendicular component, responsible for 
a smaller absorption band at shorter wave lengths. Zechmeister has called the 
latter a “cis peak.”’ Its presence is a sign of bending, its height a measure of the 
degree of bending. A single, centrally located cis linkage causes most bending and 
hence the highest czs-peak. Two cis linkages oppose each other in this regard: 
by bending the molecule in opposite directions, they straighten it somewhat, re- 
sulting in a lowered cts peak. '* 

In our earlier discussions of retinene spectra, we assumed that the absorption 
bands in the region 255-260 my represent c7s peaks. Now that the configurations 
of these isomers are established, it is clear that the true ez7s peak region, in the sense 
intended by Zechmeister, lies at about 285 mu. On this basis, as expected, the 9- 
cis isomer, being most strongly bent, has the highest absorption in this region, the 
13-c7s less, the 9,13-dicis isomer still less, and the all-trans molecule least of all. 

Absorption maxima at shorter wave lengths than the cis peak will be referred to 
here as ‘subsidiary bands.”’ Such bands appear commonly among the carotenoids. 
In terms of modes of electronic oscillation, the main band can be considered the 
fundamental, the cis peak the first overtone, and the subsidiary bands further 
overtones.'* No simple correlations can yet be drawn between structural 
configuration and the subsidiary bands. 

Hindered cis linkages are associated with distinct changes in the absorption 
spectrum. The main absorption band loses what fine structure it may possess in 
the all-trans isomer; is displaced much farther toward shorter wave lengths than 
with an unhindered cis linkage—in the type of molecule here considered, about 
twice as far; and is depressed much further in extinction. The cis peak is normal, 
continuing simply to indicate the bending of the conjugated system; but strong 
additional absorption appears in the ultraviolet, which has been ascribed to partial 
chromophores induced by the break in coplanarity. !” 





Il. THE NEO-) ISOMER: MONOCIS OR DICIS? 


The absorption spectrum of neoretinene 6 is shown as a broken line in Figure 3. 
It exhibits typical hindered czs characteristics in the very low specific extinction of 
the main band and the extraordinarily high subsidiary band. The position of the 
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main band, however, is anomalous; Ajax iS the same as in the unhindered monocis 
isomers. We shall return to this point later. 

It has been suggested that neo-b is the 11,13-dicis isomer.'* This structure is 
incompatible with its absorption spectrum. The presence of the hindered 1 1-cis 
linkage alone should displace Ajax about 14 my from the all-trans position (see 
below). In neo-b vitamin A, a terminal 13-c7s linkage, as explained above, might 
pull Amax back slightly toward longer wave lengths. This factor, however, would 
not operate in the retinene, in which the 13-c7s linkage, being internal to the alde- 
hyde group, would be fully effective; and here the combination of the 11- and 13- 
cis linkages should result in a total displacement of about 20 my from the all-drans 
position in alcohol. Instead, one observes a total displacement of 7 mu, as with a 
single, interior, unhindered cis linkage. On this basis alone, a dicis structure 
seems out of the question. 

The isomerization experiments shown in Figure 4 lead to the same conclusion. 
On mild heating (3 hours in aqueous digitonin at 70° C.), neoretinene b isomerizes 





Effect of heat 


on retinene 
isomers 


o before incubation 
eafter incubation 


° 

© 
oe) 
j 


| | 
| T 


all-trans 


389 
| 





“4 


Extincti 

















| | i | | | i | | 1 | | 
T T T T T T T T T T T T 


560 580 400 420 350 370 390 410 350 370 390 410 


Wavelength ~ mp 


Fic. 4.—Incubation of retinene isomers in aqueous digitonin for 3 hours at 70° C. Neoretinene b 
isomerizes almost completely to all-trans. All-trans and iso-a retinene show no significant changes, 
neo-a and iso-b only slight deterioration. 
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almost completely to the all-trans configuration. All-irans and iso-a retinene re- 
main unchanged, and neo-a and iso-b deteriorate slightly (isomerization would 
cause the Ena, of these isomers to go up, not down as here). That is, the hindered 
cis linkage of neoretinene b isomerizes under these conditions, but not the 9-cts 
and 13-cis linkages of iso-a, iso-b, and neo-a. If neoretinene b contained any com- 
bination of a hindered with an unhindered cis linkage, the first product of its 
isomerization would be 9-cis or 13-cis retinene, which under these conditions would 
isomerize further only with difficulty. This is not the observed result. 

The irradiation of any retinene isomer with light which it can absorb leads 
eventually to the formation of the same quasi-equilibrium mixture of isomers. 
When this reaction is carried out in such a polar solvent as alcohol or water, the 
final mixture contains a large fraction of neoretinene b.! 

If, as suggested, neoretinene 6 contained a 13-cis in addition to the hindered cis 
linkage, it should be formed more readily from 13-cis (neo-a) than from all-trans 
retinene. If, on the contrary, it contains the hindered cis linkage alone, it should 
be formed more readily from the all-trans isomer. By taking advantage of the 
specific reaction of neoretinene b with opsin to form rhodopsin, one can measure 
the rate of formation of this isomer from any of the others. When neoretinene b 
is incubated in the dark with excess opsin, it forms rhodopsin with a rise of maximal 
extinction (/,4x) Of 1.5-1.7 times. This assay yields reliable results with as little 
as 10~‘ micromole of neoretinene b.® '° 

In Figure 5 we compare by this method the rates of formation of neoretinene b 
from all-trans and 13-cis retinene. Samples of both isomers in digitonin, containing 
about 0.050 micromole of retinene per milliliter, were exposed to the focused light 
of a 160-watt microscope lamp passing through a '/;) neutral filter. At intervals 
the absorption spectra were recorded, and aliquots of 0.014 micromole of retinene 
were removed and assayed for neoretinene b. The formation of neoretinene b from 
the all-érans isomer starts at a high rate and follows an approximately hyperbolic 
course. When 13-cis retinene is the starting material, however, there is an initial 
lag, and the course of synthesis is $-shaped. At the end of the reaction both 
samples contain about 0.0034 micromole of neoretinene b, about 25 per cent of the 
total retinene in the sample. An experiment performed at a higher light intensity 
yielded the same result in shorter time. 

It is apparent that neoretinene b is formed more readily from all-trans than from 
13-cis retinene. It is hardly possible, therefore, that this isomer contains a 13-cis 
linkage. On the contrary, the lag in neoretinene b synthesis when neo-a retinene 
is used as precursor suggests that it is necessary to isomerize the 13-cis linkage to 
trans before neoretinene 6 can be formed. 

In summary, therefore, neoretinene 6 is almost certainly a monocis isomer.  All- 
trans retinene appears to be its immediate precursor as well as the first product of 
its isomerization. 


il. 1l-cis VITAMIN A AND RETINENE (NEO-C) 

Since the neo-b isomer possesses a single cis linkage in a hindered position, it 
must be either 7-czs or 11-cis (Fig. 6). The synthesis of 11-cis vitamin A makes it 
possible to choose between these alternatives. It was prepared by the following 
procedure, the details of which will be published elsewhere :'7 
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SYNTHESIS OF |l-czs VITAMIN A 
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The crystalline phenylazobenzoate and 6-anthraquinone carboxylate were pre- 
pared by the procedures of Isler et al.'* The phenylazobenzoate has the following 
properties: m.p., 99-100° C.; Amax in ethanol, 317.5 my; molar extinction (e) = 
51,900. The 6-anthraquinone carboxylate has the following properties: m.p., 
139-140° C.; Amax in ether, 313.5 mu (e = 33,000) and 255.5 my (€ = 31,300). 

When mixed in chloroform solution with saturated antimony chloride, 11-c7zs 
vitamin A yields a blue product with Amax = 620 muy, identical with that formed 
by other geometrical isomers of vitamin A. 

The absorption spectrum of 11-c/s vitamin A in ethanol is shown in Figure 7. 
It displays the typical properties associated with a hindered cis isomer. Ajax in 
hexane is at 311 my (312 my in ethyl alcohol), 14 my below that of all-trans vitamin 
A. E (1 percent, 1 cm.) in ethanol is 920; « = 26,200. These are by far the lowest 
values of Ajax and specific extinction found among the monocis isomers of vitamin A. 

For comparison, the spectrum of neo-b vitamin A, obtained by reducing crystal- 
line neoretinene b!* with potassium borohydride, is also shown in Figure 7. It is 
markedly different from the spectrum of the 1 1-czs isomer. 

A portion of 11-cs vitamin A was oxidized to the corresponding retinene by filter- 
ing through a short column of dry manganese dioxide. The retinene was chroma- 
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Fic. 5.—Formation of neoretinene b by exposing all-trans and 
neo-a (13-c7s) retinene to light in digitonin solution. The rise of 
neoretinene ) was measured by incubating aliquots with excess opsin, 
so forming rhodopsin. All-trans retinene isomerizes to neo-b at a 
high initial rate, as though its immediate precursor. From 13-cis 
retinene, neo-b forms with a lag period, presumably owing to the 
necessity to form all-trans retinene first. 
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Fie. 6.—Structural formulas of the hindered 
monocis isomers of vitamin A. Both isomers are 
twisted as well as bent at the cis linkage, owing 
to steric hindrance, in the 11-czs isomer as shown 
in Fig. 2, B, and in the 7-cis isomer as shown in 
Fig. 2, C. 
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tographed on a column of alumina, weakened by mixing with 7-10 per cent its 
weight of water. 1l-cis retinene appears in the lowermost fraction. It was 
eluted and purified by two further chromatographings of the same type. We have 
not yet succeeded in crystallizing this isomer. 

1l-cis retinene is sterically unstable, much more so than 1l-cis vitamin A or 
than any other retinene isomer. When kept in solution even at —15° C. and dark, 
it isomerizes from day to day toward the all-trans condition. It is probable that 
all our preparations were to some degree isomerized. 
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Fig. 7.—Absorption spectra of the hindered monocis isomers of 
vitamin A in ethanol. The main bands lie at 311 my (neo-c) and 
321 my (neo-b). 


The absorption spectrum of 11-cis retinene, shown in Figure 8, again is typical of 
a hindered cis isomer. Ajnax in hexane is at about 355 my, 14 my from that of the 
all-trans isomer. The absorption band is broad and skewed by a sharply rising 
absorption below 310 mu, apparently associated with a large subsidiary band. EF 
(1 per cent, 1 em.), calculated from isomerization data and comparison with the 
antimony chloride product, is 620-700. As in the case of 11-cis vitamin A, these 
are the lowest \,,., and specific extinction among the retinene isomers. The spec- 
trum of crystalline neoretinene b, also shown in Figure 8, has very different prop- 
erties. 
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A final property is of special interest. A preparation of 11l-cis retinene was in- 
cubated in the dark for one hour with excess cattle opsin. Then hydroxylamine 
was added to break up fortuitous combinations of retinene with other molecules, 
and the product was exposed to light to bleach any photosensitive pigment that 
may have formed. A trace of bleaching occurred, maximal at 498 mz, apparently 





| | | | | 
104— Retinene in hexane a 
e neo-b 
- o neo-C i, 
0.84- a 
; a as 
Q 
x 0.64- a 
E 
5 q ms 
ae 
ty OF anne” a 
0.24 Se 
Or | | | | | | | ti 











Wavelength- my 


Fig. 8.—Absorption spectra of the hindered monocis isomers of 
retinene in hexane. The main bands lie at 355 mu (neo-c) and 
363 mu (neo-b). Neo-c retinene is highly unstable, and this prepa- 
ration may have been slightly isomerized at the time these measure- 
ments were made. 


due to rhodopsin. Its Hymax Was 12 per cent of that of the retinene employed. 
Neoretinene 6 under the same circumstances would have yielded rhodopsin of 
Emax 150-170 per cent of the retinene employed.'* The trace of rhodopsin ob- 
tained in this instance probably arose from a small contamination with neoretinene 
b; for in the alumina chromatogram neo-a, neo-b, and 11-cis retinene are adsorbed 
close together, in that order of descent. We conclude that 11-czs retinene has 
little or no capacity to form rhodopsin. 
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IV. DISCUSSION 
The properties of the two hindered monocis isomers of vitamin A and retinene— 
neo-b and 11-cis—are summarized in Table 2. It is plain that these are different 
substances. We shall refer hereafter to the 11-c7s isomer as ‘‘neo-c.’’ As the only 
apparent alternative to the 11-cis configuration, neo-b must be 7-cis. 


TABLE 2 
PROPERTIES OF THE HINDERED Monocis Isomprs or ViraMIN A AND RETINENE* 


—-——7-cis (NEO-b)- - -—-l1l-cis (NEO-c)- 
Vitamin A Retinene Vitamin A Retinene 


Amax (mu) in ethanol 320(1,200)  376(857)t 312 (920) 373 (700) 
Amax (mu) in hexane 318(1,170) 353(915) 310 355 
Rhodopsin formation ay Emax Vises 4 Trace? 
1.5-1.7 
times 
Biopotency (per cent) 23} 47t <38 
* The numbers in parentheses are values of E (1 per cent, 1 em.). Rhodopsin formation was 


tested by incubating retinene isomers one hour in the dark with excess opsin. Biological potency 
is stated as percentage effectiveness (molar) for rat growth and liver storage, compared with all- 
trans vitamin A acetate (= 100). 
Dieterle and Robeson, Science, 120, 219, 1954. 
t Ames et al., Abstracts, 126th Meeting of the American Chemical Society, New York, 1954, pp. 
36C (a) and 37C (6). 
§ Oroshnik e/ al., J. Am. Chem. Soc. (in preparation). 


This structure has not yet been proved unequivocally by synthesis. Attempts 
to make 7-cis vitamin A by hydrogenating the corresponding acetylene (7,8- 
dehydrovitamin A) by a variety of methods have so far failed. The 7,8 position is 
so extremely hindered (Fig. 2, C) that the addition of two hydrogen atoms in cis 
configuration at the triple bond would require a large rearrangement of the shape 
of the molecule, including the introduction of a large twist between the ring and 
side-chain. It is presumably this that makes hydrogenation so difficult. 

Spectrum of the Neo-b Isomer.—We have noted that, for a hindered cis form, 
Amax Of the neo-b isomer is displaced abnormally little from the all-trans position. 
The spectrum is shifted no more than those of unhindered monocis isomers— 
5 my in neo-b vitamin A and 7 muy in the retinene. 

The key to this anomaly goes back to a consideration raised by Oroshnik ef 
al.2 The absorption spectrum of all-trans vitamin A, though generally considered 
typical of a conjugated pentaene, is in fact highly degraded. Compared with a 
straight-chain pentaene with comparable substituents, the vitamin A spectrum has 
lost its fine structure, is displaced some 20 my toward shorter wave lengths, and has 
an abnormally low specific extinction. The reason for this spectral degradation is 
steric hindrance between the l-methyl groups and the 8-hydrogen atom, which 
forces the B-ionylidene ring out of coplanarity with the side-chain. Consequently, 
the double bond in the ring is only partially conjugated with the side-chain and 
makes a correspondingly small contribution to the chromophore. 

For this reason a hindered cis linkage in the 7,8 position has much less effect 
than one in the 11,12 position. It only decreases further the already reduced con- 
jugation between the side-chain and the ring double bond. The displacement of 
spectrum is correspondingly small, in the present instance as small as for an un- 
hindered cis linkage elsewhere in the chain. It is only in the 7,8 position that a 
hindered c7zs linkage could have this restricted effect. 
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Formation of the Neo-b Isomer.—The cis isomers of the polyenes represent higher 
energy states than the all-trans configuration, the hindered cis isomers very con- 
siderably so (cf. Pauling’). Hindered czs forms are therefore hardly to be expected 
under conditions of true thermodynamic equilibrium. It is in this sense that one 
may interpret Pauling’s conclusion that they should scarcely occur at all. They 
probably always require special conditions for their formation, including a special 
supply of energy. It is now clear, however, that, once formed, they are much more 
stable than had been anticipated. Indeed, on occasion they may be more stable 
than unhindered c7s isomers (cf. Oroshnik and Mebane!”). 

Two methods of forming hindered c7s isomers have been described: irradiation 
with light, used initially to form neoretinene b;! and hydrogenation of the corre- 
sponding acetylene.’ The latter process, as noted above, succeeded in the syn- 
thesis of neo-c vitamin A yet failed in that of the neo-b isomer. 

On the other hand, the simple exposure of any geometrical isomer of retinene in 
alcohol to light that it can absorb—blue, violet, or near-ultraviolet-—induces a large 
formation of neoretinene b. As Figure 5 shows, at the end of such irradiation about 
25 per cent of the retinene is in thisform. For so hindered a structure this seems an 
astonishingly high yield. It illustrates an important consideration. Among all 
the possible geometrical isomers of the polyenes, a few configurations tend to form 
selectively. The factors which govern this choice are still obscure but clearly 
override the distinction between hindered and unhindered forms (cf. Zechmeister 
and Pinckard''). The neo-b isomer appears to have such a preferred configuration. 

Another condition, important physiologically, aids in forming the neo-b isomer. 
An enzyme occurs in eye tissues which catalyzes specifically the interconversion of 
all-trans and neo-b retinene. In the dark its equilibrium lies far over toward the 
all-trans configuration; but in dim light a quasi-equilibrium mixture is produced, 
containing about 35 per cent neoretinene 6.” In this instance the geometry and 
energetics of the attachment of retinene to the enzyme protein appear to act 
selectively on the 7,8 linkage. 

Distribution of the Neo-b Isomer in the Tissues.—The only conditions at present 
known to form the neo-b isomer involve light. These are important in the eye 
and may operate in the skin, but of course are not available to internal tissues. 

Before its peculiar structure and the special problems which this raises were ap- 
preciated, the neo-b isomer was assumed to occur in liver oils and to circulate in 
the blood. In the first syntheses of “rhodopsin” from liver oil vitamin A, the 
fact was intentionally disregarded that it always possessed a lower A,,ax than 
rhodopsin extracted from retinas, in the belief that at the time it was more im- 
portant to learn the conditions of synthesis of this type of pigment than to hesitate 
over what seemed a minor disparity.’ ?! It is apparent now that these early syn- 
theses had yielded mainly or entirely iso-rhodopsin, so demonstrating the presence 
of iso-a vitamin A, not of the neo-b isomer. We are now undertaking new experi- 
ments to determine whether neo-b vitamin A occurs in the liver or in the blood. 
At present it is known to occur only in the eye. 

Nutritional Significance of Neo-b Vitamin A.—Beyond its special role in vision, 
vitamin A fulfills some general function in the tissues, the nature of which is not 
understood. We know only that warm-blooded vertebrates deprived of this vita- 
min stop growing, deteriorate, and eventually die. 
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Is neo-b the active form of vitamin A in tissues other than the eye? A first ap- 
proach to answering this question is offered by the recent bio-assays of .neo-b 
vitamin A and retinene performed by Ames ef al.22. Compared with all-trans vita- 
min A acetate, neo-a is about 75 per cent active, and iso-a, iso-b, and neo-b about 
23 per cent active, in growth and liver storage assays in the rat. On the same basis, 
all-trans and neo-a retinene have about 91 per cent, iso-a and iso-b about 18 per cent, 
and neo-b retinene about 47 per cent the molar activity of all-trans vitamin A ace- 
tate.*? An independent bio-assay of neo-c vitamin A shows it to have less than 3 
per cent the activity of the all-trans isomer.” 

According to these data, neo-b has a relatively low biological potency compared 
with all-trans. Indeed, these observations imply that the all-trans isomer may be 
the form used in general by the tissues and that other isomers may be effective in 
promoting growth and liver storage primarily to the degree that they serve as pre- 
cursors of all-trans vitamin A. In this sense the positions of the unhindered isomers 
seem clear, for the 13-czs linkage is known to isomerize to trans very readily, the 
9-cis linkage only with difficulty. The position of the neo-b isomer in this regard 
needs further clarification, both its low potency in the vitamin A, which in vitro 
isomerizes to all-trans very easily, and the doubling of potency in the retinene. 
The reason for the very low activity of neo-c vitamin A is still more obscure. Be- 
fore these bio-assays can be interpreted physiologically, two further problems need 
to be solved: (1) how efficiently the various isomers are absorbed from the in- 
testine and (2) the configurations in which they are deposited in the liver and other 
tissues. 

At present we can be certain of the occurrence of the neo-b isomer only in the eye 
tissues and of its function only as the precursor of the pigments of rod and cone 
vision. 

Vv. SUMMARY 

The neo-b isomer of vitamin A and retinene, a hindered cis form, is the precursor 
of the visual pigments rhodopsin and iodopsin. The present paper shows it to be 
a monocis isomer, therefore either 7-czs or 11-cis. The 11-cis isomer has been syn- 
thesized and found to possess very different properties from neo-b. It may be called 
‘neo-c.”” The neo-b isomer must therefore be 7-cis. Its properties are in specific 
agreement with this structure. The formation of the neo-b isomer raises special 
problems chemically and physiologically. Nothing is known of its biological occur- 
rence or function outside the eye. 


* This investigation was supported in part by grants to G. W. from the Rockefeller Foundation 
and the Office of Naval Research. We are indebted to the Organic Research Laboratory of Dis- 
tillation Products Industries for a gift to G. W. of crystalline neoretinene 6. 
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ISOTOPE SEPARATION BY IONIC EXPANSION IN A MAGNETIC FIELD 
By JOSEPH SLEPIAN 
WESTINGHOUSE RESEARCH LABORATORIES, EAST PITTSBURGH, PENNSYLVANIA 
Communicated May 12, 1955 


I. Tonic Gas, Single Ion, in Stationary (Hypothetical) Container.—Consider a gas 
at a very high temperature 7’, more than 100,000° absolute, or 10 electron volts or 
more, in a stationary (hypothetical) container, which reflects the ions and electrons 
which strike it completely and elastically. The gas might be a metal, with ions of 
mass m, and electrons of mass mp. Let the ion density be n per cubic centimeter 
and the electron density also be n per cubic centimeter. 

We will have 


! ‘ymyv,2 = i mou? = '/skT, (1) 


where v,2 and 4? are the mean-square velocity in any particular direction of ions 
and electrons, respectively, and k is Boltzmann’s constant. The ratio of m, to mo 
is 20,000 or more. Hence the ratio of Vu? to V 0,2 is several hundred (see Fig. 1). 

An ion impacting and rebounding off the wall conveys a momentum 2m, to the 
wall. An electron conveys a momentum 2mu to the wall. Thus the gas conveys 
a momentum nmv,2 + nmou? to each second, or a force of 2nk7 per square centi- 
meter of wall. 

Now let the walls slowly expand from volume V to V + dV cubie centimeters. 
The work done upon the total wall is 2nk7 dV. This must be equal to the dimi- 
nution of energy of the ions and electrons inside. 


—d(3nkTV) 2nkT aV, 
—3knV dT = 2nkt dV, 


dT ,, aV ,, an 
oe = t"/s—) (?) 
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II. Jonie Gas, Single lon, Moving in (Hypothetical) Container, with Reflecting 
Side Walls.—Let the ionic gas, with accompanying density of electrons n be stream- 
ing through the (hypothetical) container. Let the temperature of the gas entering 
be T)° absolute, the initial streaming velocity be vo cm/sec, and the entering sec- 
tion be So square centimeters. 

Let the gas be emitted from the container with section S cm.?; temperature 7’, 
absolute; and velocity v» em/sec (see Fig. 2). Instead of emitting the gas from the 
container, we can instead let the gas precipitate on real slats. By insulating the 
slats, we insure that the rate of arrival of electrons equals the rate of arrival of ions 
or that vp is the same for electrons as for ions. 


(3) 


OWG.57-0-5365 





Fic. 1.—loniec gas, stationary in (hypothetical) container with re- 
flecting walls. 
We have, for the random initial motion of ions and electrons, 
1/omy,2_ = '/omyu? = '/2kT', with V 42 several hundred times V 02. (4) 


For the random final motion of the ions and electrons, we have equation (1). 

There is no exchange of energy between the charged particles and the side 
walls, since the reflection of the ions and electrons is completely elastic. The total 
energy per ion is therefore constant: 


i oMyWoo" + 3kT a ome” ot 3kT, (5) 


neglecting the energy of the electron given by '/2movoo? and 1/:mov9”. Also, we have 
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NwWooSo = NUS. (6) 


Substituting equations (3) and (6) in equation (5), 


l 9 "WY 1 ™ pasar VooSo Vs rs 
2 MWoo" + 3kT> = 2 Myo + 3kT 5 g ; (4) 
Vos 


If '/.k7 is nearly equal to '/2myvoo, then, when So/S is nearly equal to zero, vp? = 
7v0”, and the temperature becomes 


me 1/3 cy2 y y 2/ f S ” 
1S"? T = §&” T, or T= 7(@) T°. (8) 


Thus the velocity of the gas approaches a limiting value 1) = +/7v00, and the tem- 
‘ yan 2 . eae te 
perature of the gas approaches zero as S~ “*, as S becomes infinite. 


OWG. 57-0-5366 
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Mo 


Fig. 2.—Ionic gas, moving in (hypothetical) container. Reflecting 
walls (hypothetical). 


Ill. Effect of a Magnetic Field Parallel to Side Walls, with Electric Potential 
Difference between Them; Other Walls of Insulating Segments; All Walls Real.—Now 
let us trace the effects of magnetic field and electric field indicated in the title of 
this section. The motion of the ions and electrons with mean radial velocity v 
in the magnetic field causes a mean force of —eHvy on the ions and a mean force 
of eHvy on the electrons. If the electric field produced by the potential difference 
on the walls is just enough to produce a mean force +eHv» on the ions and a mean 
force —eHv, on the electrons, then the electric field will balance the effect of the 
magnetic field on the uniform motion of ions and electrons. There, 
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E,=- ” = Hyp. (9) 
r 60 

What will now be the effect of V at smaller radii? A current will flow in the 
central part from the positive plate to the negative plate. At the positive plate, 
a space charge of negatively charged electrons forms. These charges are in the elec- 
tric field which is perpendicular to the positive plate. Under these conditions, 
the negative space charge moves outward parallel to the positive surface. 

In the same way, at the negative plate, a positive space charge exists, which is 
in the electric field perpendicular to the negative plate. The net effect is that the 
positive ions move outward parallel to the negative surface. 

It follows that everywhere where r is less than the value which makes He = Hv», 
the force is outward on electrons and ions, but it is not very much greater than 
normal because of the magnitude of the space charges and their velocities; there 
will be a voltage subtracted by space charge. 

If we now terminate the container by a set of insulated slats, the moving ionized 
gas will terminate by precipitating upon these slats, if only the electric potential 
of the side plates is high enough (see Fig. 3). 


OWG. 57-0-5367 
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Fig. 3.—Ionic gas, expanding in actual radial container. Gas is of 
two isotopes, m and m2. Curved arrowed lines are current lines of 
negligible magnitude; closed curve is of considerable magnitude. 


IV. Electric Field and Potentials Associated with Magnetic Field.—To calculate 
the electric fields and potentials, we need to know the value of vo, since, with Ee = 
Hv, the mean motion of the ions and the accompanying electrons is uniformly 
outward. Of course, the random superimposed motion of the ions and electrons of 
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temperature 7’ will give a current of both ions and electrons to the plates, but the 
currents will be equal to both plates. 

If the initial radial velocity is vo and the initial distance measured circumferen- 
tially between the side plates is Lo, then 


E, = How, Vo= - 10H vo0, (10) 


is the open-circuit voltage between the plates. 

Now we raise the voltage between the radial plates until it is more than enough 
to bring the discharge out to the insulated end plates. Now the field H, at the 
insulated plates will equal the negative current reaction, Hv. The potential of the 
last insulated plate will be — LE, greater than the potential of the first plate: 

E, = Hw; V = —LHw, (11) 

where L is the circumferential distance between the two side plates. 

Since, for L/L large, S/So = L/L, then 

VL 

— (12) 
J 0 Lo Vo0 


V. Gas Consists of Two Isotopic Ions: The Problem.—Let the initial composition 
of the gas now include isotopic ions of two masses, m; and m2; a radial velocity 
Yoo; and a temperature 7p, giving a random velocity 


. ‘my? a. 1/smov»? =! mu? = 1/2 kT. (13) 


Now it is obvious that we cannot have simultaneous solutions of the two equations 


Cy 2/3 
" 9 ‘ 7 Al f 9 € 7 VooSo , 
'/omywvoo? + 3kT,) =} omy” + 3kT Fox *) (7 ) 
vo S 
and 
VooS 2/3 
/ 917 . ‘ 000 whan 
. 2MWoo" + 3kT, = ! gM”? = 3kT ( =) “ ; 
Vo A 


for the same value, vp’ = vo”, for the two different mass isotopes, m; and mo. 
If the gas were all of isotope m, the velocity ve would take on a value vo’, which 
satisfies equation (7’) and 
E,’ = Hw’; V’ = —LAw’, (11’) 
the first equation holding at all the end plates and the last equation holding at the 
last end plate. 
If the gas were all of isotope me, the velocity vp would take on a value v9”, which 
satisfies equation (7”) and 
E,” = Hoo": ibe = — LH," " (11’’) 
for the field at all end plates and the potential at the last end plate. 
If m, < me, then vo’ > vo”, and EF,’ > E,”; —V’>-V". 
If ions of mass m; and m: are both present, then the 2, may be different for the 


different end plates, but we must have, for each one, 


E,’ > E,> E,’, -V’>-V>  —V,’. (11’”") 
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However, the temperature which the gas has at any end plate also would be 
different if the gas were made up entirely of isotope m; or of isotope m:. If the gas 
were all isotope m, the temperature at each end plate would be 


io - VooSo "a , 
T’ = r(“=>) (14’) 


If the ions were all of isotope ms, the temperature would be 


) Si Va 
7 = 7( ) (14’’) 


Vo" S 


The actual temperatures of the gas at the different end plates will all be different, 
but the temperatures 7’ will satisfy 
Pe FT", (14’"") 
VI. Gas Consists of Two Isotopic Ions: Problem Solved.—The problem set by 
the preceding section finds a ready solution if we remember that the assumption 
was made that the azimuthal composition of the ion stream remains constant for 
the one-ion-isotope case. For two isotopes, we do not need the azimuthal con- 
stancy of the ion stream, and it is readily seen that it is actually disturbed by the 
forces which set themselves up here. 
KE,’ and FE,” and V’ and V” are respectively different for the case of the two ions 
taken separately. Letting w; and w, be the velocity of the ions perpendicular to 
the ion stream, we have 


w, = 0 when ions of m only are present; (15’) 
W. = 0 when ions of ms. only are present. (15’") 


But each end plate can have but one electric potential, and it must take a value V, 
intermediate between V’ and V”, that is, 


-~V’> =-V > -9". eer) 


The lighter ion, m, must be retarded radially by the interactions of the electrons 
which are affected by the other ion, m:. The other, heavier ion, m2, must be 
speeded up radially by the interactions of the electrons which are affected by the 
lighter ion, m. 

The lighter ion, m, is then acted upon by an additional radial force coming from 
the additional interactions of electrons with ions m2. This force is a radial force 
and is directed radially inward. Corresponding to this force, the velocity w; is 
directed to the left in the figures and is toward the negatively charged side plate. 

The heavier ion, m2, however, is acted upon by an additional radial force coming 
from the additional interaction of electrons with ions m. This force is a radial 
force and is directed radially outward. Corresponding to this force, the velocity we. 
is directed to the right in the figures and is toward the positively charged side plate. 

The difference between w; and we will be given by 


H(w,; — w) = —V' + V" =A Liu’ — 00") (16) 


or 
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Ww, — We = L(v0’ ~ vo") = Livo we (17) 


cee | 


where m. = m, + dm. 

The velocity w; varies from 0 next to the left-hand, negatively charged plate to 
+Luo(dm/2m) next to the right-hand, positively charged plate; w. varies from 
—Lv(dm/2m) next to the left-hand, negatively charged plate to 0 next to the 
right-hand, positively charged plate. 

It is clear that these motions will make the end plate next to the negatively 
charged side plate have an excess of its share of the m, or lighter, ion, and will 
make the end plate next to the positively charged side plate have an excess of its 
share of the me, or heavier, ion. 

It is also clear that there will be an opposing flow of ions due to the effect of the 
temperature upon the ions of the two masses. The excess of ions of mass m, in the 
end plates to the left will cause a thermal motion of ions of mass m, to the right. 
The excess of ions of mass m, in the end plates to the right will cause a thermal mo- 
tion of ions of mass m, to the left. Fortunately, the temperature 7 falls off as 
(vo0So/voS), as S increases. 

Thus it is clear that the enrichment of the ions among the end plates varies per 
unit length nearly as the first power of the total length and that the enrichment over 
the whole length of end plates varies nearly as the square of the whole length of end 
plates. 

VII. Initial Conditions.—A convenient source of ions and electrons is an electric 
are of short length, '/s inch, parallel to the magnetic field. The negative electrode 
is a cored carbon such as is manufactured by the National Carbon Company. It 
is hardly consumed at all. The positive electrode is of the metal whose isotopes 
are wanted. Some of the time we used solid uranium. Sometimes we used a mix- 
ture of uranium with uranium oxide. 

With an are of 10 amp., we obtained a positive-ion current of 500 ma. to 1 amp. 
We have tested an are of 100 amp., obtaining a positive-ion current of 10amp. No 
limitations of vacuum were encountered. With the metal anodes, the neutrals were 
much less than 10 per cent of the number of ions. 

The radially disposed side plates should be placed as close as possible to the are. 
We set the beginning of the side plates one inch from the are. We adjusted the 
magnetic field so that when the upper and lower one-inch circles around the are 
electrodes were excited by more than 20 volts, the radially disposed side plates 
with zero potential difference between them floated at 20 volts negative to the 
are anode. About 6,000 gauss sufficed. The enrichment with —7 volts potential 
difference between the two radial plates, was with uranium, about one-half of 
1 per cent. 

The are voltage was 20 to less than 100 volts. The voltage between the radial 
side plates would take negligible power. We obtained grams of uranium deposit 
per kilowatt of arc energy. The ion random energy was about 5 volts in each direc- 
tion. '/omv,2 = '/2k7. The mean ion energy was also about 5 volts. We ob- 
tained an ion current of about half an ampere. The vacuum vessel was about 
three feet high and four feet in diameter. However, we used about one-foot- 
diameter central portion. 








GRAVITY ANOMALIES AT CONTINENTAL MARGINS 
By J. LAMAR WorZEL AND G. LYNN SHURBET 


LAMONT GEOLOGICAL OBSERVATORY, * COLUMBIA UNIVERSITY, PALISADES, NEW YORK 


Communicated by M. Ewing, April 27, 1955 
Introduction.—The oceans are the normal part of the earth’s crust and the con- 


tinents the anomalous part, 64 per cent of the earth’s surface being covered by 
oceans deeper than 1,000 fathoms. 


The M-discontinuity between the crust and 


the mantle lies 9-12 km. below sea level in ocean areas and 30—40 km. below se: 
level beneath the continents. 


At the continental margins the crustal thickness 
changes by a factor of about 5. There are few data on this transition. It is 
difficult to investigate by explosion seismology because of the thickness of the 
layers, the complex nature of the interfaces, layer slopes, changes of lithology, etc. 
Gravity surveys at present give the only data which cover this transition zone in 
detail. 


Observations.—Gravity observations at 104 stations were made on board U.S.S 
“Tusk” at the locations shown in Figure | 
Meinesz 


pendulum apparatus. ! 
with an 


The observations were made with a 
The data for Browne corrections were obtained 
auxiliary long-period pendulum apparatus? loaned by Professor Vening 
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Meinesz. Modifications in the apparatus were made and are described by Worzel 
and Ewing.’ 

Seven gravity sections are shown in Figures 2-8, the distance from a point near 
the shore line being given in nautical miles and kilometers. Part A of each figure 
shows the topographic section and the seismic data available. Part B shows the 
assumed structures and densities used for the gravity calculations. Part C shows a 
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PORTLAND SECTION, USS TUSK (SS 426) 
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WOODS HOLE SECTION, USS TUSK (SS 426) 
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CAPE HENRY SECTION, USS TUSK (SS 426) 
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CAPE HATTERAS SECTION, USS TUSK (SS426) 
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comparison of the anomaly curve computed from Part B and the observed free-air 
gravity anomalies. Part D shows the generalized structure deduced from the 
gravity and seismic evidence. 

The gravity calculations were made by two-dimensional analysis of a section 
approximated by rectangular blocks. The values used for specific gravity were 1.03 
for sea water, 3.27 for the mantle, and 2.30 for the sediments. The specific gravity 
2.84 was chosen for the oceanic and continental crustal rocks. Together with the 
other densities and thicknesses chosen, this value indicates isostatic equilibrium be- 
tween the mean continental column and the mean oceanic columns.* Where seismic 
depth determinations were available, they were adopted; otherwise, the depths of 
interfaces were adjusted until the gravity data could be adequately fitted. In 
sections where the seismic data were inadequate, guidance was obtained from near- 
by sections. ; 

On several sections it was desirable to use land gravity observations to complete 
the gravity sections. Bouguer anomalies were used for the land observations and 
free-air anomalies for the sea observations. Although named differently, these 
anomalies are strictly comparable. The Bouguer anomaly on land depends on 
deviations of density from the standard continental column, and the free-air 
anomaly at sea depends on deviations of actual density from densities which would 
put the column in isostatic balance with the standard continental column. An 
additional contribution to the free-air anomaly arises from “edge effects.”’ 

Figure 2 shows the Mount Desert section. Only two seismic stations are avail- 
able.» © There is little control for the sedimentary thickness beneath the con- 
tinental rise and slope. However, in view of the other sections to be discussed later 
and unless improbably narrow and sharp fluctuations of the M-discontinuity are 
introduced, the sedimentary section deduced here must be approximately correct. 
It is obvious that the transition from the crust beneath the continents to the crust 
beneath the ocean occurs within about 200 km., starting approximately at the 
northern boundary of Georges Bank. There is structure on the basement surface 
at the continental margin beneath approximately the 1,000-fathom curve which 
is required in order to fit the gravity data, unless the line of section follows a sub- 
marine canyon in this vicinity. 

Figure 3 shows the Portland section. The seismic data’ for the Gulf of Maine in 
this section show less than 0.3 km. of sediments and water, which was too thin to 
show on the scale of our drawing. The computed gravity curve fits well with the 
gravity data. Note the dip in the curve at approximately 175 km. from Cape 
Elizabeth, which corresponds to the deeper water in the central Gulf of Maine. 
One can choose a sedimentary wedge thickening seaward from about the northern 
edge of Georges Banks, with the greatest thickness occurring near the southern edge 
of the section approximately 550 km. south of Cape Elizabeth. The transition from 
the continental crustal thicknesses to the oceanic thicknesses is somewhat steeper 
here than in the previous section. Again there is structure beneath the continental 
slope. Some fluctuation of the M-discontinuity must be allowed in order to fit 
the gravity data on the inner part of the Gulf of Maine. Any gradual change in 
the thickness of the sedimentary wedge could be compensated in gravity effect by a 
gradual shift in the M-discontinuity. The section chosen conforms with neighbor- 
ing sections. The Gulf of Maine is a continental block that is flooded. 
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Figure 4 shows the Woods Hole section. Here we have seismic measurements’ 
of the sedimentary thickness to the basement surface from shore to about 120 km. 
beyond the continental slope. Since we cannot move these boundaries, most of 
the adjustment for gravity calculations had to be done on the M-discontinuity. 
The fit of the computed curve to the observed points is satisfactory. G.P. Woollard 
(personal communication) provided the gravity values for the land part of the 
section. The northern end of the profile crosses the eastern end of a rather large 
gravity feature. We did not attempt to fit these data, as we had too little evidence 
of its structure within our section to do it justice. The transition of the M-dis- 
continuity from land to sea is more rapid here than in the previous sections (which 
had little seismic control), although tapering off on the seaward end more gradually. 
The structure on the basement surfaces which we had to add in the previous sections 
is shown in the seismic data. The seismic data confirm that the greatest sedi- 
mentary thickness occurs on the oceanic side of the continental slope. This sedi- 
mentary wedge obviously thickens to about 350 km. from Gay Head and then thins 
seaward to a seismic station 100 km. beyond the end of the section. The bulk of 
the sediment lies beneath the continental rise. 

Figure 5 shows the New York section. There is good seismic data’ for the base- 
ment surface beneath the continental shelf. Beyond the continental slope there is one 
section at a distance of about 280 km. from Fire Island Light and a second seismic 
section 530 km. southeast of Fire Island Light just off the edge of the diagram. 
The fit of the computed curve with the observed data is quite satisfactory with the 
sections shown in section B. It was necessary to put some structure on the M- 
discontinuity surface beneath the continental part of the block in order to fit the 
curve adequately. Two great accumulations of sediment are found. The thick- 
ness is about 5 km. on the shelf and about 6 km. at 300 km. southeast of Fire Island 
Light. The rise of the M-discontinuity from typical depths starts beneath the 
first flexure of the basement, and typical oceanic crust is found at a distance of 
about 200 km. The greatest sedimentary thickness beneath the continental rise is 
found where the Hudson Canyon delta was reported by Ericson et al.8 

Figure 6 shows the Cape May section. There is considerable seismic detail” ° 
beneath the continental shelf, at the base of the continental slope, and seaward at 
about 420 km. The fit of the computed gravity curve with the observed data is 
quite satisfactory. There is less structure on the basement surface than for pre- 
vious sections. The sediments thicken gradually seaward, achieve their greatest 
thickness beneath the continental slope, and then thin gradually seaward. A large 
part of this section cuts across the Hudson Canyon delta on the continental rise. 
The M-discontinuity rises more gradually beneath this section than the previous 
ones, starting its rise approximately at the shore line and achieving its typically 
oceanic depth at about 250 km. from the shore line. This is the greatest accumula- 
tion of sediments indicated on any of the sections. 

Figure 7 shows the Cape Henry section. The seismic data” from the fall line to 
the 100-fathom curve is shown at the top. An additional seismic section 400 km. 
east of Cape Henry Light is available. The computed gravity curve fits quite 
satisfactorily with the observed values. Note the steepening of the gravity curve 
at about the 1,000-fathom curve. No gravity observations could be made in the 
interval from 0 to 90 km. from Cape Henry because of the water depth and the 
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very busy ship channel. The sediment thickens seaward again, reaching its great- 
est thickness at the foot of the continental slope or the top of the continental rise 
and then thins gradually across the continental rise. The M-discontinuity 
rises quite rapidly in this area from approximately 50 km. from Cape Henry Light 
to a normal sea depth at about 200 km. from Cape Henry Light. 

Figure 8 shows the Cape Hatteras section. The data inshore is made available 
by Skeels.!!. There is no additional data until just beyond the end of this section 
at 520 km. from Cape Hatteras Light. The fit of the computed curve with the ob- 
served anomalies is quite satisfactory. Note how extremely steep the anomaly 
curve is near the 1,000-fathom curve. The region from 30 km. west to 30 km. east 
of Cape Hatteras Light was not observed because the water was too shallow. The 
sediments thicken seaward, achieving their maximum thickness at the base of the 
continental slope and then thins gradually seaward. Again the maximum 
thickness of sediments occurs beneath the continental rise. The M-discontinuity 
has the sharpest rise here of any of the sections, rising from its continental depth to 
a typically oceanic depth in a distance of only about 80 km. 

Discussion and Conclusions.—We conclude that the true edge of this continent 
occurs at about the 1,000-fathom curve. The maximum sedimentary thickness is 
found near the base of the continental slope, and there is a significant amount of 
sediment across the whole continental rise. By far the greatest volume of sedi- 
ments is found on the ocean side of the continental slope, probably accounting for 
the existence of the continental rise. This distribution of sediments probably re- 
sults from the action of turbidity currents. 

If all the sediment were removed and these areas remained closely in isostatic 
balance, as they are at present, there would be a continental shelf floored with 
basement rocks, a continental slope floored with basement rocks, no continental 
rise, and a deep oceanic area floored by nearly level basement rocks. This must rep- 
resent an earlier stage in the development of this coast. 

Sedimentation, mostly from the continent, must have produced the present 
structure. As sedimentation continues, the oceanic crust is depressed closely in 
isostatic equilibrium until the sediment surface reaches sea level. The continental 
slope moves seaward over this thickened sedimentary section as the top of the sedi- 
ment approaches sea level. All these sections represent intermediate stages of this 
process, with the New York and Cape May sections the farthest advanced, probably 
owing to the much larger supply of sediments to this area. The Gulf Coast geo- 
syncline represents a much more advanced stage of this process.‘ Of course, if 
orogeny occurs within the region, this process is interrupted. 

The steepness of the M-discontinuity varies considerably and is not simply re- 
lated to the near-surface structure at the continental margin. The continental 
crust thins fairly abruptly in about 200 km. to the oceanic crustal thicknesses. 
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THE SEISMIC SURFACE PULSE 
By C. L. PekErRIS 
DEPARTMENT OF APPLIED MATHEMATICS, WEIZMANN INSTITUTE, REHOVOT, ISRAEL 


1. Introduction.—The problem under investigation is to determine the motion 
of the surface of a uniform elastic half-space produced by the application at the 
surface of a point pressure pulse varying with time like the Heaviside unit function. 
The original formulation of the problem is due to Lamb,' who synthesized the 
solution for the pulse from the periodic solution. Lamb’s method is, however, 
very intricate. In a previous publication? the author gave an exact and closed 
expression for the vertical component of displacement for the case when the pressure 
pulse varies like the Heaviside unit function H(t). The derivation of this result, 
as well as the solution for the horizontal displacement, are given in this paper. 
The seismic pulse problem was treated nearly simultaneously by Cagniard,*® and 
more recently by Pinney‘ and Dix. Because of the complexity of the analysis, 
it was thought worth while to reproduce in this and a subsequent publication the 
original solution for the surface source and the buried source. 

2. Formal Solution.—In this section we derive a formal solution for the problem 
of the motion produced by a seismic source buried below the surface in a uniform 
elastic half-space, when the time variation of the pulse is H(t). The solution for 
the surface source will then be obtained by letting the depth of source H approach 
zero. Referring to Figure 1, we choose a cylindrical system of co-ordinates with 
origin at the level of the source and the surface situated at z = —H. Quantities 
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referring to the portions of the medium below and above the source will be desig- 
nated by the subscripts 1 and 2, respectively. 

The space variation of the source will be specified by the condition that at the 
level of the source (2 = 0) the surface integral of the applied stress (p,.1 — D222) 
shall be equal to Z (negative): 


2fo° (pe — Px) dr = Z. (1) 
One representation of this point source is 


(Dez ~~ P22) = lim fle, r), 
e—>0 


4 


: Z - ee Ze * a, — 8/2 
fle, r) = { Jo(éerje * EdE = (e+ rr)” (2) 
Qa 0 2r 


The vertical displacement 


Z=-H SURFACE w and the horizontal displace- 


ment q can be derived (Lamb,! 


@ p. 29) from 





qd = Or + Xrzy 
Z=O—— * SOURCE me + kn ~ bx (3) 
| oO where the subscripts denote 


partial differentiation, and the 
potentials @ and x satisfy the 
wave equations for dilatational 


Fic. 1.—A point source is situated ata depth H below the and equivoluminal motion, re- 
surface of a uniform elastic half-space. 





spectively : 
V*o — h*o = 0, V*x — k*x = 0, (4) 
a . A+2 
h? = s k? = oe gS, c,* = —— = 3c’. (5) 
C5" ce p p 


Here c, denotes the velocity of compressional waves, c the velocity of shear waves, 
and p denotes the operator 0/dt. A and u are the elastic constants of the medium, 
which in the sequel we shall assume to be equal. 

The normal stress p,, and the shear stress p,, are given by 


Pz = Ah? + 2u(o.. + Xzzz k*x,), (6) 
) 7 € 9 Lod 
Prez = BD (2¢. + 2X22 bate kx). (7) 
or 
It is understood that after we have obtained the operational representation of 


w(r, 2, p), say, the actual w(r, 2, t) will be obtained by performing the integration 
over the Bromwich contour 


l peti (em 
w(r, z,t) = 5 J ” w(r, 2, p) dp. (8) 
“7 a-ia2 
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Appropriate solutions of equations (4) in regions | and 2 are 
d: = Ae~**Jo(tr), x, = Be~™*Jy(ér), (9) 
go = (Ce~** + De'™)Jo(ér), x, = (Ee ** 4 Fe**) Jo(ér), (10) 
where 
ka = (+h)? ke = (8 +k)", (11) 


and the functions are to be integrated over a suitable path in the &plane. This 
path, as well as the functions (of € and p) A, B, C, D, E, and F, are determined 
from the following six boundary conditions. At the level z = 0 the required con- 
tinuity of g, w, and p,, yields 


(6+ x21 = (O+ xz)2, (12) 
(b. + X02 — h*x)a = (be + Xe — k?x)2, (13) 
(26, + 2xs — k?x)1 = (2¢, + 2x5 = k?x)s, (14) 

while the source condition (1) is met, in the limit of vanishing e, by putting 
Z “ 
(Dzz1 — Da) = 9 Jo(Eryé (15) 

2a 

and integrating with respect to — from zero to infinity. At the surface z = —H, 


Pz2 and p,,, aS given by equations (6) and (7), must vanish. Solving these six 
boundary equations for A, B, C, D, E, and F, we get 


Zé 


7 = — D= —-keF, B=E+F, A=C-D, 
12mukBh? D, (16) 

kBF 
etc = { [(2e2 + k?)? + 4k*t%aBle**" — 4¢2(2¢2 + k2)e~*P#} (17) 


> 
itek ; {4k2(28? + k2)aBe~**# — [(2E% + k*)? + 4h%%aB]e~*#}, (18) 


M==[(2e? + k?)? — 4k?Ea8]. (19) 


Substituting these values into formulas (9) and (10), we can derive from equations 
(3) the expressions for the operational representation of the displacements every- 
where in the medium. The values of the latter at the surface are 


Zk 
OO te J, Jo(er)é[— (28? + ke *# 4 2e%——#OM (i) dé, 
2mu Jo M 
(20) 


s f- ! 
q(p) = dr, 9, p) = —- va» Ji(ér)§2[—2k2 ape *e# 4 (2g? 4 k2)e*OH] x 


1 
(*) dé. (21) 
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3. Vertical Displacement.—The expression for the operational representation of 
the vertical displacement in the case of a surface source can be obtained from 
formula (20) by letting the depth of source H approach zero: 


Zk ) 
w(p) = _ i) Jo(eryéa[(2E? + hk)? — 4k%%a8]—! dé. (22) 
0 


By writing & = kx, this expression simplifies to 


Zk 
w(p) = (= \ween, (23) 
N = fo Jo(krx)xm(x) ax, (24) 
where 
Qa 

mr) = ‘ e : ; 

[(2x? + 1)? — 42%a8) ] 
lates 2 1/ 2 rats t Pp Or 
QM £47 y. B= Vz? +1, i=". (25) 


and k appears only in the external factor and in Jo(kra). 

The interpretation of the external factor is simply 0/dt. In order to interpret 
the integral, we transform it into the form [ exp [—&Af(r, x) ]g(x) dz, with f(r, x) 
real. The interpretation of the latter is fg(x)H[t — (1/e)f(r, x)] dx, H denoting 
the Heaviside unit function. This can be accomplished by transforming the path 
of integration in equation (24) into one along the imaginary axis v in the complex 
z-plane, when Jo(krx) goes over into Ko(krv), the interpretation of which follows 
from the representation 


Ky(krv) = f° exp (—krv cosh 6) dé. (26) 

Now m(x) has branch points at x = + (i/V 3), +7, and poles at 
2 = & sy, vy = ¥V3 + V3. (27) 
We shall cut up the z-plane in such a manner as to Jeave m(x) regular to the right 


of the imaginary axis, in order to be able to apply to the integral in ecuation (24) 
the theorem of Bateman and Pekeris,® leading to 


- 2 - 
N= i) Jo(kra)am(2) dx = — ( )if Ko(krv)vm(iv) dv, (28) 
0 Tv 0 


I denoting the imaginary part. This is shown in Figure 2, where the values of 
a and 8 in the various portions of the imaginary axis are also indicated. Near 
the pole 7y the path is indented by a semicircle passed in the positive direction. 
It is easy to see that the integral over the semicircle is real and therefore does not 
contribute to equation (28). We must, however, take the principal] value of the 
integral in equation (28) once we pass the pole. Since, moreover, Jm(iwv) = 0 for 
ef 1, V3, we obtain 
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Tv [4/3 


9 xo 
N=- (*) p { _ Ko(krv)vm(iv) de, (29) 
1/ 


where P denotes the principal value. 












































tee oe 4 ] 
4 
| | 
B a V | | 
ie iV S 3V3"V3 | | 
Viv iW SF | | 
| | IN3 
| | Vi-v¥ “y | | 
fr | | X+ PLANE , X=U+iV 
| | bel os 5 U-—- lo 
| Irv’ a 
| | —-—— 
| V3 








| | 
| 7 nen eee Se 
“VA -iVEq EVENS 7 


| a) | i 


Fia. 2.—The complex z-plane is cut up by branch lines (thick) at 
the branch points v = +(1/ V3) andv = +1.0. The poles are situ- 


ated atv = #3/2/3 + /3. a and 6 denote the values of Vil + x? 


and of 1 + 2? on the various portions of the v-axis. 











Now the interpretation of Ao(krv) is 
n “ —rpv cosh 6 
Ko(krv) = J exp ( ) dé —> 0,v> 1, 
0 c 


— cosh! (") v<rt, (30) 
» 





where + = (ct/r), and — denotes operational interpretation. Hence 
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1 fo l 
kKo(krv) > ( ) cosh (*) = — (¢*? — p3)~ ”, o< fF, 
c \ot v r 


Q, v>r, (31) 


| 
kN — 0. rT < = 
V3 


2 | | 
kN—> — ( ) P| vm(iv)(r2 — v2)~'” do, tT> 5. “CZ) 
Tr 1/V3 V3 


Using the values of a, 8 shown in Figure 2 to evaluate Jm(iv), and substituting in 
equation (23), we get 


] 
T< , 
V3 
w(r) = 
ul V3 
Z 


Z ] 
extn), <r<il, (8) 
ur 


w( 7) 


JIG + G2(r)], 7>1, (34) 


“(2 
(2 


where 


T Vy? — 1/1 — Qy2)? dy 
Gr) = P f SR ns allel 3. (35) 
V3 Vr? — v°(3 — 24v? + 5604 — 320) 


viVv? — 1 [40 — 4/3] dv 


Gxt) = P f Ps — (36) 
1 Vr? — vX(3 — 2402 + 56v4 — 32v4) 


The integrals G,(r) and G2(r) can be evaluated in closed form by partial fraction 
decomposition. Writing, in equation (35), 
vp? = 1/3; + w* sin? 6; / (37) 


the integral is transformed into 


l x/2 
Gi(r) = 96 pf dé 
‘ 0 


a 12 + a eed cis) Seecenaete in cen ES re = 3 \ (38) 
\ 7 ['/io + w? sin? 6] (—b + w? sin? 6) (¢ + w? sin? 6))- 
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B se 5 / : Cue =, 
V3 


/ 2. (39) 
Using the result 


he ase T is He 
0” d0(a? +? sin? @)-! = — (a? + w)~”, 
<a 
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P So’? do(—8? + w? sin? 6)-' = 0, B<w 
us 9 9\ —'/2 
= Ore”, ta 


we obtain 


Chr) oa 

a V3 ‘ V3V3 +5 N3V3-—5 
a P — oo T Y; 
Vere NG + V3) — 22 Vr — 148 — V3) 


/ 


- V3 V3V3-5 | 


G(r) = —(-6 + ~ — = 2 t>vy. (42) 
96 | Vrr—'/, Vr? - 1/,(3 — V3)\ 
Similarly, using the substitution 
v2 = | + sin’ 8, o = 7? — ]l, (43) 


in equation (36), we obtain 
Go(r) = '/4 fo’? do[—3 — 3-(3 + 4a? sin? 6)-' — (1 + V3). 
(1 — V3 + 4a? sin? @)-! — (1 — V3)(1 + V3 + 4a? sin? 0)-*] 
x | V3 V3V3 45 
= —< : -+ - 
96 } Vr-—', V'/,(3 + V3) — 7 


\3V3—5 i 
<P, 9 TSN 
Vr? — 1/,(3 — V3) 


/ / 3 
V3 V3V3 —-: 
~£)-4+—2* + a r>y. (44) 
96 | it Gaal 2, Ve — 143 — V3) 
It follows that 
1 
w(r) = 0 <-— =, (45) 
(tr T V3 ( 





V3 V3V3 + 5 V3V3 - 5 
Z 6 — eo _— oe + — — ’ 
Ws) oe. Fe ee ee 
Jet R92 cere 


—=-<1r<1, (46) 
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sy V3 5 Saeapee 
gee N3Vv3 + : l<r<y ='V34+ V3, (47) 


w(r) = — ——)J6 - 
16aur [3 V3 
woe 
3Z 
w(r) = — t> vy. (48) 
Saur 


4. The Horizontal Displacement.—The operational representation q(p) for the 
horizontal displacement at the surface in case of a surface source can be obtained 
from formula (21) by putting H = 0: 


Z © Ji(ér)&?[(2E? + hk?) — 2k? aB] dé 
gp) = — = ae age epee (49) 
2ru Jo [(2&? + k*)? — 4k*?a8] 
which simplifies again, by the substitution ¢ = kz, to 
; Z ) 
) = ; 50) 
uP 2mu (3 “ 
Q = f° Jo(krx)n(x)zx dz, (51) 
| (22? + 1) — 2V 22°41 V2? + 1), f 
n(x) = qiesbaey (52) 


(Qa? + 1)? — 4r?°V 2? + 1 Ve + '/5 


By cutting up the complex z-plane in the manner shown in Figure 2, we may 
again apply the Bateman-Pekeris theorem to equation (51), leading to 


@ y) eel 
Q= f Jo(kraz)n(x)a dx = — (*) rf Ko(krv)n(w)v dv — '/4Ko(ky), (53) 
0 T 0 


where the last term arises from the residues at the poles x = +iy. Now 


' I 
In(w) = 0, <-— 
V3 
rs 9 f / 2 ee . 
eee ere 6) l <o<i, (54) 
(3 — 24v? + 56v4 — 32r°) ° 3 v/3 
= 0 v>l 
so that 
9 Lf Ko(krv) Vv? — 1/3 V1 — v°(12 — 2402)v dv Ky(ky) eas 
= - - — — - ee pi ee sag: =) ~ mae — - 5 Y5) 
TJ1/V3 (3 — 24? + 5604 — 320°) 4 


In order to obtain the interpretation of 0Q/dr, we differentiate equation (30) 
with respect to r and get 


or 
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y rV 72 — v? 


v<t, (56) 
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which, when substituted into equation (50), yields 


1 
q(r) = 0, T< Va’ 
Z l 
g(r) = — ontar tR,(r), Va << I, 
Z 
g(r) = — (, . ) pa), l<r<y, (57) 
2r7ur 
Z Z 9 9 —/3 
qQir7) = — tRo(r) + (7? — y*) i > +4; 
2r7ur Sarur 
where 
5 oVy? — l/s tf dnive v?(12 — 24v?) dv 
Ri(r) = f eer ee (58) 
V3 V 72 — vX(3 — 240? + 564 — 3208) 
1 yy? — 1/, V1 — v2(12 — 240%) dv 
Ry(s) « { NS) eee (59) 
1/V3 Vr? — v2(3 — 24v? + 56v* — 320) 


By partial fraction decomposition, Ri(7) and R2(r) can be expressed in terms of 
elliptic integrals, as follows: Using the substitution (37) in equation (58), Ri(r) 
is transformed into 


ae 2 Se [3 ap 8 2 Mi ani 
-™ 8 V2Jo Vi—ksin?0]- (1 + 8k? sin? @) 
6 —4V3) (6 + 4V3)_ 


[1 — (12V3 — 20)k? sin? 6] [1 + (12V3 + 20)k? sin? 4] 








[3 y 
- {6K(k) — 1811(8k2, k) + (6 — 4V3)N[—(12V3 — 20)k2, k] + 
(6 + 4V3)M[(12V3 + 20)k2, k)}, (60) 
where 

k? = 1/,(3r% — 1), (61) 

/2 dé /2 d6 
K(k) = sheneies- Teo Oe, a f SDNY nec troy eT: ar POTy 
0 V1 — k?sin? 6 0 (1+ nsin? AVI — k? sin? 6 


(62) 


II(n, k) can be expressed in terms of incomplete elliptic integrals,’ for which tables 
are available. Similarly, the substitution 
] 2 sin? 6 
y? = + aaa (63) 


transforms the integral in equation (59) into 
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L pre 10 9 
R.(7) = — , 7 ; \ vo tas + 

do Jo V1 — x? sin26@ | (8 sin? @ + 1) 
(3 + V3) (3 — V3) 
oe ae = : 
(8 sin?9-—5—3V3) (8sin?6 — 5 + 3V3)\ 

= — 7. [3K(w) — 9118, «) - (2V3 — 3)M[—(12V3 — 20), x] + 
Ww 


2V3 + 3)M[(12V3 + 20), x]}, (64) 


2 
2 = 65 
; 37? — | 7 


Substituting these into equations (57), we finally get 


g(r) = 0, t< : 
V3 
q(t) = V3/erZ {6K(k) — 1811(8k?, k) + (6 — 4V3)N[—(12V3 — 20)k2, k] + 
1677ur 
(6 + 4V3)N(12V3 + 20)k2, ky}, wees = 
g(t) = V8/oreZ {6K(x) — 1811(8, x) + (6 — 4V3)N[—(12V3 — 20), «] + 
1677ur 
(6 + 4V3)N[(12V3 + 20), «]}, peec y=, ¥8 + V3, 
g(r) = V8/srxZ {6K (x) — 1811(8, x) + (6 — 413) [— (1213 — 20), «] + 
167°ur 


Zr 


(6 +- 4V3)1[(2V3 + 20), «]} + : 
SaurV 7? — 7’ 


r>vy_ (66) 

5. Discussion of Results—The vertical displacement at the surface w(7r) due to 
the application of a point surface pressure pulse of the form H(t) was computed from 
equations (45)—(48), and is plotted in Figure 3. The applied pressure being down- 
ward, the steady-state vertical displacement is downward. The initial displace- 
ment, however, at the time of the arrival of the P-wave, is upward. We note 
that the shear wave is marked only by a discontinuity in slope of w(r). The dis- 
placement becomes infinite at the time of the arrival of the Rayleigh wave, and 
subsequently it reverts to the steady-state value. 

Similarly, the horizontal displacement g(r) was computed from equations (66), 
and is plotted in Figure 4. The ultimate displacement is inward, but the initial 
displacement at the P-epoch is outward. The arrival of the shear wave is marked 
only by a change of slope of g(7). which is even less marked than in w(r). Near 
the time of the arrival of the Rayleigh wave the infinity in ¢(7), as well as in w(r), 
varies with distance like 1/+/r, which is characteristic of surface waves. 
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Fic. 3.—Vertical displacement at the surface w(t) due to the ap- 
plication of a Heaviside unit pressure pulse H(t) at the surface. 
w(t) = -—(Z/arur)G(r), +r = (est/r). Z (negative) is the surface 
integral of the applied pressure. w(t) was computed from equations 
(45)-(48). P denotes arrival time of the compressional wave, S of the 
shear wave, and FR of the Rayleigh wave. 
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Fic. 4.—Horizontal displacement at the surface q(t) due to the 
application of a Heaviside unit pressure pulse H(t) at the surface. 
q(t) = —(Z/apr)E(r); + = (est/r) Z (negative) is the surface in- 
tegral of the applied pressure. g(t) was computed from equation 
(66). P, S, and R denote times of arrival of compressional, shear, 
and Rayleigh waves, respectively. 


6. Summary.—Exact and closed expressions are derived for both the hori- 
zontal displacement g(t) and the vertical displacement w(t) of the surface of a 
uniform elastic half-space due to the application at the surface of a point pressure 
pulse varying with time like the Heaviside unit function H(t). The applied pres- 
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sure pulse is specified in equation (2) and is such that its integral over the surface 
is finite. Both w(t) and q(¢) turn out to be, in this case of a surface pulse, of the 
form (1/r)f(r), where + = (ct/r), c denoting the shear velocity. The solution for 
w(t) is given in equations (45)—(48) and is plotted in Figure 3. The solution for 
g(t) is given in equation (66) and is plotted in Figure 4. Both w(é) and g(t) become 
infinite at the time of arrival of the Rayleigh wave, but the arriva] of the shear 
wave is marked only by a change in slope of the displacements. 


1H. Lamb, Phil. Trans. Roy. Soc., ser. A, 203, 1, 1904. 

2C. L. Perkeris, these PROCEEDINGS, 26, 433, 1940. 

3L. Cagniard, Réflexion et réfraction des ondes seismiques progressives (Paris: Gauthier-Villars, 
1939). 

4E. Pinney, Bull. Seis. Soc. Amer., 44, 571, 1954. 

5 C. H. Dix, Geophysics, 19, 722, 1954. 

6 H. Bateman and C. L. Perkeris, J. Opt. Soc. Amer., 35, 655, 1945. 

7P. F. Byrd and M. D. Friedman, Handbook of Elliptic Integrals (Berlin: J. Springer, 1954). 
p. 225. 


ON INVOLUTORIAL ALGEBRAS 
By A. A. ALBERT 
UNIVERSITY OF CHICAGO* 
Communicated May 18, 1955 


The main result of this note is the proof of a conjecture of A. Weil.! 
Let % be an n-dimensional associative algebra over a real field §. It will be con- 


venient to regard Y% as consisting of the set of all vectors x = (&,..., &) with 
co-ordinates £; in § and the product xy as being defined by 
zy = zk,, (1) 
where F, is an n-rowed square matrix whose elements are linear forms in the co- 
ordinates n; of y = (m,...,,). Since Y is associative, we have 
Res = R,R,. (2) 


A linear transformation x —~ 2x* of Yf over § is said to define an involution over § 
of 9 if 
(zy)* = y*z*, (a*)* =z, (3) 


for every x and y of Y. Assume that 2 > 2* is an involution of %& over § and that 
there exists a linear functional 6(z) on YI to § such that 


5(zx*) > 0 (4) 


for every nonzero x of A. Then we have the following result. 

THEOREM 1. The algebra Wis semisimple. 

For the radical 2 of % is the set of all properly nilpotent elements z of YW. If z 
is in N, we have (zx*)” = 0 for every x of A, [(zx*)"]* = (xz*)" = 0, and 2* is in 
N. Then y = z + 2* is in N for every z of N and y* = y. If y ¥ 0 for some 
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zin N, then y* = w ¥ O for some k, w = w*, w? = 0, ww* = 0, whereas 6(ww*) > 
0. Hence z* = —z for every z of N, w = zz* = —22 = w* isin N, w = —w* = 
w*, and w = 0, 6(w) > Ois impossible. Then z = 0, R = 0, Y% is semisimple. 

If a is any element of %, the algebra R.= § [a], of all polynomials in a with co- 
efficients in §, is isomorphic to Ry = F[R,]. Then the minimum function of R, is 
the minimum function 4(£) of a, and the characteristic roots of R, are the roots of 
¢(é) = 0. Evidently & is isomorphic to the difference algebra §[E] — (¢), where 
() is the principal ideal defined by ¢ = 4(€). 

THEOREM 2. Leta = a*. Then the complex characteristic roots of R, are all real. 

For a = a* implies that 8 = § [a] is an algebra such that 2 — 2* is the identity 
automorphism over § of B. Also, 6(x) is a linear functional on 8 to § and 


5(x?) > 0 (5) 
for every nonzero x of 8B. By the argument above, the characteristic roots of 
R, depend only on the minimum function of a, and this is the same when a is con- 
sidered as an element of 8 as when a is considered as an element of %. This re- 
duces our study to the case where A = Fla]. By Theorem 1, this algebra has a 
unity quantity e. 

Since 6(x) is linear, the function 6(xy) is a symmetric bilinear form on A = § [a]. 
Then 
d(zy) = xSy’, (6) 
where S is a symmetric matrix and y’ is the transpose of y = (m,...,%.). By rela- 
tion (5), the matrix S is positive definite. Then 
d(x) = d(xe) = b(er) = xSe’ = eSx’. (7) 
It follows that 
d(xy) = (ry)Se’ = rR,Se’ = xSy’, (8) 
and so Sy’ = R,Se’, 
ys = eSR;,. 
We now use the fact that % is commutative to write ySR, = eSR)R, = eS(R,R,)’ 
= eSR!, = eSR), = (yx)S = yR,S. This yields the fundamental relation 


RS = SR’ = (R,S)’. (10) 


Since S is positive definite, we may write S = PP’, where P is a real nonsingular 
matrix. Then R,PP’ = PP’R,, P-'R,P = P’'R{(P’)-' = (P-'R,P)'. We 
have proved that P~'R,P is symmetric. Its characteristic roots are all real, so 
that the characteristic roots of the similar matrix R, must all be real. 

We shall say that a is positive if 


a = a* £0, d(zax*) > 0, (11) 


for every x of AM. Ifa = \oB,a,-2,* ¥ 0 for B; > 0, then zax* = YBywra,*r* = 
SBwyiyi* and 6(xax*) = >8,8(yy.*) > 0 if at least one B; + 0. A. Weil con- 
jectured that if § = ® is the field of all rational numbers, then a is positive if and 
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only if a = >08;-x,;-2;*. We shall derive a somewhat stronger form of this result. 
We first prove 

THEOREM 3. Let a be positive. Then the nonzero characteristic roots of Rq are 
all positive. ‘ 

For we have already shown that we may restrict our attention to the case where 
% = §la]. By hypothesis, 6(xax) > 0 for every x of YA, and so we see that 5(zar) = 
zaSx’ = xrR,Szxr' > 0 and the symmetric matrix R,S of equation (10) is positive 
semidefinite. If x is a characteristic vector of R, and a@ is a characteristic root 
such that xR, = az, we have 6(xar) = arSzx’. But 2Sz’ > 0, zR,Sz’ > 0 for all 
real nonzero vectors 2;, and so azSx’ > 0, a > 0 as desired. 

We now derive the final result. 

THroremM 4. Let % be an involutorial algebra over the field N of all rational num- 
bers, and let 6(x) be a functional satisfying relation (4). Then an element a = a* of 
YW is positive if and only if a = ay? + xe? + x3? + ay? for x, = x,* in Ra]. 

For ® [a] is semisimple and so may be written as the direct sum R[a] = Rai] @ 
... ®Ria,] for algebraic fields R[a;] over R such thata = a, +...+a. If 

4 4 t 
a, = ), 2,?thena = >> 2,2, wherez; = >> 2. Also, (z,az,) = 6(z,a,2,) > 0 
j=! j=l i=1 
for every x; of R{a;], 6(27,°) > 0 for every x; ¥ 0 of Rla;]. Thus we have reduced 
the proof of our theorem to the case where A = R{a] is a totally real algebraic num- 
ber field. Since a is totally positive, the quadratic form f(y) = yi? + ye? + ys? + 
ys? — ays’ has the property that all algebraic number conjugates are indefinite. 
By a theorem of H. Hasse,’ the form f(y) is a null form, that is, y:? + yo? + ys? + 
ya” = ays? for y; not all zero and in N{a]. But ys ¥ 0, since otherwise y:? +>. . . 
+ ys = 0 for real y; not all zero. Hence a = 2,2 + ao? + 23? + 24°, with x; = 
yilys) —?. 

* This paper was sponsored in part by the Office of Ordnance Research, United States Army, 
under Contract No. DA-11-022-ORD-1571. 

1 The conjecture is connected with Weil’s attempt to generalize the results of H. Morikawa, 
“On Abelian Varieties,’’ Nagoya Math. J., 6, 151-170, 1953. A special case of the theorem of this 
note is proved there, but the present proof is more general and of a more algebraic nature. 

2 J. f. Math., 152, 113-130, 1923. 


FUNCTIONAL EQUATIONS IN THE THEORY OF DYNAMIC 
PROGRAMMING. IT. NONLINEAR DIFFERENTIAL EQUATIONS 
By RicHarp BELLMAN 
RAND CORPORATION, SANTA MONICA, CALIFORNIA 


Communicated by Einar Hille, May 19, 1955 


1. Jntroduction.—In the first paper of this series! (ef. also other papers*) we 
discussed the existence and uniqueness of solutions of the functional equation 


f(P) = Max f(T(P, Q)), (1.1) 
Q 
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the basic functional equation of the theory of dynamic programming. In this paper 
we shall present the first of a series of investigations of various continuous analogues 
of the foregoing equation. Our results here concern the nonlinear system 


dx; 
= Max f(x, t; q), 2(0) = c, i(-E 24m (1.2) 
dt q 


which we shall write in the vector form 


= = Max f(a, t; q), x(0) =e. (1.3) 
dt q 

The quantities x, f, and c are n-dimensional vectors, while g is an m-dimensional 

vector whose components may range over a finite or denumerable set, or continuum 

of values. 

In subsequent papers we shall present analogous results for partial differential 
equations of parabolic type, for the first-order partial differential equations ob- 
tained from the calculus of variations,* * and for other classes of functional equa- 
tions as well. A detailed exposition of the results of this note and a discussion of 
their relation to an interesting class of problems in the calculus of variations will 
appear elsewhere. 

Equations of the above type appear in the study of continuous decision proc- 
esses; some particular examples have been treated in unpublished notes by Cope- 
land and Darling. In many ways these equations constitute natural generaliza- 
tions of linear systems and thus play an important role in analysis, apart from their 
probabilistic applications. | Furthermore, as we shall see below, some important 
types of equations, such as the Riccati equation, can be written in this form. 

2. A General Existence and Uniqueness Theorem.—Let us first state a general 
result. 

THEOREM |. Consider the system in equations (1.2). Let us assume that f(x, q, t) 
satisfies the following conditions: 

a) ||f(a,t;q) — f(y, @\| < k(q, d/|2 — yl|, for x, y, t in the region R defined by the 
inequalities lle — ell <a, ly —e) <a,0 <t < ty, for some nonzero constants a and 
lo. 

b) The maximum value of f(x, t; q) is assumed for a function q; = qi(x, t) which is 
uniformly bounded in R. 

c) The function k(q, 1) is uniformly bounded for uniformly bounded q and for tin R. 


Then there is a unique solution of equation (1.3) for0 < t < t, where t, is some non- 
zero constant less than to. This solution may be obtained as the limii of the sequence 
defined by 


% = C, 
Lass = C + fi Maz f(xn, 8; q) ds. (2.2) 
q 


There are many interesting questions concerning the structure of the solution x 
and the policy function g which we do not have space to discuss here. 

3. Quasi-linear Systems.—As usual, we can obtain'a much more complete result 
for systems of linear type. 
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THEOREM 2. Consider the system 


la, N 
ages Max | ot ag t+ >¥ a(t; axs| 20) = ¢,, o=i-2...,%, (ot) 
=1 


d qd 
where we assume that 


Maz |ai(t;q|, Max |b(t,q| < fd, (3.2) 
¢ q 


2 


with f(t) integrable over any finite t-interval. 

Then there is a unique solution for t > 0 which may be found by the same method 
of successive approximations described above. 

In many cases we can determine the asymptotic behavior of the solution (see a 
forthcoming paper® for a discussion of this problem for difference equations). Re- 
sults of similar type may be established for differential-difference equations and 
more general types of integro-differentia] equations. 

4. Approximation in Policy Space-—The approximation technique used above 
was the classical one. Let us now describe another type. Given q(¢), for each 
index 7, we may compute an initial approximation, 2°, using the equation 


0 


d. . 
s = A(t; qo)x® + b(qo)z, 490) =. (4.1) 


To obtain a better approximation, we determine q,(¢), again for each 7, by the 
N 


condition that q,(t) maximize the expression b a(t; dae + bi(t, o| Having 
j=l 
determined q,(¢) in this way, we compute a second approximation, x!, using equation 
(4.1), with q replaced by q. We now continue in this fashion, computing a se- 
quence of policy functions | q,} and a sequence of approximations {x"}. This type 
of approximation we call ‘approximation in policy space’. Concerning this 
process, we can establish the following result: 
THEOREM 3. /f the linear system 


d. 
= =A(t,ga+f(t),  2x(0) =0 (4.2) 


has, for any fixed q, the property that f > 0 for 1 > 0 implies that x > 0 for i > 0, then 
approximation in policy space yields monotone convergence to the solution. 

In particular, there is always convergence when x ts a scalar variable. 

If A(q) is a constant matrix for each allowable q, then a necessary and sufficient 
condition for the above condition to hold is a,; > 0, for 7 ¥ j; if A(t, q) is variable, 
then a sufficient condition is a,;(t) > 0, for7 # 7.° 

5. The Riccati Equation.—An important tool in the qualitative and quantitative 
theory of the second-order linear differential equation vw” + a(t)u’ + db(t)u = 0 
is the transformation of the equation into the form 


vp’ = —v? — a(t)v — O6(b), (5.1) 
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the Riccati equation, by means of the substitution vu = exp (fv dt). We now ob- 
serve that the Riccati equation may be written in the form 


v’ = Min [q? + 2quv(t) — a(t)v — b(t], (5.2) 
q 


proach to the study of the solutions of equation (5.1). There are several other 
classes of equations of both analytic and physical interest to which this quasi- 
| linearization may be applied. This will be discussed elsewhere. 

i 

} 





} 
an equation of the above general type. This representation furnishes a new ap- 
} 
U 


i 'R. Bellman, “Some Functional Equations in the Theory of Dynamic Programming,” these 
i PROCEEDINGS, 39, 1077-1082, 1953. 
| 2 R. Bellman, An Introduction to the Theory of Dynamic Programming (RAND Monograph 
R-245, 1953); “Some Functional Equations in the Theory of Dynamic Programming. I. Func- 
tions of Points and Point Transformations,” T'’rans. Am. Math. Soc. (to appear). 
3 R. Bellman, “Dynamic Programming and a New Formalism in the Calculus of Variations,’ 
| these PROCEEDINGS, 40, 231-235, 1954. 
4R. Bellman, ‘“Monotone Convergence in Dynamic Programming and the Calculus of Vari- 
ations,’’ these PROCEEDINGS, 40, 1073-1075, 1954. 

5 R. Bellman, ‘‘Quasi-linear Equations”’ (to appear). 

6 R. Bellman, I. Glicksberg and O. Gross, ‘‘SSome Variational Problems Occurring in the Theory 
of Dynamic Programming,”’ Rend. Palermo (to appear ). 





STATIONARY SPACE-TIME IN GENERAL RELATIVITY 


By 8S. BocHNER 
PRINCETON UNIVERSITY 
Communicated May 20, 1955 


In the notation of the recent book of A. Lichnerowicz! we take a V4 with an 
everywhere regular hyperbolic local line element which, somewhat differently from 
the book, we denote by 


3 
>, bag Gt* dx’. (1) 
a, 8=0 
} We assume that V, is stationary, as there defined, and this implies that V4 is the 


topological product of Vs; with R: (— © <t< @) and that on V; the line element 
(1) induces certain data which on V; are linked by the relation 


] 
At = —tR + 4H, (2) 


with which we will be concerned exclusively. Underlying it is a positive-definite 
Riemannian metric, 





3 
Lo gy dx' da’, (3) 


\ 
1} i,j=1 








which is induced by expression (1). The quantity ¢ is a scalar function which is 
positive throughout, 
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é> 0, (4) 
and Aé is the Laplacean g”é,, ;; R is another scalar, which is denoted in Lichnero- 
wicez with indices because it is induced by the Ricci tensor Rg pertaining to expres- 
sion (1) and which is such that, if 


Rog = 0 
throughout, then also 
R=0 (5) 
thoughout; and, finally, 
H? = HyH" > 0, (6) 


where H ;; is a certain skew tensor. 
Definition 1: We will say that V; is empty if on V3 we have 


H = 0, £ = const., (7) 
and we will say that it is semi-empty if only 
H =.0. (8) 


Remark 1: The space is empty if and only if expression (1) is locally Euclidean, 
but this will not concern us. Also, relation (2), always on a V3, also occurs in some 
field theories with suitable interpretation of the symbols, whether the theory be 
based on a V, or on V;, k > 5, and relation (7) or (8) always implies a certain man- 
ner of emptiness. 

In what follows, we are assuming that V3 is C?, g;; is C!, is C?, Ris C°, and H,,; is 
c*. 

Important criteria for the occurrence of relation (7), or perhaps only (8), are 
due to, among others,' Levi-Civita, R. Serini, Ch. Racine, Lichnerowicz, Einstein 
and Pauli, and Yvonne Fourés and Lichnerowicz (for a recent result in a related 
but different direction see a paper by A. H. Taub’), and it is the purpose of the pres- 
ent paper to comment on one of these criteria by way of a generalization. The 
generalization as such has been given before* in our work on “Curvature and Betti 
numbers” but will be restated for the present context. 

THEOREM 1. Let V3becompact. (2) If we assume 


—Sv, ER dv > 0, (9) 
then V4 is semi-empty, and equality holds: 
—Jfy, ER dv = 0. (10) 
(a2) If, more precisely, 
—Rk > 0, (11) 


then V4 is empty, and, in fact, equation (5) holds. 
Proof: If we integrate equation (2), we obtain 


l 
i] Atdv = — i) ER dv + ri HH? dv. (12) 
Vs Vs Vs 
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; On the other hand, since V3 is compact, we have* 
Sv, AE dv = 0, (13) 


and relations (6), (9), (12), and (13) imply (8) and (10) as claimed. If, more pre- 
cisely, relation (11) holds, then (2) implies AE > 0, and the desired conclusion — = 





} 

| 

{ const. follows now either from the formula‘ 

j Sin & + Cdk dv = —So.g" Baes dv (14) 

| or from the following local theorem due to E. Hopf. 

LemMMA. In a co-ordinate neighborhood V, if a function &(P) of class C? satisfies 

the inequality AE > 0 and if there is a fixed point P® in V such that &P) < &P®) 

| everywhere in V, then we must have &(P) = &(P°) everywhere in V. 

We are now turning to comment on Theorem 1, in case V; is not (necessarily) 
compact. 


Assumption 1: There is given on V3; a group [':(y) of one-to-one continuous 
transformations (with fixed points, perhaps) having the following property. There 
{ exists a compact set S° in V3 such that, corresponding to any point P® in V3, there 
is an element y° in I such that the transformed point y°P® lies in S°, 

Definition 2: Let f(P) be a continuous function on V3, either a scalar or, more 
generally, any tensor density. We call it ‘almost periodic” (or almost automor- 
phic), relative to the given group I, if every infinite sequence of elements {| y,} in 
I contains an infinite subsequence } y,,{, Pp = 1, 2,..., having the following prop- 
erty. Corresponding to any point, there is a co-ordinate neighborhood in which 
the sequence of functions f(y,,/?) converges uniformly in each component, and, if 
we denote the limit function by /*(P), then, reversely, the sequence of functions 
f*(¥np~'P) converges, in the same fashion, toward the function f(P) from which we 
have started. 

Remark 2: The italicized part of definition 2 makes it unnecessary to require 
that (the space V; shall have uniform structure and that) the sequences shall con- 
verge uniformly in all of V3. 

It is easily seen that sums, products, and contractions are again almost periodic 
and that for an almost periodic scalar ¢(P?) the number 





| 
| 
| 


m = supp ¢(P) (15) 


is finite (and also the inf.). 
THEOREM 2. Under assumption 1, if the functions 


Ji; g", g, Ei, Ei. 55 Rk, H;;, (16) 


are all almost periodic with respect to T, and if we have R < 0, then V4 is empty. 

Proof: Relation (11) implies Ag > 0. Now, if the supremum (15) is attained, 
that is, if there is a point P® for which &(P°) = m, then the lemma implies ¢ = const. 
everywhere. But relations (2), (6), and (11) then also imply H = 0 as claimed. 

If, however, it is not attained, then there exists a sequence of points {P,} for 
which &(P,,) > m. By requirement | there exist elements y, ¢ [ such that y,P, > 
S°, But S°is compact, and if we apply definition 2 judiciously, we may assume that 
the sequence | y,} has the following properties simultaneously. We have y,?,, > 


P°, where P® is a point of S°, and the 7 sequences (see expressions [16]), 
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Jis(¥nP rn), g"(Y¥nP nr); E(val’s); ai? 2% 


all converge uniformly in the neighborhood of every point of V3. If we denote the 
limit functions by asterisks, g;;*, g*"’, &*, R*, etc., then we again have 


l 
A*e* = —£*R* 4 ; Fad | ha 


and the function &* does attain its supremum. Therefore, ¢* is constant; but, by 
the reverse part of definition 1, the function € is likewise a constant, and this com- 
pletes the proof of the theorem. 

If V; is the Euclidean £; and I the group of all translations, then assumption | 
is fulfilled, with S being any one single point of 23. Also, in this case a function ¢(2), 
x = (z', x*, x*), is almost periodic in the sense of definition 2 if it is almost periodic 
in the sense of Bohr,' that is, if and only if it is a limit, uniformly in £3, of finite 
sums 


t(Arz! + Aer? + Asx) 
rd Ay rodrs © 
r 


for arbitrary real numbers Aj, Ae, As. ‘Theorem 2 now implies the following: 
THEOREM 3. Assume that V; ts the Euclidean FE; and that we are given one global 
line element (3) on it. If the individual functions 


r Ogi; 0 Gij fe) 0? 
qg@ (x), Gis(X), J * J * gE é <a R, Hi, 
are all almost periodic in the sense of Bohr, and if R < 0, then V4 is empty space. 


An almost periodic function f(x) on a group has a mean M{f(x)} which on £3 can 
be evaluated as 


1 T T T 
Mi f(x)} = tim (27): Jf se x?, x*) dx! dx? dx’, 


say. Now, for a domain D in £; with sufficiently smooth boundary B, Stokes’s 


formula gives 
of 
/ . (= — di 2 ( 7 
[ AE-dt ® a Oe (17) 


but if Dis the cube —T < x' < 7,7 = 1, 2,3, then, for T—> o~, 


Sz |dw| = 0(T?), (18) 
and 
dé 
dn o(1), (19) 
and therefore 
M{ de} = 0. (20) 


However, relation (2) implies 








SS 


Se es ee 
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1 
M{ de} = —M{ER} +4 MEH, (21) 
and 
f2>0 and Mf =0 imply f =0. (22) 


Hence the following conclusion. 
THEOREM 4. Jf in Theorem 3 we replace relation (11) by the weaker assumption 


—M{éER} > 0, (23) 


then V3 ts semi-empty. 

Remark 3: If in Theorems 3 or 4 the line element (3) is not global but only local, 
then we must assume that functions (16) are almost periodic in the sense of Bohr, 
and the conclusions remain in force. For Theorem 3 this is a direct consequence 
of Theorem 2, and for Theorem 4 this follows from the fact that, if we put 


Aé-dv = Atv/g dx dx? dz’, 


then the integrand is almost periodic in the sense of Bohr, and similarly for the other 
integrands in equation (21). 

A meaningful extension of Theorem 4 to a V3; with a general group I is not easy 
to obtain; but we will make a statement, a somewhat technical one, for the case in 
which the topological quotient V3;/T is itself a (compact) manifold, which we will 
denote by U3. In this case the group T is countable, and we enumerate its elements 
as a sequence vy’, y', y’,..., in which 7’ is the identity, and we decompose V3; into 
a sequence of disjoint “fundamental domains” U°, U!, U?,..., which are such that 
y" transforms U° into U”. They are not really “domains,” that is, open sets, but 
their closures are (compact) closures of domains. From relation (2) we obtain 


l 
f Aé dv = -f ER dv + f & H? dv, (24) 
wy uy@™ 4 uy” 
and we note that for any function f(P) on V3 we have 
Sum {(P) de(P) = Su f(y"P) dv(y"P). (25) 


Now, the integrand on the right side is the product of dx! dx*® dx* with f(y"P) 
(g(y"P))'”, and, if the last function is almost periodic in our sense, then it is almost 
periodic on IT in the sense of von Neumann, uniformly for P in U®. Such a fune- 
tion has, again, an almost periodic mean M} } again with property (22), and 
equation (24) implies 


1 
M, ' f At avt 2 1 f ER ao +=M, { f eH? aw. (26) 
Ue Ue 4 ue 


Hence the following conclusion. 
TuHeoreM 5. [f it is known that, in analogy to equation (13), we have 


M,\ Suv At do} = 0, (27) 
then the assumption 
MS ER dv} < 0 


implies that V4 is semi-empty. 








490 MATHEMATICS: BOTT AND SAMELSON Proc. N. ALS. 


This theorem is useless unless equation (27) can be established, and, if one wishes 
to prove this by imitating the procedure used for Theorem 4, then the following 
features have to be verified, subject to assumptions. It is possible to arrange the 
ordering | y”} in such a manner that for a sequence of integers 7, > © we have 

T. 
M\f(P)} = lim , > f(y"P) (28) 
n—>o nm=t1 
uniformly for Pin U®. Also, the union of sets 
Tn 
> U 
m=0 
constitutes a piece of space whose actual boundary consists of several sufficiently 
smooth surfaces whose two-dimensional! induced surface area is 0(7',) asn > ©. 

Remark 4: Finally, we note that, with the aid of the mean (27), we may even 
obtain the statement on emptiness itself without requiring almost-periodicity at all. 
Assume that for a certain mean (28) we may go over from equation (25) to equation 
(26); that relation (22) holds; and that we have not only equation (27) but also 


My Soo & AE dv} = —M, Sing! Eé;-dd} 
in analogy to equation (14); then —R > 0 implies that V4 is empty. 

1 Theories relativistes de la gravitation et de l’électromagnétisme (Paris: Masson & Cie, 1955). 

2“Empty Space-Times Admitting a Three Parameter Group of Motions,’”? Ann. Math., 53, 
472-490, 1951. 

3 “Vector Fields and Ricci Curvature,’’ Bull. Am. Math. Soc., 52, 776-797, 1946, esp. pp. 782- 
83; K. Yano and 8. Bochner, Curvature and Betti Numbers (Princeton, N.J.: Princeton Uni- 
versity Press, 1953), pp. 182-185. 

4 “Remark on the Theorem of Green,”’ Duke J., 3, 334-338, 1937. 

5 “Beitrige zur Theorie der fastperiodischen Funktionen. IJ. Funktionen mehnerer Vari- 
ablen,”’ Math. Ann., 96, 383-409, 1926. 


THE COHOMOLOGY RING OF G/T 
By R. Borr anp H. SAMELSON 
UNIVERSITY OF MICHIGAN 
Communicated by C. Chevalley, May 12, 1955 


1. In this note we describe a method for constructing the cohomology ring 
H*(G/T) with integral coefficients, for an arbitrary compact Lie group G, with 
maximal torus 7’. This refines results due to A. Borel,' who found the cohomology 
ring of G/T mod p under the assumption that G has no p-torsion. We make strong 
use of the fact that the cohomology ring of G/T with rational coefficients is gener- 
ated by H*(G/T). On the other hand, our procedure demonstrates anew that 
G/T (as well as G/C, where C is the centralizer of any torus of G@) is torsion-free.’ 

Our procedure is as follows: We construct a manifold B, whose cohomology ring 
we can determine in terms of the diagram of G; in particular, there is no torsion; 
dim B = dim G/T; and we construct a map f of B into G/T. 
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Proposition |. The map f of B into G/T has degree one. 

Coro.Luary. G/T has no torsion (since B has no torsion). 

THeorEM I. H*(G/T) is isomorphic to the smallest subring of H*(B) which is 
additively a direct summand and which contains f*H?(G/T). 

This is a consequence of Proposition I and the fact that H*(G/T) is rationally 
generated by H?(G/T). 

To determine H*(G/T), we have then to describe H*(B) and to identify 
f*(H?(G/T)); this will be done in Proposition IT. 

2. We recall briefly the concepts connected with the diagram of G; we may and 
shall assume G to be semisimple. JL is the Lie algebra of G, 7’ a maximal torus, H 
the Lie algebra of 7; dim 7 = /, dim G = n = 1+ 2m. The diagram > consists 
of the elements of H whose images in G under the exponential map are singular, 
i.e., have centralizers of dim > /; it consists of m families of parallel equidistant 
planes, which partition H into compact convex cells. The dimension of the cen- 
tralizer corresponding to x « H is/ + twice the number of planes of = going through 
x. H carries a natural Euclidean metric, given by the Killing form; we denote 
the inner product by a dot. The lattice A is obtained by reflecting the origin o in 
all planes of =. Let =’ be the subset of © through 0. We denote the planes in 
>’ by Pi, ..., Pm in the order in which they meet some directed straight-line seg- 
ment \ in H, which meets them all and does not go through 0. Each P; determines 
an element 6; « A, obtained by reflecting o in the next plane parallel to P; and such 
that 6; > 0. Under the exponential map H — 7, the segment from o to 6; maps 
into a simple closed curve (closed 1-parameter group and geodesic) in 7, whose 
homology class in H,(7') we denote by y;. Let }t,} be a subset of {yi}, which is a 
base for H,(7’) (corresponding to a fundamental system of roots). This defines 
integers a;,,2 = 1,...,m,a=1,..., l, by 


. 2Bigta 


Finally, we recall the fact that the numbers a;; = (20;-6;)/(0;-0;) are integers, the 
Cartan integers of G. 

3. We come now to the construction of B and f. Let C;,7 = 1,...,m, be the 
component of the identity of the centralizer in G of the image, under the exponential 
map, of P;. We form the products W* = (C; X... XK Cy, k = 1,...,m. Since 
T ¢ C; for all 7, we can let 7’ = TX... X T (k factors) operate on W* on the 
right by the law 


(C1, mar rae Cr) (4 i a ee l,) = (Cyd), t leolo, lo leats oe. wae d.—1—'e,4;). 


(This skew operation was suggested to us from Morse theory considerations; it rep- 
resents roughly variations of segments between successive singular points on a geo- 
desic in G.) This defines W* as a principal 7*-bundle; let B* be the corresponding 
base space. One verifies easily that the projection W*+! —~ W* (by suppressing the 
last factor) induces a map of B*+! — B*, which is a fibermap with a 2-sphere as 
fiber and which has a cross-section. This makes it possible, following Gysin,* to 
compute *(B**!) in terms of 1*(B*) and the characteristic class of an associated 
l-sphere bundle; it is this class which brings the homology classes y; and the in- 
tegers a, into our argument. The space B promised above is B”. Let f: W" > 
Gand f’: 1" — T be defined by f(c, . - . , Gm) = CrC2. . - Cm S(t, . . ~ 5 tm) = tm 
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If we think of G as a principal T-bundle, with base G/T, then f defines a bundle 
map of W” into G over f’ and induces a map f of B” = B into G/T; this is the map 
f of Proposition I. 

Proposition II. With a suitable base }x1,..., Xm} of H(B), the algebra H*(B) 


becomes Z (x1, . . . , Xm|/I, where the numerator is the polynomial ring in the variables 
Li, ..., Um over the integers Z and where I is the ideal generated by the elements 
b~1 
Pr = t+ >> antes, eo haa m. 
I 


Further, f*(H?(G@/T)) ts generated by the elements 
m 


Ze = ) Bets, pets. iet):, 1. 
1 


4. The main step in the proof of Proposition I is the observation that the point 
T of G/T has exactly one inverse image point under f or that the inverse image of 7’ 
under f is precisely the set 7” ¢ W”. To prove this, we argue as follows: Let p; 
be the point of intersection of \ with P;; let po, resp. p*, be the initial, resp. terminal, 
point of \. We map W” into the space [ of polygons in L, starting at po, by the 
following law: (¢, . . . , ¢m) goes into the polygon { pop:, Ad e:(pyp2), Ad cyco(pops), 

., Ad eee... Cm(Dmp*)}. One verifies that this ordered collection of line segments 
in L is a polygon from pp to Ad cco. . . Cm(p*). Now the orbit of p* under Ad G can 
be identified with G/T, and with this identification the map (c, ... , Cm) ~ Ad 
(C2 . . . Cm(p*) is the composition of f with the projection G—~> G/T. Since Ad G 
leaves the Killing form invariant, all these polygons have the same length as X. 
But the polygon corresponding to (¢), . . . , Cm) has a corner at p; unless c; € 7’; there- 
fore, the inverse image of p* is 7” c¢ W™”, as was to be shown. 

5. As an example of our method, consider the case of G:. Here, for a suitable 
\, H*(B°) is described by the relations 


x," = 0, 
Le? = DT 
x3" = (a1 + 322)23, 


R 
Cd 
I 


— (Xe + X3)X4, 
Xs” = (2 ace; ame 324)Xs, 
Xe? = (X%1 + Lo — Yq — Is)Ae, 


while f*H*(G/T) is generated by | a, 6}, with 
a= X% + 3x. + 223 + 3x4 + 4, 
B = x — 2%. — 3 — Xs + NX. 
An admittedly tedious computation yields first of all that 
a’? + 36? + 306 = 0, 


so that a? is divisible by 3 and H*(G./T) is rationally spanned by products of the 
form af* and g*. Further relations are a® = 6° = 0; moreover, a6 = 2x22... r% 
(= twice the fundamental cocycle of B), and 8° is the first element divisible by 2. 
A basis for H*(G./T) is then given by the elements 


{1, a, B, aB, B?, a6, B°/2, a6*/2, 64/2, aB'/2, B°/2, a*/2}. 
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A routine computation with the spectral sequence of the fibering (Gz, 7, G2/T) 
confirms the result of Borel‘ that H*(G.) has a cyclic group of order two as torsion 
group in dimensions 6 and 9. 

6. Let K be a connected subgroup of maximal rank in G. The inclusions 
T ¢ K c G induce a fibering f,; G/T ~ G/K of G/T over G/K with fiber K/T. 
If P,,,..., P., are the planes through o in Y’ which belong to K, set X; = %,%q,-- - 
a, in H*(B). Our imbedding f* of H*(G/T) — H*(B) now defines a homomor- 
phism x: H*(G/T) > Z by x (u) = coefficient of X, in the expansion of f*(u) in 
terms of the multilinear base in the z;’s. 

By Poincaré duality, x, can be realized by a definite element k « H*(G/T) once 
G/T is oriented and H?"(G/T) thus identified with Z. 

Proposition III. f,* maps a fundamental cocycle of G/K onto k. 

Corotiary. If k ¢ H*(B) is not divisible by a prime p, then G/K has no p-torsion. 

This proposition follows from the fact that the image under f,* of the funda- 
mental class of G/K is the Poincaré dual (in G/T) of the fiber K/T. 

1 A. Borel, “Sur la cohomologie des espaces fibrés principaux et des espaces homogénes de 
groupes de Lie compacts,’”’ Ann. Math., 57, 115-207, 1953. 

2 R. Bott, “On Torsion in Lie Groups,” these ProcEEpINGs, 40, 586-588, 1954; A. Borel, 
“Kahlerian Coset Spaces of Semisimple Lie Groups,” these ProcegpinGs, 40, 1147-1151, 1954. 

3. W. Gysin, “Zur Homologie-Theorie der Faserungen und Abbildungen von Mannigfaltigkei- 
ten,”’ Comm. Math. Helv., 14, 61-122, 1941. 


4A. Borel, “Sur l’homologie et la cohomologie des groupes de Lie compacts connexes,’’ Am. 
J. Math., 76, 273-342, 1954. 


JEOMETRIC ASPECTS OF CURRENTS AND DISTRIBUTIONS 
By James EELLs, Jr. 
INSTITUTE FOR ADVANCED STUDY 
Communicated by H. Whitney, April 28, 1955 


The theory of currents! on a manifold M is a calculus of abstract linear func- 
tionals, based on a generalized Stokes’s theorem. These functionals are defined 
on a certain topological linear space of exterior differential forms on M; analytical 
operations on the currents are defined by duality (i.e., by their effect on the differ- 
ential forms). 

In the present note we outline a procedure which permits a direct definition and 
geometric interpretation of the currents and of operations on currents. In the 
case of the currents of degree zero (i.e., the distributions) the definitions can be 
modified to give, for example, a theory of distributions on locally compact groups. 

The essential parts of the material summarized here are contained in the author’s 
doctoral dissertation, written under the direction of Professor Hassler Whitney. 
The present work rests heavily on the ideas developed by Whitney in his forth- 
coming book.? 

1. Let M be an ~-smooth (i.e., class C” ) n-dimensional Riemannian manifold; 
actually, our results are independent of the Riemannian structure. Let C,(M) de- 
note the linear space (over the real number field) of all ©-smooth singular r-chains 
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of M. Given a tangent vector v at p « M, we let |v| denote its length; if v(p) is 
an ©-smooth vector field with compact support on M, we let |v| = sup {|o(p)|: 
peM}. If ois a singular r-simplex and » is a vector field (both ©-smooth) on M, 
we define 7,0 as the translate of o by v (each point p of o is moved the distance 
|v(p)| along the trajectory determined by v(p)). 7,0 is always defined if | | is small 
enough; we can extend the definition of 7’, to an endomorphism of C,(/). 

For each m = 0, 1, 2, . . . we define the seminorm ),, on C,(M) as the largest 


seminorm A satisfying a, b, co, ... , Gm: 


a) \(c) < Jol, for all r-simplexes ¢ on M, where lol, denotes the Lebesgue 
r-dimensional measure of ¢ on M; 

b) (Or) < | r|,41for any (r + 1)-simplex 7 on VV, where 07 denotes the bound- 
ary r-chain of 7; 

badass shan 


rF+1)...@7@+k+ 1) 


Pall 
Cx) Mil, =. ¥ sd T)c) < | F\r 
for all r-simplexes o and all vector fields v;, where (7, — To) (7,, — Ted) = Tro 


— T,,0 — T,,0 + a, ete. 
Using the linear functionals on C,(M) which are bounded in X,,, we can establish 


the following formula: 
Given A ¢ C.(M), we consider all expressions of A of the form 


A= Yao; + >) Or; a par i: Pe ‘ioe T 0) 00; +- p thas + > Ris I, "9 (Tras — To) Omi: 


Then, minimizing over all such expressions for A, we have 


Am(A) o inf {5 a, lol, + >|: py ee + Ray 





| | Ue | 
Ao; | | vo:| | oo: | Omi 


lr 4 re 7 | v0; | ae | mil | omi| 7) 


r+] ied). + + 1) pan 


+> 


The case m = 0 was discovered by Whitney.* 

We can use the theorem of Section 2 below to show that A,, is a norm on C,(.V/); 
let €,"(M) denote the completion of C,(V) inX,,. From the definition of \,, we ob- 
serve that \o(A) 2 Ai(A) 2 ... forall A €C,(); we can identify €,"(M) as a linear 
subspace of ©,""*'(.M) for m = 0,1, 2,.... Let €,(M) = VU moo6,"(M). 

Proposition. For any A ¢ C,(M) we have Xm41/0A) < j,,(A). Hence we can 
extend 0 to be a bounded mapping of ©,"(.M) into cmt'(M). 

This can be proved with the help of formula (3) below. 

2. Let ©" "(M) denote the conjugate space of ©,"(M); let (CW) = Aa pio 
¢&” "(M). In the following theorem the term “derivative” refers to the Lie deriv- 
ative V,é of an ©-smooth r-form ~ with respect to the vector field v. 

THEOREM. Given X ¢ G’(M), there is a unique ©-smooth r-form — on M such that 
all derivatives of all orders <m are uniformly bounded (by a constant depending on m), 


and for any r-simplex o on M we have 
X-o =f, é. (2) 


Conversely, every such r-form on M defines, by formula (2), an element of ©'(M). 
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The construction of from a given X is fundamental in Whitney’s theory see 
Whitney,” * where ~ is shown to exist under very weak hypotheses). In order to 
establish its differentiability and boundedness properties, we first show that if X «¢ 
¢” "(M), then the “mth Lipschitz constant” 

| => 
) |X°1,9(7., — To)o| 
¥n( X) me sup ; | | a } | > 


(all o and v,; ¥ Of 


| | | | 
0} - ++ | Um Clr 


is finite. It follows that & has bounded Lipschitz constants of orders < m, which fact 
in turn implies that £ has derivatives of orders <m, each satisfying a Lipschitz 
condition of order 0 on VM. 

From this theorem it follows that if X « €’(M), then its co-boundary 6X (de- 
fined by 6X-A = X-OA for all A € C,(.M)) is in @’+!(M) and corresponds to the ex- 
terior differential dé of — Also, if X « ©” '(M), then 


Y m—1(6X ) < (r + 1)%,(X). (3) 


3. Incase r = 0 and J is an open subset of Euclidean n-space, we can identify 
S (VM) with the space of summable distributions’ of L. Schwartz. This can be seen 
from the following representation theorem or by studying the duality of the corre- 
sponding topological linear spaces. In general, €,(./) coincides with the space 
of summable r-currents. (Actually, de Rham treats the r-currents on M as (n — r)- 
forms with distributions as coefficients; we prefer the contravariant formulation. ) 
The exterior differential of a current corresponds (except for sign) to the boundary 
of the associated chain. 

THEOREM. Given A e¢ &,(.\), there exists a finite set of summable r-vector fields 
ae a, and vector fields v;;0n M such that 


X+A =Sy Li Vos, . Vous &(p)-an(p) dp (4) 


for all X ¢« &'(M), using the correspondence (2); the representation (4) is not unique. 
Conversely, every such representation defines an element of ©,(M). 

4. The cup and cap products of algebraic topology can be defined in the spaces 
("(M), ©.(M) and correspond through formulas (2) and (4) to the exterior and in- 
terior products of the r-forms and s-vector fields on M. In particular, if r = s = 0, 
then the cap product corresponds to the “‘ multiplication of distributions.’’4 
5. If M X N denotes the Cartesian product of two manifolds M and N, then 
LY p+q=-rSp(M) ® ©,(N) is contained in ©,(M X N); in fact, the tensor product @ 
is a bounded operation in each \ ,,. 

In particular, if M is a Lie group, we can use the tensor product to define the 
convolution product A * B ¢ G,,,(.M) of A € ©&,(M) and Be &(M). If A and B are 
singular chains on M, then 


X-(A *B) =f, fe t(zy) for all X ¢ G'**(M). (5) 


The operation * is bounded in X,, and is associative. By taking for B an ©-smooth 
function with compact support, we find that A * B is an ©-smooth summable r- 
vector field on M, for any A € €,(M). 

If A = aand B = B are summable r- and s-vector fields on M, then @ * 8 corre- 
sponds to the exterior product a V 6 as follows: 
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X+(a* B) = firxmi(zy)+[a(x) V B(y)] dx dy forall Xe @'t*(M). (6) 


Formula (6) and Section 4 relate the convolution product (which coincides with 
the Pontrjagin product® on the homology classes) in the homology ring of a Lie 
group to the cup product in the cohomology ring. 

6. With suitable modifications of the definition of \,,, we can define families of 
seminorms on (,(/) which will give other spaces of currents on M. 

7. In the results just outlined, an essential role is played by orientation proper- 
ties ifr >0. Incase r = 0, the theory takes on a simpler form and is applicable to 
other classes of spaces. For instance, if @ is a locally compact group satisfying the 
first axiom of countability, then we can define an invariant metric on G; formulas 
for \, can be given as in Section 1 (condition b is not needed), replacing vector 
fields by group elements and setting |p o = lforallpeG. Wecan use the known 
structure of G to show that \,, isa norm. In this way we can develop a theory of 
distributions on G similar to, but more general than, that of J. Riss. The theory 
does not require the explicit notion of differentiation. 





1 Due to de Rham (see G. de Rham, Variétés différentiables [Paris: Hermann & Cie, 1955]). 
2H. Whitney, Geometric Integration Theory (in preparation ). 

3H. Whitney, Prec. Intern. Congr. Math., 1, 245-256, 1950. 

4L. Schwartz, Théorie des distributions, Vols. I and II (Paris: Hermann & Cie, 1950-1951). 

5 H. Hopf, Ann. Math., 42, 22-52, 1941. 

6 J, Riss, Acta Math., 89, 45-105, 1953. 


GENERAL FACTORS OF AUTOMORPHY 
By R. C. GuNNING 
PRINCETON UNIVERSITY 


Communicated by S. Bochner 


Let D be a simply connected analytic manifold, and let T be a group of analytic 
homeomorphisms of D onto itself, transformations with fixed points being permitted, 
such that (1) the orbit space D/T is compact; (2) each point z e D has a co-ordinate 
neighborhood U such that TU a U + ¢for any T «T implies TU = U; (8) for each 
point z « D the set Tz has no limit points in D; for each pair of points 2, 22 € D 
either 'z, = I'ze, or there exist neighborhoods U;, U2 of 21, 2. such that TU; an TU, = 
¢; (4) there exists a Kaehler metric on D which is invariant under the group I. 

Definition: A set of factors of automorphy for the group T on D, as defined by 
Bochner,' is a collection of functions | vr(z)} which are holomorphic on D, are in- 
dexed by the elements 7 « I, and satisfy vsr(z) = vs(Tz)vr(z). We say that two 
sets of factors of automorphy { ur(z)} and { vr(z)} are equivalent if there exist com- 
plex constants er and a holomorphic nonvanishing function h(z) on D such that 
ur(z) = ervr(z) h(Tz)/h(z2). 

Consider any presentation + F/R, where F is a free group on a set of gener- 
ators {4,,..., ¢:}. For each ¢; select an arbitrary branch of the logarithms to de- 
fine the single-valued holomorphic function o;,(z) = log vr,(z), t; corresponding to 
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the transformation 7’; of !. Extend these functions to a set | ,(z)} indexed by all 
t « F and satisfying o,,(z) = o,(Tz) + o,(z). Then, for any t e R, ¢(t) = 1/2x7 
o,(Z) is an integer constant; this defines a homomorphism ¢:R — Z, where Z is the 
additive group of the integers. The restriction of this homomorphism to the sub- 
group R n [F, F] of RF is called the character class of the factors {| vr(z)}; it is 
independent of the choice of the functions o7,(2). 

THEOREM |. T'wo sets of factors of automorphy are equivalent if and only if they 
have the same character class. 

THEOREM 2. Jf IT’ contains no transformations with fixed points and Cousin’s 
second problem is solvable on D, then, whenever a meromorphic function f(z) on D 
satisfies {(Tz) = vr(z)f(z), the divisor of f(z) considered as a two-cocycle of D/T belongs 
to the cohomology class determined by the character class of the factors | vr(z)} under the 
Hopf isomorphism? R n [F, F\/(R, fF] = H2(D/T; Z). (Remark: This theorem 
can be adjusted to include a result for groups admitting transformations with 
fixed points.) 

Applying Theorem | when D is the entire complex affine space and [ a group of 
translations, the following theorem of Appell* results: 

THEOREM 3. In the complex mulliperiodic case all factors of automorphy are of the 
form vr(z) = exp (> ar’2, + br)h(Tz)/h(z), where ar” and br are complex constants 


and h(z) is a nonvanishing entire function. 

On the other hand, when D is as in Poincaré’s work, the following theorem re- 
sults: 

THEOREM 4. Jf D is a bounded subdomain of the complex affine space and 
all closed T-invariant differential forms of type (1, 1) on D with integer periods on the 
two-cycles of Theorem 2 are multiples of the Ricci curvature form (which is always 
true in one complex variable, for instance), then all factors of automorphy are of the 
form vr(z) = er (J r(z))™™ h(Tz)/h(z), where cr are complex constants, J7(z) are the 
complex Jacobian determinants of the transformations T, and h(z) is holomorphic and 
nonvanishing in D. 

These fractional powers of the Jacobians in one complex variable have been 
studied by H. Petersson,‘ and this theorem shows that in a sense they are the only 
factors of automorphy that need to be envisaged in many contexts. 

As a generalization of the previous theory, the concept of a set of matricial factors 
_ of automorphy can also be introduced; these are matrices {M7(z)} of holomorphic 
functions indexed by elements 7’ « T and satisfying Mgr(z) = Ms(Tz)Mr(z). Two 
such sets of factors |Mr(z)} and |Nr(z)} are called equivalent if there is a non- 
singular holomorphic matrix function H(z) on D such that H(Tz) = Mp(z) H(z) 
Nr~‘(z).. The simplest results are obtained in one complex variable. A set of ma- 
tricial factors | M7(z)} is called irreducible if there exists no nonsingular holomorphic 
matrix function F(z) such that the factors F(7'2z)M7(z)F-'(z) are of the form 


feds prt 
0 Cr(z) ‘ 


THEOREM 5. In one variable all irreducible matricial factors of automorphy are 
equivalent to factors of the form Crmr(z), where Cr are constant matrices and mr(z) are 
scalar factors of automorphy. 
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INTRODUCTION AND THE DIGITALIZATION 


The formal theory of symbolic logic, as developed by Boole,! Peirce,? Jevons,’ 
and Schréder‘ before 1900 and by Whitehead, Russell, Hilbert, Gédel, and others 
after the turn of the century, is universally recognized for its fundamental im- 
portance as a cornerstone of scientific method and thought. Yet it has remained 
largely a philosophic and esoteric realm of study, and there have existed no methods 
of actual straightforward computation for direct large-scale applications to specific 
types of realistic problems. The purpose of this paper is not to present any new 
results in the formal theory of symbolic logic, but rather to offer a new system of 
digitalized computational methods in the propositional calculus for application to 
the great number of practical nonnumerical problems that so frequently occur in 
science, industry, and government. The object was to formulate a logical 
“arithmetic” of extreme simplicity, which would provide, for this realm of non- 
numerical problems, systematic methods of solution as simple, straightforward, 
and versatile as those of numerical analysis are for problems of a numerical nature. 
The theory of Boolean equations’ becomes a special instance of the more general 
methods presented in this paper. 

It is perhaps not generally realized how wide a range of problems can be attacked 
and solved by methods of symbolic logic and the related Boolean algebra, which is 
the calculus of sets and classes. Direct application of logic can always be an aid 
to deductive reasoning—such as determining consequences of given premises, rules, 
or axioms and making hypotheses or theorems from which the given premises or 
factual relationships can be deduced. Besides the use of logical propositional 
methods in problems concerned with sentences, such as analysis of military in- 
formation reports and legal and insurance documents, there appear to be even more 
important applications to fields of operations research, biology, medicine, design 
of experiments, etc., where the utility of symbolic logic per se is not as immediately 
evident. Now that the results of the propositional calculus of symbolic logic are 
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made computationally feasible by these new digitalized techniques, it is hoped that 
they can be directed toward successful solution of many such problems—two 
examples of problems in biochemistry being given at the end of this paper. The 
digitalization of the methods in terms of the binary number system provides easy 
mechanization on existing high-speed electronic computers or on a special electronic 
digital logic machine. 

This paper will contain no proofs, although the foundations of the methods are 
briefly indicated. Also, it is assumed that the reader is familiar with the methods 
of symbolic logic’ and Boolean algebra’ as given in the elementary reference texts. 


REVIEW OF DEFINITIONS AND NOTATION 

Propositions are sentences that can be called either true or false and are sym- 
bolized by A, or X,. Propositions can be combined by the propositional operations 
yr 4 and “~”, where “A,-A,” (read “A, and A,’’) is true when both A, and 
A, are true; “A, + A,” (read A, or A,) is the inclusive or; and “A,” is the negation 
of A, (read ‘“‘not A,”’). (The analogous set or class interpretation of the symbols 
as intersection, union, and complementation, respectively, is easily recognized.) 
An uncombined proposition is called an elementary element, while a combination of 
propositions by the propositional operations is called a combined element or Boolean 
function, this latter sometimes being denoted by f,(Ai, A», ...). 

It is important to distinguish between two different meanings of truth. Some 
combined elements are true due to the form of the combination, independent of 
the truth values of the component parts, as, for example, X, + X, and (X,-X,) + 
(X, + X,)-(X, + X,). Such elements are called tautologically true, or simply 
tautologies, and are denoted by J. (The intrinsically false element is denoted by 
O. Hence I = 0.) On the other hand, a certain combined element, not a tautol- 
ogy, may be called true under the circumstances of a problem, this being called 
factual (sometimes empirical) truth. (We will not distinguish between the language 
and the metalanguage but will leave the distinction to be ascertained from the 
context.) 

Two combined elements that are emphasized in all the original as well as in the 
modern works in logic are X,-X, + X,-X, and X, + X,. The former is read 
“X, is equivalent to X,’’ (in truth value), abbreviated by X, = X,; and when 
X, = X,, then f(X,, Ai, Ao, ...) = f(Xs, Ai, Ao, ...) for any f. The latter is 
read “X, implies X,” or “if X, then X,”’, abbreviated by X, ~ X,. Note that 
this combination is true if X, and X, are both false, which might seem contrary 
to usual usage in language; but the proposition “If you can do that trick, then 
I’ll eat my hat” actually has this meaning. In general, the difficulties in formulat- 
ing a problem in symbolic logic are no more formidable than the precisely analogous 
task in mathematics and depend primarily on how well defined the problem is to 
begin with. However, only the computational aspects of already formulated 
problems are considered in this paper. 





THE BASIS AND THE DIGITALIZATION 
To every proposition A, will be associated a binary number #4, called the 
designation number of A,; in particular, this holds for the elementary elements. 
The set of designation numbers for the elementary elements is called a basis, and 
one such basis for a system of three elementary elements A;, A», A; is: 
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87654321 
#A, = 01010101 
#A. = 00110011 
#As = 00001111 


where the small numbers above the columns number the positions of the bits. In 
general, it can be shown that in a system of n elementary elements the designation 
numbers will have 2” positions. A basis is distinguished by the fact that the 


columns represent the 2” possible combinations of 0, 1 taken n at a time. Thus 


there are 2"! bases; the basis shown will be used throughout this paper due merely 
to its visual simplicity of pattern. The relation of a basis to the familiar truth 
tables is evident. It will be convenient to denote the bit in the 7th column and 
rth row of a basis by b[A;, ..., Ay]!. For example, in our basis b[A;, A», A3]3 = 1. 

The propositional operations are interpreted in terms of the designation numbers 
in a manner reminiscent of vector addition as follows, where #A! represents the bit 
0 or 1 in the 7th position of #A,: 


#(A, + A,)t = fAL+ #A' according to the rule 1 +1=1+0=0+41 = 1, 
0+0=0; 
#(A,-A,)! = #A'-#A' according to the rule 1-0 = 0-1 = 0-0 = 0, 1-1 = 1; 


and #A‘ = 0, 1 when #A‘ = 1, 0, respectively. 


The operation ‘“‘+”’ is called logical addition, ‘‘+”’ logical multiplication, and ‘~” 
inversion. Thus the designation number of a combined element with respect to 
a particular basis is obtained by performing the operations indicated by the combi- 
nation. For example, with respect to the basis given above, #(A; + As-A;3) = 
10101010 + (00110011)-(00001111) = 10101010 + 00000011 = 10101011. The 
set of all the 22° possible designation numbers forms a representation of a free 
Boolean algebra. (Note that for elementary elements, #4! = b[Ai, ..., A;]}, of 
course. ) 

For the reverse process—namely, given the designation number, to find its 
propositional representation—note that every product of all the elementary 
elements or their negations, called an elementary product, has a single unit, i.e., 
#(A,+Ao+A;) = 00000001, 4(A;+A2-As3) = 00000010, etc., for each of the 2% such 
products. Hence a propositional representation called the disjunctive normal 
form is obtained by forming the sum of those elementary products corresponding 
to the units of the designation number. For example, 01100010 = #(Ai-A2-A; + 
A,+ AoA; + A;+A2:A;). Other propositional representations can be systematically 
generated, such as the simplest-sum-of-products form, simplest-product-of-sums 
form, and the conjunctive normal form.*® 

Important observations are that #/ consists only of units, i.e., #7‘ = 1 for all 7. 
Hence 40‘ = O for all 7. Also, Xi = X°! if and only if #X, = #X, for all 7. In 
addition, X, > X, if and only if #X, has units in at least the same positions as does 
#X,; ie., if 4X! = 0, then 4X‘ = 0, 1, but if #X! = 1, then 4X} = 1. 


COMPUTATIONAL METHODS 


CONSTRAINTS 
The concept of constraints is fundamental for the computational methods. 
Constraints are logical relations between the elementary elements which are to be 
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considered as tautologically true, due to intrinsic or stated circumstances of a 
problem. In other words, constraints are factually true combined elements 
which are to remain true throughout a problem and therefore, with respect to this 
problem, can be considered as tautologies. A systematic computational method 
will be given that automatically guarantees the maintenance of any desired con- 
straints throughout the solution of a problem. After constraints have been 
applied, a constrained propositional calculus results, which can be interpreted in 
terms of a nonfree Boolean algebra. 

In particular, if H is an automorphism of a Boolean algebra B, such that ‘+’’, 
“.”’ and “~” are preserved [i.e., for X,, X,, «eB, H(X, + X,) = H(X,) + A(X;), 
H(X,-X,) = H(X,)-H(X,), and H(A,) = H(A,)], and if F(A,, As, ...) € B repre- 
sents a constraint, then we desire that H(F) = /. Under this automorphism, 
equivalence classes are formed by the set |F} together with the sets | X},, where 
(F + X,)«€ \ F} for all X, « B, this being the kernel of H, and for X, ¢ { F}, X,€ {XxX} , 
if H(X,-F) = X;. The equivalence classes form the elements of the nonfree 
Boolean quotient algebra B/F, when the operations between the classes are defined 
by {X}, + {X}, = {X, + X,}, {X},-{X}. = {X,-X,}, and {X,} = {X,}. If 
#F has u units, then there will be 2°"-” distinct elements in each equivalence 
class and 2" equivalence classes all together. This quotient algebra automatically 
satisfies the constraints and represents a constrained propositional calculus. 

Given the constraints F;, F2, ..., F,, the computational method consists in first 





w=PT 

forming II (#F,) and then reducing the reference basis to include only those 
w=! “wer 

columns j for which II (#F}) = 1 (i.e., those columns corresponding to the unit 
1 


“ 
positions of the product of the constraints). All computations with respect to 
this reduced basis will automatically have 4F, = #7 (u = 1, 2, ..., r), as desired. 
For example, the constrained basis b,[A:, A2, A;] that will make A;, As, A; mutually 
exclusive and exhaustive as is often desired in problems concerning classification of 
subject matter, namely, A;-A2 = O, Ai: A3 = O, Az? As = O, Ay + Ao + Az = 7, 
becomes 


764 
#A, = 100 
#As = 010 
#A; = 001 


For this constraint, A;, As, A; are called the components of the three-component 
proposition A and form an example of an n-component proposition. 
ANTECEDENCE AND CONSEQUENCE SOLUTIONS 

The majority of logical problems involve given or accepted premises, hypotheses, 
rules, or other logical relationships which are the given equations and essentially 
comprise the statement of the problem. There are two types of solutions to a set 
of given equations, the antecedence solutions and the consequence solutions. Ante- 
cedence solutions are hypotheses or theories from which the given equations can 
be deduced; consequence solutions can be deduced from the given equations. In 
other words, the truth of the antecedence solutions is sufficient for the truth of the 
given equations, but the truth of the consequence solutions is necessary for the 
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truth of the given equations. If the given equations are true, then the consequent 
solutions are true, while the antecedence solutions may or may not be true; but 
the truth of the given equations can be deduced from the hypotheses embodied in 
the antecedence solutions, this latter being the method for theory construction. 

However, to produce just one antecedence or one consequence solution to given 
equations is usually trivial; hence logical problems always require solutions of a 
specified form or solutions involving only certain specified elementary elements or 
both. Occasionally solutions with the required properties do not exist, and the 
problem is extended to determine under what conditions such solutions do exist 
for the given equations. (It is of historical interest to note that the problem of 
solution to equations as posed by Schréder and others is merely a special case of the 
general theory presented here, being antecedence solutions of the particular form 

¢ -, fs.) 

More specifically, a complete set of antecedence solutions for given equations 
is any set of constraints F that makes each of the equations true in the constrained 
system B/F (i.e., that maps each equation into J). All combined elements that 
become true (are mapped into J) in the constrained system B/E, the constraints 
being the given set of equations EL, are consequence solutions. With these defini- 
tions in mind, the problem is summarized in Figure 1. The given equations /), 

., £, are combined elements or functions of the elementary elements A,, 

A;, X,, ..., Xx. The specified forms F;, ..., F, for the antecedence or consequent 


Condition for 








Condition for Consequence Solutions 
anes Solutions la A } 
‘A ad cE 
ANTECEDENCE GIVEN CONSEQUEN 
SOLUTIONS EQUATIONS SOLUTIONS 
of forms | of forms 
| Fatty Bie M yee My | yA gate A ee MD Fyltyere so Kye My) 
| F>( ) EA( ) F,( ) } 
| 
where #, = 6,(Ay,°°°.A)) 
Fia. 1 
solutions are, in general, known functions of X;, ..., Xx and of fi, ..., f;, where 
the f; are as yet unknown functions of A,, ..., Aj. The object is to ditaruine 


explicitly the functions f; or, if they do not exist, conditions for their existence. 
Summarizing, then, computational methods for determining the existence of 
antecedence or consequence solutions of specified forms and, if they exist, for 
generating all such solutions, or, if not, for determining conditions of existence, 
will be presented. 

The fundamental formulas for determining antecedence and consequence solu- 
tions are the following: 

ANTECEDENCE CONSEQUENCE 


(Fy) © (Ex) = (Ryd) (Px) ® (Ex) = (Ry) 





where (Fy), (£,;). (Rj;) are Boolean matrices (of elements 0, 1 only), and @ repre- 
sents logical matrix multiplication [i.e., if (Pix) ® (Qe) = (S;), then u PrQE = 


S;, with logical multiplication and addition]. The “~” represents inversion of 
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each individual element of the matrix. The formation of (/;,) and (F,), and the 
interpretation of (R;;) as the desired solutions, proceed as follows. 

The Formation of Matrices (E,;) and (F ,).—Two bases enter into this calcula- 
tion namely, b[A,, ..., A;, Xi, ..., Xx] and b[fi, ..., fj, Xi, ..., Xx], where the f;’s 
are considered as elementary elements and both bases have the usual pattern, 

n=r 
the last rows being formed by the X,’s. Then form II (##,) with respect to 
v=s w=l 
b[A,, ..., Ai, Xi, ..., Xx] and Tl (#F,) with respect to b[f,, ..., fj, Xi, ..., 
y= 1 uw=r 


X,]. Now separate the positions of the designation number II (#£,) into 


m » 
2* successive groups of positions, with 2' positions in each group. Index the 
groups from right to left by k (k = 1, ..., 2") and the positions in each group by 
?(¢=1,..., 2"). The rows of the matrix (/;,;) are simply these respective groups 


of positions. Similarly, I (#F,) is separated into 2" groups, this time of 2! 
y=] 

positions, the groups being indexed by / and the positions in each group by 7 (7 = 1, 

..., 2’). The columns of the matrix (7) are simply these respective groups of 


positions. 
In more mathematical terms, the formula for an element F,; of the matrix (E;,;) 
a @Ff : “wr 
isE,, = TT (##®) for (p — 1) = (k — 1)2'+ ( — 1). Similarly, (Fy) = 
uw=1 : u=r w=! 
(#F%) for (q — 1) = (k — 1)2?+ (j — 1). (Note, of course, that (FE?) = 
= vy=s8 n= 1 


uw=r 
#( Il £,)? and similarly for TT (#¥F%)). 

p= 1 »=1 

The Interpretation of (R;;)—Two bases are consulted in this process, namely, 
b[fi, ..., fj] and 6[A,, ..., Aj], both written in the usual pattern. The desired 
solution is computed by means of the result array, which consists of an as yet 
empty array with 2' columns and j rows; the columns are indexed by 7 from right 
to left, the rows corresponding to fi, ..., f; assigned from top to bottom. Now, 
having computed (R;,;) by the proper formula, consider those pairs of indices J, 2, 
for which R;,; = 1: place the jth column of b[f,, ..., f;] in the 7th column of the 
result array. The rows of the result array thus filled are the designation numbers 
of the corresponding f, ..., f; with respect to the basis b[A,, ..., A;] and hence 
give the desired explicit functions fi(Ai, ..., Aj, ..., fj(Ai, ..., Aj. 

In more mathematical terms, #f,,(Ai, ..., Ai)’ = b[f, ..., fj], for those 7 for 
which R;, = 1, and m = 1, ..., j, where #f,,(Ai, ..., A;) is interpreted with respect 
to b[A;, ..., Ail. 

Multiple Sets of Solutions.—The occurrence of several columns of b[fi, ..., fi] 
in a single result-array column indicates multiple solutions, for each such column 
determines a different set of solutions. Hence the total number of different sets 
of solutions is the product of the number of columns of 6[f,, ..., fj] in each array 
column. If a column of the result array has no basis column in it, then no solutions 
exist at all. 

Conditions for the Existence of Solutions.—Consider the row vector (V;) com- 
prising 2! units. If no solution exists, compute the row vector (C;) from the formula 
(V;) ® (Rj) = (C). Then (C;) is the designation number of the desired condition 
for the existence of solutions, as referred to 6[A;, ..., Aj]. In order to impose 
this condition new bases constrained by the condition must now replace b[A,, .... 
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A;, Xi, ..., X,] and b[A,, ..., Aj], and the solution is computed exactly as above, 
except that for ‘2’ now read “‘u’’, where u is the number of units in (C,). 

A Simple Example.—Consider the especially simple equation of historical interest 
discussed at length by Schréder, namely, A;X + AsX = O, to be solved for ante- 
cedence solutions of the form X = f(A;, As). The computations proceed as follows: 


b{f, X]: b[A,, As, X]: 
#f = 0101 #A, = 01010101 
#X = 0011 #A. = 00110011 
#f = X) = 1001 X= 00001111 
j 2121 #(A\X + AX = 0) = 11001010 
} sole ie 


whence (F.) = (6?) and (By;) = (j6%). Substituting in the antecedence formulas, 


(Fn) @ (Ex) = (32) @ (O31) = (298) = (Ry), whence (Ry) = (188). To 
interpret (2;,), 
b[f]: Result array: 
ee i 4 |3/2 1 | 
4f = 01 wlot 10} 1| | none | 


Thus no solution exists. The condition for the existence of a solution is computed 
by (V;) @ (Rj) = (2 ® (090) = (1110) to be interpreted with respect to 
b[ Ay, Ao]: #41 = 0101, $A. = 0011, whence 1110 = #(A; + AQ), i.e., Ai + Ay = / 
in order that A:X + AX = O have a solution. 

The constrained basis b[A,, Az, X] becomes #A; = 010 010, #42 = 001 001, 
#X = 000 111, whence #(A,\X + AX = O) = 110101 and (E,,) = (102). Then 
#(R;;) = (499), and the result array is 

i] 3 |2)1 


#\0,1/0/1 


which, interpreted in terms of the constrained bases b,[Ai, Az]: #41 = 010, #A2 = 
001, is #f = 001 = 4A. and 4#f = 101 = #A;. Hence, under the proper conditions, 
the given equation has two solutions, X = A. and X = A). 

Finding Antecedence and Consequence Solutions Involving Only Specified Elementary 
Elements.—Let A;, ..., A; be the specified elementary elements, all the rest being 
included in the X;, ..., Xy. Then all the computations proceed as above. For 
example, consider the problem of determining consequence solutions of A,X + 

AX = O involving only A; and A». As above, (E;;) = (jo90), but the desired 
form is simply f, and so, since #f = 0101 as referred to bases b[f, X], (F.) = (92). 
Thus (Fx) @ (Exi) = (60 )) ® (jor0 = (d000 = (Rj). Therefore, (Rj:) = (Trt) 
whence 
Peay 1] 4 j2j2) 
| 


#f = 01 #1}i}1{o,1| 


or f = A, + Az,andf =J. In other words, A; + As = J or, equivalently, Ai- A» = 
O is the desired nontrivial solution. 

Special Theorems.—Solutions of the form f,, = X,,,m = 1, ..., k are particularly 
important and oceur frequently. For antecedence solutions of this form, (Fx) 
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is the identity matrix, and hence (E;;) = (R;;) (as was seen in the above example), 
enabling (R;;) to be written down directly. In addition, it can be shown that if, 
to a set of given equations, both antecedence solutions and consequence solutions 
of this form exist, then there is only one such set of solutions (i.e., no multiple 
solutions), and it is the same for both cases. Also, if only one consequence (ante- 
cedence) solution of this form exists, then it is also an antecedence (consequence) 
solution of this form. 


LOGICAL DEPENDENCE (CONDITIONS ) AND LOGICAL INDEPENDENCE 

Conditions on the Solutions.—Often it is desired to find only those solutions 
fi, ..., fj which satisfy given conditions or constraints G(fi, ..., f;) or, in other 
words, only those fi, ..., f; for which the given combined element or function 
G(fi, ...,f;) = I. The computational procedure is similar to calculating solutions 
with conditions on the A,, ..., Aj, except that now bases constrained by G(fi, 
..., fj) replace b[fi, ..., fj, Xi, ..., Xx], and b[fi, ..., fj] (where #G(fi, ..., fi) 
is determined with respect to b[fi, ..., fj]); and for “2” now read “w’’, where 
u is the number of units in 4G(fi, ..., fj). All sets of solutions thus calculated will 
automatically satisfy the given conditions. 

Logically Independent Solutions.—Since f; = fi(Ai, ..., Ai), ..., fj = fi(As, --., 
A;), it can happen that these equations imply a logical relation between fi, . . ., fj. 
For example, if f; = Ai-Az2 and fe = A; + Ag, then fi > fe. However, often a 
prime requirement of sets of solutions for a problem is that they be logically inde- 
pendent, that is, no relations of the form G(fi, ..., f;) = J exist other than tautol- 
ogies. Then fi, ..., f; will be logically independent when the result array has at 
least 2' different columns. Three cases arise: j = i,j <i,j >i. Ifj = i, then 
there can be no repeated.columns in the result array for independent solutions. 
For j < i, 2! different columns must appear in the array for independence. For 
j > i there can be no independent set of solutions. 

Determining the Logical Dependence between a Given Set of Solutions.—Given a 
set of solutions, it is often desired to find the logical dependence between them, 
if any exists at all. The computational method is to consider the result array for 


this given set of solutions fi, ..., fj. Rewrite the array, omitting any duplicate 
columns, and adjoin to this (on the right) those columns necessary to make the 
completed array a basis, say b’[f,, ..., fj], not necessarily of the usual pattern. 


Then the designation number with zeros in the positions corresponding to the 
adjoined columns, units elsewhere, produces the desired conditions between fi, 
_., fj when referred to b’[fi, ..., fj]. 


EXAMPLES IN BIOCHEMISTRY 


ENZYME CHEMISTRY 

In biochemistry an enzyme often cannot be isolated easily, and therefore the 
experiments performed often will involve several other enzymes. Consequently 
a combination of reactions must be observed. In such a complex situation the 
ordinary simple logical analysis usually used in experimental science is found 
inadequate. In these cases the logical computational methods presented above 
‘an be extremely useful for evaluating experimental results as well as for planning 
future experiments to yield the maximum information, 
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Suppose that a chemist were studying enzymes, A, A», A; in relation to reactions 
X,, X2, X;. Suppose that he has done the following four experiments: (£,)—A 
solution containing none of A;, As, A; produced reaction X» but not X;, and not X3. 
(£.)—The solution contained A; and either A» or A3, or both; he could not be 
sure. The reaction was neither X, nor X; and X; together. (#;)—The solution 
had A» but not A;, or did not have A» but had A;. Reactions X,; and X» occurred, 
or reaction X; occurred but X, did not. (£,)—The solution was obtained from a 
source that had A; and A, or A, or both, or else had neither A; nor A;; the solution 
turns color, which means X, does not take place or both X, and X; do. 

The antecedence problem is: What theories about the enzymes associated with 
each reaction will explain the experimental results? The consequence problem is: 
What combinations of enzymes are necessary for each reaction to take place? 

The experiments and their designation numbers with respect to the usual basis 
are as follows: 

EF: A,:A2:A; = Xe kek: 
O111 1111 O111 1111 1000 0000 O111 1111 O11] 1111 OLL1 1111 O1L1 1111 O111 1111 
E,: Ay+(Ay + As) = Xo*(X1°X3) 


0001 0101 0001 0101 1110 1010 1110 1010 0001 0101 1110 1010 1110 1010 1110 1010 


E;: Ay: A, + Ay: Az = X,-X.+ X,°X3 


1101 0001 1101 0001 1101 0001 0010 1110 0010 1110 1101 0001 0010 1110 0010 1110 


Ex: A3*(A1 + Ae) + A, As - Xi + X2°X; 


1010 0111 0101 1000 1010 0111 0101 1000 1010 0111 0101 1000 1010 0111 1010 0111 


aw=4 
Thus I (#h,) = 0000 0001 0001 0000 1000 0000 0000 1000 0000 0100 0100 0000 0010 0010 0010 0010 

w= 1 
For the antecedence problem the form of the solution is X, = fi, Xx = fo, X3 = fs. 
tecall for antecedence solutions of this special form that (£,,;) = (R,,). 


Hence 

8765 4321 
0000 0001 
0001 0000 
1000 0000 
0000 1000 


0000 0100 
0100 0000 
0010 0010 
0010 0010 


>I OO 


(Axi) = (Rj) = 


mh e oF 


from which the solutions are read: 


bi fi, fo, fs}: Result array: 
j 87654321 Ne) ee ee 6 5 | 4/3 2 1 
#f; = 01010101 #1.}0}1)1,0)/1)1/0]1,0/0 
#f2 = 00110011 #11/0/1,1/0/1/0]1,11]-0| 
#f; = 00001111 #e10])1/1,1}0}]0]1)] 1,1 | 0 
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Hence there are four sets of antecedence solutions of the form X, = f,, namely, 
A, A; + A;-As, A3*(A + A2) + A,: Az As, 
A\ ° A\ ’ 

Ay: Ag + Ay:A 


X, = fi, where f; 
X>2 = fo, where fo 


Xe = fs, where fs = A, Ae + Ay: Ae, 


II 
II 


As(Ai + A») + A, As, AyA; + A: A»: As; 
A ’ Aj; 
A, Ae os A, Ae ; A, As ++ A,: Ao. 


The first two sets are logically independent. The last two sets of solutions imply 

that the reactions are not independent but that X.-X; —~ Xi, and X,-X2-X; = O, 

respectively. For example, we write the designation numbers of the last solution, 

cross off the repeating column, adjoin the missing column and obtain: 
#(A\+A; + A,+Ay+A;3) = 0101 10601, 
#4, = 1010 104031, 
#(Ai+ Ay + Ay As) 0110 01%0i1, 
1111 1110= #(firferfs = 0), or X,+Xq:X3 = O. 


For the consequence problem, the form of the solution and corresponding designa- 
tion numbers with respect to b[fi, fo, fs, X1, X2, X3] are as follows: 


#(X, — f)) 1111 1111 0101 0101 1111 1111 O101 O101 1111 1111 0101 0101 1111 1111 0101 0101 
#(X2 — fo) 1111 1111 1111 1111 OO11 OO11 OO11 OO11 1111 1111 1111 1111 0011 0011 0011 0011 
#(X; — fs) 1111 1111 1111 1111 1111 1111 1111 1111 0000 1111 0000 1111 0000 1111 0000 1111 
y= 4 
I] (#F,,) = 1111 1111 0101 0101 0011 0011 0001 0001 0000 1111 0000 0101 0000 0011 0000 0001 
y» = ! 
hence 
1000 0000 0000 0001 
1100 0000 0001 0001 
1010 0000 1000 0001 
: 1111 0000 ak : = 1001 1001 
€ 7 - = F y 2 aye) a = 
(Fx) ioe’ “shi and (F,) ® (E;:) (R;,), whence (R;;) aoe aaen 
1100 1100 0101 0101 
1010 1010 1000 0101 
1111 1111 1111 1111 


whence the result array: 
ee ee 7 


Ce Se 6 5 x a ee ee 
#f, | 1010 11 | 1 | 


1111 | 11 | 1010 | 1 
1010 | 11 | 1100 1 | 
1100 | 10 | 1111 | 1 


#f. | 1111 10 
#f, | 1100 | 11 


Since the most significant set of solutions for this form is the minimum set, i.e., the 
solutions with the most zeros, we choose 


#f, = O111 1010 = #[A3-(A, + Av) + AiAs], Xy— Ase (Ay + Ao) + Are As; 
#f2 = 1010 1010 = #A,, or X,— A; 
#f = 01 10 01 10 #[A,- Ay oe Ay Ago], As ae A,:Ae a AAs. 
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THE PROTEIN DECODING PROBLEM 

According to present views,’ hereditary informations are stored in chromosomes 
in the form of long sequences of four different units (known as_ bases) 
constituting the molecules of nucleic acid. In the process of protein synthesis, 
these polynucleotide sequences must be translated into the corresponding poly- 
peptide sequences of twenty different amino acids. Thus there must exist a 
mechanism which carries through such translations in a unique way, producing 
a long word written in a twenty-letter alphabet for each long number written in 
the four-digital system. Since twenty different triplets can be formed (with 
disregard to order) out of four different elements, it is inviting to associate each 
amino acid in the resulting polypeptide chain with a group of three bases in the 
original polynucleotide sequence. One can make different assumptions as to the 
type of correlation involved. For example, we can assume that while one amino 
acid is defined by three neighboring bases, its neighbor is given by the next three 
bases: 


Or one can make the assumption of overlapping triplets, in which case each base is 
participating in the choice of more than one amino acid. For example, we may 
have the scheme 

—1—3—4-—-3—2—3 -38— 


(eeepc 


u—_—" 


Oe oe 


or the more restrictive scheme 


The possibility is also not excluded that because of twisting of the chain, non- 
neighboring bases participate in the choice, as in the following: 


oe a a a 
eS 


E wN M 





In all overlapping translation codes, one should expect a stronger or weaker inter- 
symbol correlation, i.e., some restriction in the choice of neighboring amino acids. 
Thus, comparing these restrictions with actually observed sequences of amino acids 
in protein chains, one should be able to find the correct coding procedure and the 
one-to-one correspondence between the twenty amino acids and the twenty hypo- 
thetical triplets of bases. Since there are 20! = 2.3 X 10" such one-to-one corre- 
spondences (which equals the number of seconds in the age of the universe), the 
straightforward test of all possibilities is out of the question, and it is desirable to 
find a systematic and feasible procedure for solving the problem (with the aid of 
high-speed electronic computing machines). 

Use of the Logical Computational Methods Developed Above.—To illustrate the 
method of solution, it is best to consider a much simpler case which makes more 
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transparent the computational methods involved and can easily be extended to 
the complex problem given above. Suppose that there were only three kinds of 
amino acids and only three types of base triplets allowed. Let the rules of combina- 
tion involve groups of three consecutive triplets, and let the first position of the group 
be denoted by A, the second by B, and the third by C. In addition, the subscript 
1, 2, or 3 will denote which one of the three allowed triplets is in the position, i.e., 
A, stands for “The first position has triplet 2”. Hence A,, B,, and C,, for s = 1, 
2, 3, are each three component propositions. Similarly, let X,, Y,, Z,, where r = 1, 
2, 3, represent the positions of one of the three amino acids, considered in consecu- 
tive groups of three. Let the hypothetical conditions on the sequences of triplets 
be: 


A, By > C3 + Cr A, By > Ci + C2 
A,-B; > C2 + C3 A,B; > C; + C0; 
Az? By > C, + C3 Ao By > C, + C2 
Ao B, > C, + C; Ao Bo > Ci + Co 
A3-B, > C) + Cy A2B; > C2 + C; 
A;3°B,—> C, + C2 A;-B,—> C2 + C; 
Ay By > Ci + C3 A;:B; > C; + C; 


Suppose that the experimental results gave the following conditions on the amino 
acids: 
X3: Y3 ony Ze oa Zs, Xi: Y; — Zi a Zo, Xo Ys ce Z, + Z>. 

The problem is to find antecedence solutions of the form X, = f,, Y, = g,, and 
Z, = h, for r = 1, 2, 3. Introducing the following shorthand: 1 represents the 
1 0 ¢.. Sey 8 
column 0, 2 the column 1, and 3 the column 9, since A,, B,, C, (and X,, Y,, Z,) are 


three component propositions, find for b[A,, B,, C,| (or b[f,, g,, h-]): 


i (or j): 272625 242322 212019 181716 151413 121110 987 654 321 
#A,(or#f,) =123 123 123 123 123 123 123 123 123 
#B, (or #97) = 111 222 333 111 222 333 111 222 333 
#C, (or #4.) =~ 1111111311 222 82322 222 333 333 333 


However, observe that we really know all possible answers, because there must be 
a one-to-one correspondence between the components of X, and those of A,; and, 
in addition, the same correspondence must hold between Y, and B,, and Z, and C,. 
Hence the 3! possible solutions are as follows: 


Solution number: 1 2 3 t 5 6 

X, = Aj, Ao, Ai, As, As, As 

Xo = Ao, Ai, Az, Ai, As, As 
1;, 


X3 _ A;, é A», Ao, A, A, 
Letting P = 121233, Q = 213312, R = 332121, the result array for the solutions 
with respect to b{A,, B,, C;| can be written as follows: 
a 27 26 25 24 23:22 212019 181716 151413 121110 9 8 7 65 4 321 
#f, PQR PQR PQR PQR PQR PQR PQR PQR PQR 
#9, PPP QQQ RRR PPP QQQ RRR PPP QQQ RRR 
#h, PPP PPP PPP QQQ QQQ QQQ RRR RRR RRR 
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where this is to be interpreted as 3! = 6 solutions (not 6” solutions according to 
the notation used above), where each solution corresponds to the same respective 
column for 7 = 1, ..., 27. Hence (R;;) can be constructed, where we put the 
number of each solution in the proper Rj, element. Thus we obtain Figure 2. 


U 
27 26 25 WA23 16 15 14 13 

(=k) 27,3 
26 
25 
24 
23 
22 
21 


Fig. 2 


However, this matrix (R;,) would be the result with no conditions on A,, B,, 
C;, X,, Y,, Z,. Adding the conditions to b[A,, B,, C;, X,, Y,, Z,] amounts to 
eliminating some columns (corresponding to the zeros of the conditions on A,, 
B,, C,) and some rows (corresponding to zeros of those on X,, Y,, Z,) as indicated 
in Figure 2. Now only those solutions remain that are contributed to by every 
column. Hence, solution 1 is eliminated, for it does not appear in column 19 or 
column 9, and so forth. Only solution 4 is not eliminated and hence, under the 
given and experimental conditions, it is the only valid hypothesis. 

Now we are able to sketch briefly the feasible method of solving this problem 
by means of computers.'° First, there is a method for putting the n! permutations 
of 1, ..., n into one-to-one correspondence with 1, ..., n!, so that the correspond- 























a, LLL EEO —V— — 
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ence can be easily accomplished in either direction. Thus these integers will 
number the solutions uniquely. Second, there is also a systematic procedure for 
determining the solution number, given 7 and 7 (= k). So choose some allowed 
column 7%, and find the solution numbers for all 7 that have not been eliminated 
due to constraints. Solution numbers not found by this step cannot be valid 
solutions, for they do not appear in this column iz. Third, there is a procedure 
for determining j, given 7 and the solution number. So for all other allowed 
columns 7, and for each solution number just obtained, ury to determine 7. Those 
solutions for which no allowed 7 can be found for every allowed 7 cannot be valid 
solutions. The remaining solutions are valid. 

Of course, for 3! = 6 solutions it is easiest to try them all. However, the feasible 
method just presented is intended for choosing one out of the 20! possible solutions 
of the original problem. The (R;,;) matrix in this case is 8,000 X 8,000, but, when 
a few hundred conditions are applied to A,, B,, C,, X,, Y;, Z,, it will reduce to about 
5,000 X 5,000. It is estimated that the complete calculation should take a com- 
puter no more than a hundred hours. On the other hand, to try 20! solutions, 
a computer put to work in the days of the Roman Empire, at a rate of one million 
solutions per second, 24 hours a day, all year round, would not yet be close to 
finishing the job. 


The author gratefully acknowledges the assistance and guidance of Dr. George 
Gamow, particularly with respect to the general description of the Protein Decoding 
Problem. The kind encouragement of Dr. Nicholas M. Smith and the Operations 
Research Office is especially appreciated. 
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tions Research Office and is part of an investigation of methods for applying symbolic logic to 
operations research problems. 
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THE BEHAVIOR OF A COMPLEX FUNCTION AT INFINITY* 
By Kart MENGER 
ILLINOIS INSTITUTE OF TECHNOLOGY 
Communicated by Lars V. Ahlfors, May 11, 1955 


Traditionally, the behavior at © of a complex function f is defined as the be- 
havior at 0 of the function obtained by substituting the — 1st power into f. This 
definition adequately describes the class of values f(z) for large | z| . For instance, 
the range near ~ of the identity function 7 (whose value for any z is j(z) = z) coin- 
cides with the range near 0 of the function 7~'. But that definition does not in any 
way describe the structural behavior of f near «, reflected in properties of the 
class of pairs (z, f(z)) for large | z| . In fact, the association of the value f(z) with 
z may, by the substitution of 7~', completely change its character. For instance, 
the derivative of 7 near © is the constant! function /, whereas that of j7~! near 0 
goes even faster to © than does j~'. 

As a result of their emphasis on the mere functional range near ©, analysts had 
to assign to ~ a role differing from that of other points on the Riemann sphere, R, 
much more profoundly than the arithmetical peculiarity of © seems to warrant. 
In what follows, the function-theoretical discrepancies between © and the other 
points on R will be reduced to a minimum. Inevitable are only differences which 


result from the fact that quotients of the form (f(z) — f(@))/(a — ©) remain 
undefined, which, in turn, is a direct consequence of the arithmetical peculiarity of 
co, expressed in the formulaa + © = o foranya # o, 

Let f be a function whose domain is an open subset of R and whose values belong 
to R. For any point, a, in Dom f or on its boundary (a # © or = ©), I will 


define the exterior derivative of f at a as the element b of R (if there is such an ele- 
ment) which has the property 


lim [(f(2n’) Le f(Zn))/(Sn" sr Zn) | a b 


for any two sequences {z,} and {z,’} of elements 4a and # © for which lim z, = a 
and lim (z,’' — z,) = 0. The number b, which is a complex analogue of the slope 
of a “paratingent’”’ in the sense of Bouligand,? will be denoted by Ef(a). More 
generally, 

E‘f(a) = lim [(E*~"f(an’) a E*~"f(an))/(an’ ail Zn) |, 
where E°f = f. 


If ais in Dom f and #~=, then E*f(a) exists and is equal to the classical k-th- 
order derivative D“f(a) if and only if f is regular at a (7.e., differentiable at a and in 
a neighborhood of a). At ©, one finds, for instance, Ej(@) = 1, Ej?(o) = o, 
E"j"(o) = n! for any positive integer n, and EXj"() = 0 for any negative integer 
n. If f(z) = (az + b)/(ez + d), where ad — be # 0, then Ef(~) = 0 or a/d, 
according to whether c #~ 0 or = 0. Except for their ranges, these functions “‘be- 
have” at ~ as they “behave” at 0. 

If a point a (finite or ~) belongs to Dom f or to its boundary, then it will be called 
a point of extensibility of f if lim, f exists on R. I will say that, at such a point, 
f is steep if E*f(a) = © for almost all k > 0 and flat if E*f(a) = 0 for almost all 
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k > 0 without E*f being constant. Points of extensibility at which f is neither 
steep nor flat may be called ordinary; and f will be said to be ordinary in an open 
set if and only if the latter consists of ordinary points of f. 

If the point a belongs to Dom f or is an isolated boundary point, then f is steep 
at a if and only if a # © and a isa pole of f, whereas f is flat at a if and only if 
a = o and f is not a polynomial in the neighborhood of ~. If Dom f is a wedge 
containing the real axis and possessing a sufficiently pointed cusp at 0, and if f 
assumes the value e~'/* for any z in that wedge, then 0 is a nonisolated boundary 
point at which f is flat. If Dom f is the region between the real axis and a curve 
having this axis as an asymptote of sufficiently high order, and if f assumes the value 
sin \/z for any z in that region, then Ef() = 0, even though is a point of inex- 
tensibility of f. 

In many celebrated theorems, the traditional distinction between © and finite 
points on R, as well as between the interior and the exterior of simple closed curves, 
becomes unnecessary. For instance, Cauchy’s integral theorem can be formulated 
as follows: If f is ordinary in either complementary region of a rectifiable, simple 
closed curve T' and satisfies proper conditions on I, then fy f = 0. (For, if f is 
ordinary in the complementary region containing ©, then f is a polynomial in that 
region.) Similarly, Cauchy’s integral formula follows from 























14 i * Sz) — Fla) ds = 0; briefly, ts f se Ja = 0, 
i Z—a j — ja 
i] for any point a not on land 4 ~. 

i | Suppose that a is an isolated point on the boundary of Dom f, and let T be a 

Ne curve such that the complementary region of I containing a consists, with the 

ie possible exception of a, of ordinary points of f. Define the residue of f at a as 

‘ (i/2mi) f; f, where [ is traversed once in such a way that the domain containing 

i a is on the left side of Tf. Then, regardless of whether a # © or = &, one finds 
Vi Res, f = 2z7ic_1, where c_; is the coefficient of z~! in that Laurent expansion of 
| f(z) about a which is “closest” to a. An expansion about an isolated singularity 
at © that is closest to © is, of course, a series >> c,(z — a)" that converges outside 

\ some circle (punctured at ©) about some center a. 

i A If |f| is bounded on the sphere R, punctured anywhere, then f is constant. 

i If, in a punctured circle about a, |f| has either an upper bound <@ or a lower 

(4 bound >0 and no upper bound, then a is a point of extensibility of f, and the value 
4 i at a of the extension is ¥ ~ in the first case and = ~ in the second. 

t\ ( The ideas here developed can readily be extended to multifunctions.’ 

A | * Part of the work on this paper was done as part of Project TB 2-0001(888), sponsored by the 
| | Office of Ordnance Research. In writing this paper, the author has profited from several sugges- 
Hi | tions made by Professor R. Struble. 
| 


1 In this paper, numbers (0, 1, i; only exception +) are denoted in roman type; functions (e.g., 
the constant function / and j, assuming the values /(z) = 1 and j(z) = z for any 2), in italics. 
Correspondingly, roman letters (z, a, .. . ) are used as numerical variables and f as a function 
variable. (Cf. K. Menger, Calculus: A Modern Approach [Boston, 1955].) 

2 Bouligand, Introduction a la géométrie infiniiésimale directe (Paris, 1932). 

3K. Menger, “A Simple Definition of Analytic Functions and General Multifunctions,” these 
PROCEEDINGS, 40, 819-821, 1954. 








ON THE GENERATION OF THE IRREDUCIBLE REPRESENTATIONS 
OF THE SYMMETRIC GROUP 
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Communicated May 8, 1955 


We point out in the present note that the various irreducible representations of 
the symmetric group S,, on m symbols may all be generated by symmetrizing 
Kronecker powers of one of these irreducible representations. If p —{ 4 mi " ") 
is any element of S,,, we obtain an m-dimensional representation of S,, by associat- 
ing with p the permutation matrix P = {ep,,..., pm}, Where ey, j = 1,..., m, is 
the m X 1 matrix all of whose elements are zero save the jth, which is 1. This is 
the representation furnished by the permutations of the cosets of S,,-; in S,», under 
left translations of S,,, and its character is the number a; of unary cycles in the class 
of S,, at which the character is being evaluated. It is reducible, being the sum 
of the identity representation T.,,, of S,, and of the irreducible representation 
P¢m—1, 1), Whose character is a, — 1, of S,,. We denote T'(,-1, 1 simply by T and 
observe that T @ 10} = Tm) andr ® {1} = T. The Kronecker square  X IT of 
I is the sum of [ @ }2} and Pr @ {1°}; the Kronecker cube Tr X TI X T of T is 
(T @ {3}) + + 2(r o {21} ) + (7 @ }15!); the Kronecker fourth power! XT Xq 
r xX Toflris(r ® {4}) + 377 @ {31}) + 27 @ {24) + 3(1r @ {214) + (FT @ 
114}), and so on (the coefficients 1, 2, 1 in the expression for T X TX TI’ being 
the dimensions of the irreducible representations of S;, and the coefficients 1, 3, 2, 3, 1 
in the expression for l & I X I’ X I being the dimensions of the irreducible repre- 
sentations of S;, and so on). The principal result which we wish to communicate 
in the present note is that © @ {1*}, k = 1, 2,..., m — 1, is the irreducible repre- 
sentation of S,, whose characteristic is }m — - 1*} and whose dimension is the 
binomial coefficient (";'). If we agree that {1°} = {0}, we see that S,, has m 
irreducible representations whose dimensions are the binomial coefficients corre- 
sponding to the exponent m — 1; form = 2 and m = 3 these furnish all the irreduci- 
ble representations of S,,, the dimensions being 1, 1 and 1, 2, 1, respectively. For 

= 4 they furnish four of the five irreducible representations, namely, those whose 
dimensions are 1, 3, 3, 1, the remaining irreducible representation being 2-dimen- 
sional. We indicate that T @ {1*} is the irreducible representation of S,, whose 
characteristic is {m — k, 1*} by the notation T @ {1*} :{1'}. 

It follows readily from the relation Fr @ {1°}:{1°} that Fr @ {2}:{2} + {1} + 
{O}, ie., that [ @ {2} is the sum of the three irreducible representations of S,, 
whose characteristics are {m — 2, 2}, {m — 1, 1}, and {ml}, respectively. Thus 


lin-2,2 = (1 @ 12h) — (Fr ® {1}) — (7 @ {0}), 
l\(m-2, 12) = ing r @ | 17} 


Similarly, it follows from the relation © @ {1*}:{15} that 


Tr @ {21}: 21} + fo} + {12} 4 
r @ {3}:{3} + f2} + (14) +21 ew fo}. 
514 
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From these and the previous results we read that 
lim—a,3) = (I ® {3}) — (Fr ® {2}) -(r @ 117}) — (Fr @ f1}), 
Pima, 21) = (F @ [21j) — (PF ® [2}) — (FP @ 14) + (F ®@ {O}), 
Dims, 13) = F @ {13}. 
Similarly, it follows from the relation T @ {14} :} 14} that 
Pr @ {217}: {212} + {2a} + {13} + {1 
r@ {2+} :{24) + {3} + {21} + 2f2} + {1} + fo}, 
P® $31j 2431} + 13) + 221) + [1 + 2k2} + 3hl} + 2a 
r @ {4}:{4} + {3} + {a1} + 3fat + {12} + 3f1} + 2{0} 
and from these and the previous results we read that 
Pim 4, 4) ro (T ® 14}) —- r @ {3}) r @ {21}) . (r @ {2}) 2% 
(r @ {17}) + (F @ {1)), 
Pim—s.31) = (T ® {31} — (F ®@ {3}) — 277 @ {21}) — (7 @ {14) + 


(Tr @ {2}) + (Fr @ {1}), 
lim—4,2) = (T @ {2%) — (Fr @ {3}) — (Fr @ {21}) + ar @ {1%}) i 
2(r @ {1}), 
U'im—4, 21) = (T @ {21°}) — (F @ {21}) — (r @ {13}) + (1 @ {2}) — 1 @ {O}), 
Pim—4,1) =T @ {14, 


and so on. We call attention to the fact that the right-hand sides of these equa- 
tions are, formally, independent of m; however, m enters them through T = ['yn—1, ». 


-_—_ 


* Present address: 6202 Sycamore Road, Baltimore 12, Md. 


ON THE ANALYSIS OF THE KRONECKER PRODUCT OF IRREDUCIBLE 
REPRESENTATIONS OF 8S, 
By Francis D. MurNAGHAN 
INSTITUTO TECNOLOGICO DE AERONAUTICA, SAO JOSE DOS CAMPOS, BRASIL* 
Communicated May 25, 1955 


We considered, ! some seventeen years ago, the problem of analyzing the Kronecker 
product of any two irreducible representations of S,, the symmetric group on n 
symbols, into its irreducible components and have given formulas furnishing 
I'(n — p, (u)) X T(n — q, (v)), where (u) and (v) are partitions of p and q, respec- 
tively, for p + q < 6 and for p = 3,q = 4. Furthermore, we availed ourselves of 
these formulas to furnish tables giving the analysis of the Kronecker product of 
any two irreducible representations of S, for n < 8. Recent papers by Gamba 
and Radicati? and by Robinson and Taulbee* have again directed our attention to 
this question. If pis any nonnegative integer <n and (u,..., u;) 18 & partition of 
p for which u, < n — p, we denote simply by (yu) the irreducible representation 
I'(n — p, (u)) of S, which is associated with the partition (n — p, (u)) of n, and 
we say that the irreducible representation (u) of S,, is of depth p; thus the identity 
representation (0) is the one and only irreducible representation of S, of depth 0, 
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and (1) = I'(n — 1, 1) is the one and only irreducible representation of S, of depth 
1. If (u) and (vy) are two irreducible representations of S, of depths p and q > p, 
respectively, the greatest depth of the irreducible components of their Kronecker 
product (u) X (v) is p + q, and the least depth of these irreducible components 
>q — p. Denoting, as usual, by {\} the character of the irreducible representa- 
tion of the full linear group which is associated with the partition (A, dx) of 
any nonnegative integer m, |X! {n’} is a linear combination, with nonnegative in- 
tegral coefficients, of characters il. where (A,”, Ao”, .. . ) isa partition of m + m’, 
and we denote by ({A} {\’}) the same linear combination of the irreducible repre- 
sentations (A”), of depth m + m’, of S,. For example, since {1} yt, where (7, 
., %) is a partition of g, = |», + 1, m,..., mit... t{n,...,m +1) 4+ 
v,, 1}, ({1} {v}) is the linear combination (» + 1,..., %) +...+ 
~>% +1) + ( vz, 1) of irreducible representations, of depth q + 1, 

of S,. The first result which we announce in the present note is as follows: 

The irreducible components of maximum depth, namely p + gq, of (u) X (») are 
furnished by the expression (|! {v}). For example, the irreducible components of 
depth 2 of (1) X (1) are (2) + (1°). This result, for the special cases p = 1, 2, 3, 
is implicit in formulas given by Makar. 

If the partition (Aj, ..., A,) of any nonnegative integer m has k nonzero parts, 
{\} may be written as the determinant of a k-dimensional matrix whose elements 
aie the complete symmetric functions p,; of the characteristic numbers of a typical 
element of the full linear group, the diagonal elements of the matrix being py,, .. . , 
Pr. Tf (Y’ \’,,) is any partition of any nonnegative integer m’ into k’ nonzero 
parts, we understand by {\/A’} the determinant of the k-dimensional matrix ob- 
tained by reducing the subscripts of the p’s in the first, second, and so on, columns 
of the matrix whose determinant is {A} by Ay’, deo’, .. ., respectively (it being under- 
stood that {A/A’} = Oif k’ > kh). It follows readily that |\/\’} = 0 if m’ > m or 
if, when m’ < m, Ay’ > Ay; moreover, if \y’ = Ay, {A N= { (Ae, Tr, ee 
d,-’)} and so on. It is easy to see that {A/1} = {A, — 1, ro, ..., Af Hee + 

\, — 1}, and we understand by ({X/1}) the sum of the irreducible represen- 
tations (A; — 1 A, — 1), of depth m — 1, of S,. Similarly, 
we understand by (j\/1} {1}) the sum of irreducible representations, of depth 
m, of S, which is furnished by the expression ({A; — 1,..., Act {1} + ...4+ 

A, — 1} }1}). In this notation the second result which we here announce 
is as follows: 

The irreducible components of depth p + q — 1 of (u) X (v) are furnished by 
the expression ({u/1} }v/1} }1!). For example, the irreducible components, of 
depth g, of (1) X (v) are furnished by the expression ({v/1} }|1}); in particular, 
the irreducible component, of depth 1, of (1) & (1) is (1). 

Our next result is as follows: 

The irreducible components, of depth p + q — 2, of (u) 
by the expression ({u/1} {v/1} + [{u/2} fy/2t + fu/1?} | 
tw/1*} + }u/1*} }v/2}] 114). 

In particular, the irreducible components, of depth g — 1, of (1) X (vy) are fur- 
nished by the expression ({»/1}), since both {1/2} and {1/12} vanish. Since (1) X 
(v) hav no irreducible component of depth <q — 1, this completes the analysis of 
(1) X (») for all irreducible representations (v) of S,: 


xX (v) are furnished 
12+] {2} + [{u/2} 


p/ 1! 
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1+ fv/t} {a} + fap fp). 


bite of (1) X (1) = (0) + (1) + ahi (12);(1) X¥ (22) = (1) + (2)+ 1%) + 
(3) + (21); (1) & (1?) = (1) + (2) + (1?) + (21) + (14), andsoon. Each of these 
is a universal formula valid for every n; a we evaluate (1) (1°), for example, for 
n = 4 the term (21) = I'(n — 3, 21) vanishes, and we obtain the analysis ['(31) X 
r(21?) = F(3l) + (2?) + F(21*) + T(14) 

We content ourselves with giving the formula which furnishes the irreducible 
components of depth p + ¢ — 3 in (uw) X (y) as the formulas become increasingly 
complicated as the depth diminishes. — It is 


({[fm/2} + fu/E?f) Uho/2f + fo/V flat + [ta/8t te/3} + fu/21f fv/21f + 
fu/1®} {v/18}] {3} + 14/3} | {v/21} + {u/21} {»/3} + 
{u/1*} {v/21} + fu/21} fv/13}] {21} + [{u/3} fo/13} + {u/1*} {v/3}] {14})- 


If (u) is (2) or (12), this reduces to ({}v/2} + }v/12!]}1}), since the three expressions 
| u/3}, | 4/21}, }u/1*} vanish. We complete the analysis of |u} X {v} when (u) is (2) 
or (1?) with the observation that the irreducible components of depth g — 2 in (2) X 
(v) are furnished by the expression (| »/2}), and the irreducible components of depth 
q — 2 in (1?) X (y) are furnished by the expression ({»/12!); thus 


(2) X ( ({o/2} + [hv/2} + to/l? fF] Uap + to/lf tafe + 
{»/2} fo} + foie} {13 + {v/1} {1}? + {>} {2}), 


(12) X (v) = ({v/12} + [fo/2} + fot} ] tat + fof} fa} + 
fo/* 42} + fo/2} 10} + fot} Lae + toh Up). 


When (x) is of depth >2, we think it simplest to proceed by a method which we 
illustrate by the example (1*) X (6), whose analysis, in the case n = 15, is given 
by Makar.‘ The terms of depth 9 are (712) + (61%), since {15} {6} = {717} + 
| 615} ; the terms of depth 8 are (71) + (62) Mio ke + (521) " (51%), since | 1*/1} 

fiet, {6/1} = {5}, and {12} {5} {1} = (tel! 1512 teh = {71} + {62} + 
312} + 1521! + {515!; the terms of depth 7 are ae + (52) 3 3(512) + (421) 4. 
3) einee {1' 2} = 0, {19/17} = {1}, {6/2} = {4}, {6/12} = 0, and {12} {5} + 

{4} {et = fort + {512} + {1} [f51} + {412}] = 2f61} + {52} + ist at 4 
2°} . | 415! ; the terms of depth 6 are (6) + 2(51) + (42) + 2(412) + (31°), since 
| 18/3} and }1*/21} are both zero, while {1*/1*} is the character {0} of the aeas 
ere ae of the full linear group and {1} {4} {1} + {3} {15} = fa} [f5} + 
f41f] + {417} + {315} = 6} + 2{51} + {42} + 2}412} + {318}. (1%) x (6) 
does not contain any irreducible representation of depth 3, since none of the products 
ret {3h fast | fort / 15} | 15} contains {6}; similarly, (1°) X (6) does not contain 
any irreducible representation of depth 4, since none of the products {12} {3} {1}, 
112} {21} fi, {12k {15 {1} contains {6}. We see, in the same way, that the only 
irreducible representations of depth 5 which can appear in (1*) & (6) are (41) and 
(317), and so we know that 


bo 


we re, a 
— —_— 


(1%) x (6) = C4 (41) + C312(317) 
(31%) + 2(61) + (52) 
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Setting n = 2, we obtain cs: = 1, since (6) = I'(—4, 6) vanishes, as do all the 
terms on the right save (317) and (31°), which are —I'(2) and I'(2), respectively. 
Setting n = 3, we obtain, similarly, cq = 1. 

* Present address: 6202 Sycamore Road, Baltimore 12, Md. 
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AN APPLICATION OF ALGEBRAIC TOPOLOGY TO NUMERICAL 
ANALYSIS: ON THE EXISTENCE OF A SOLUTION TO THE NETWORK 
PROBLEM* 

By J. P. Roru 
UNIVERSITY OF CALIFORNIA, BERKELEY 
Communicated by G. C. Evans, May 18, 1955 


The electrical network problem was first treated comprehensively by Kirchhoff! 
in 1847. A complete proof for the existence of a solution was given by Hermann 
Weyl? in 1923 for the case of a purely resistive network when sources of electro- 
motive force are placed in series with the branches. This proof was elaborated 
by Eckmann® in 1945. A condition for the existence of a solution given by Synge,‘ 
attempting to cover a more general case, was unfortunately incorrect, as simple 
counterexamples show. In this paper we show the existence and uniqueness of a 
solution to the network problem under a condition which amounts to assuming that 
the dissipative power is positive definite. Since this condition is satisfied by any 
physically realizable network, it may be said that this result covers the general 
physical case (steady state). Actually, we state the electrical network problem in 
a purely algebraic-topological way. This is of interest since “electrical” networks 
are used to solve certain ordinary and partial differential equations. The results 
of this paper are essential to the proof of the validity of Kron’s method of tearing,® 
established in a second paper.® Our first task is to describe a mathematical model 
for the quantities and relations which exist in an electrical network. 

Let K be an electrical network: We shall consider K as an oriented one-dimen- 
sional complex. <A set of currents flowing through the branches of K may be con- 
sidered as the assigning of a complex number to each branch (‘“‘coil,” in Kron’s 
terminology). Hence such a set of currents will be treated as a vector or oriented 
l-chain. The space of such sets of branch currents thus coincides with the group 
('(K) of oriented 1-chains over the coefficient field of complex numbers. A mesh 
current, the current flowing around an oriented closed loop, corresponds to a 1- 
cycle, but, as we shall see, it is more appropriate to identify it with an element of 
the first homology group H'(K) of K. in fact, the space of such loop currents 





Sn a 


a 


ee 
a = 
ee 


ee eee 
a 











aoa 

















Vou. 41, 1955 MATHEMATICS: J. P. ROTH 519 
(they are also called ‘mesh currents”) is isomorphic with H'(K). With these 
identifications, the transformation 7 = C2’ used by Kron’ is the natural homomor- 
phism’ of H'(K) into C'(K). Here 2’, 7 are column vectors representing elements 
of H'(K), C'(K) with respect to bases chosen for these spaces. We recall, also, 
that C'(K) may be regarded as the relative homology group of K modulo its zero 
skeleton. Physically this transformation may be considered as an expression of 
Kirchhoff’s current law. It can also be seen that the space of node currents (node = 
vertex) is isomorphic to the group P°(A) of bounding zero chains. In Kron’s nota- 
tion the boundary operator is expressed by the equation J’ = A,J, where J and I’ 
are column vectors representing vectors in C'(A) and P°(K) with respect to bases 
chosen for these spaces. We shall use the notation more customary in topology, 
I’ = 0J, 0 being the boundary operator. In this context the boundary operator may 
be interpreted as an alternative expression of Kirchhoff’s current law, equivalent 
to the above statement of the same law. We may combine the mappings C and 
O to form the following “reduced” homology sequence of the pair (K, K°), K® being 
the zero skeleton of K: 


c 0 
Om H'(K) > C'(K) > PK) > 0. 
The result OC = 0 (or A,C = Q) is a consequence of the exactness of this sequence 
and is referred to by Kron as “the orthogonality condition.”’ Of course this result 
can also be proved directly very easily. 

The dual set of relations involved in considering the spaces of voltages and poten- 
tials will next be identified with the corresponding “reduced” cohomology sequence. 
We start at the opposite end of the sequence. First we identify the space of node 
potentials with the group Co(K) of zero cochains of K. The space of potential 
differences (node-pair potentials) coincides with the subgroup Pe(K) of Co(K), dual 
to P°(K), isomorphic to Co(A) modulo the subgroup of 0-cocycles. Po(K) is selected 
as the image of P°(K) under the following isomorphism: Let C° = )Cajo° be an 
oriented 0-chain of K, the o,° being oriented 0-cells of K, the summation running over 
the set of all 0-cells of A, and the a; being complex numbers. Let ¢(C°) be the 
0-cochain f such that f(o;) = a; The transformation ¢ is an isomorphism of C°(K) 
on ()(K). Let Po(K) be the image of P°(K) under ¢. Also, if {b,} is a basis for 
P°(K), we shall call the set { ¢(b,)} of images the “same” basis for P)(K). Now 
the space of branch voltages (“coil voltages” in Kron’s terminology) may be identi- 
fied with the group of 1-cochains of K, C\(K). The coboundary operator E = 6K’ 
(or E = AEF’, in Kron’s notation), where FL, /’ represent elements of C)(K), Po(K), 
can be considered as a statement of Kirchhoff’s voltage law. Finally, the space 
of mesh emfs coincides with the first cohomology group H,(K). Again, if we pick 
fixed bases for H'(K) and C'(K) and write a matrix equation for the natural homo- 
morphism (' and then use the “‘same”’ basis to express the dual homomorphism of 
C,(K) onto H,(K), then the mapping takes the form used by Kron, C,V = é. This 
mapping is another equivalent statement of Kirchhoff’s voltage law. These two 
mappings yield for us the following reduced cohomology sequence: 

Ct 6 
0+ H,(K) — C\(K) — P(K) = 0. 
By exactness, C,6 = 0 (or C,A = 0); this condition is termed by Kron ‘the ortho- 
gonality of mesh emfs and of node-pair potentials,” 
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The Twisted Isomorphism.—When steady-state conditions prevail, the differ- 
ential equations describing the various branch currents and branch voltages in the 
network may be reduced to a matrix equation V = LJ, a single equation for the 
entire network, where J is a column vector representing a 1-chain with respect to a 
basis chosen for C'(K) and V is a column vector representing a 1-cochain with re- 
spect to the “same” basis for C\(K). Finally, L, called the “impedance matrix,”’ 
is determined solely by the electromagnetic properties of the various coils (their 
resistances, mutual inductances, ete.) ; it is essential to realize that the transforma- 
tion ZL is independent of the manner in which the branches are hooked together, 
i.e., independent of the topological structure of K. In Kron’s notation L is de- 
noted by Z. 

The electrical network problem may be stated in physical terms as follows: 
Given a set of coils with given electromagnetic properties determining a trans- 
formation L, hooked together in a prescribed way to determine a network A and, 
given a set e’ of sources of emf impressed in series with the meshes (e’ is an element 
of H,(K)) and a set J’ of currents impressed from outside the system on the nodes 
of K (J’ is an element of P°(K)), the problem is to find the current flowing through 
and the potential drop across each coil—these subject to Kirchhoff’s laws. 

The “‘electrical”’ network problem may be stated in purely mathematical terms as 
follows: Given a complex K (determining the transformations ( and 0) and the 
matrix L, given e’ e H,(K), I’ « P°(K), the problem is to find V ¢ ()\(K) and J ¢ 
('(K) such that 1' V = LJ,2°C,V =e’ and3° oJ = 1’. 

Whether or not a solution to the network problem exists depends, of course, upon 
the nature of L. We shall say that a matrix L and the transformation represented 
by L is power definite if L + L, is positive definite. The dissipative power 1/2(V,J 
+ 1,J) is positive definite if and only if L is power definite. Hence, for any phys- 
ically realizable network, this condition is satisfied. 

THeorEM. J/f L is power definite, then the network problem has one and only one 
solution. 

The problem may be readily visualized by means of the following 
work diagram’: 


‘ 


‘algebraic net- 


¢ ° : 
0 ———> H'(K) —— C'(K) —— P*(k) —— 0 


L’ L Y x 
Cr | 5 
0 —— H,(K) ——— C,(K) ———- P,(K) -——— 0 


Here Y is the inverse of L, L’ = C,LC, Y’ = OY6. We will show that under the 
conditions on L the mappings L’ and Y’ are one to one, that is, L’ and Y’ are non- 
singular. We will then exhibit the solution. For suppose that L’7’ = 0, and let 7 = 
Ci’. Then i,’L’7’, which equals 7,11, is zero, and so, too, is i,L,7 zero. Thus i,(L + 
L,)i = 0, and so, by the power-definiteness assumption on L, i must be zero. But 
since, by exactness, C is an isomorphism into, 7’ must also be zero. Thus the kernel 
of L’ is zero, and, since the spaces have the same dimension, L’ must be one to one. 
By a similar argument, Y’ is also shown to have an inverse. (By exactly the same 
argument one may show that L’ and Y’ have inverses if 1 — Z, is assumed to be 
definite. ) 
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One may easily verify that J = C(L’)—'e’ + Y6(Y’)-U’ and V = LJ will con- 
stitute a solution. Next we consider the question of uniqueness. 

Suppose, then, that V, J and V%*, /J* are solutions to the given problem. Let 
V4 = V — V* and J/# = J — J*. Then 0J/# = 0, and so, by the exactness of the 
homology sequence, there is an 7# in H'(K) such that J/# = Ci#. But V# = LJ#, 
and, since C,V# = 0, we have C,LCi# 0. By the above proof, 7# and hence J# 
must be zero, sothat J = J/* and V = J 

Of course, to solve the network problem, to obtain the general solution, it is not 
necessary to invert both L’ and Y’. To solve, for example, by means of inverting 
L’ only, we proceed as follows: Let C'(K) be written as the direct sum of the 
image ./'(K) of H'(K) under C (./'(K) is the group of cycles) plus another group 
P'(K) (ef. companion paper,® Sec. II). Then J, an element of C'(K), can be 
written as J = 7+ /,7¢€ M'(K),/ ¢€ P'(K). Thus I is chosen such that 0J = I’ 
and may be considered as known. Then V = L(¢+J/). Nowe’ =€C,V =C,L(G+ 
I). But 7 may be written as the image of an element of H'(K), 7 = C7’, 7’ « H'(K). 
Thus C,LCi’ = e’ — C,LI, so that 7’ = (L’)—'(e’ — C,LI). Thus J is determined, 


7# 


J = Ci’ + I, and V = LJ, so that the solution is obtained by inverting L’ only. 


Similarly, one may obtain the solution by inverting Y’. We also have the follow- 
ing consequence of the above theorem. 

CoroLuary. Jf a network problem is physically realizable, then there exists one 
and only one solution. 


* Presented before the American Mathematical Society, December 27, 1954, under the title 
“An Application of Algebraic Topology to Numerical Analysis: Kron’s Method I.” 
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THE ROLES OF THE OPTIC VESICLE AND OTHER HEAD TISSUES 
IN LENS INDUCTION 


By ANTONE G. JACOBSON 
DEPARTMENT OF BIOLOGICAL SCIENCES, STANFORD UNIVERSITY 
Communicated by V. C. Twitty, May 12, 1955 


The induction of the amphibian lens, although one of the most extensively in- 
vestigated phenomena in amphibian development, has continued to present per- 
plexing if not paradoxical features. The early work of Spemann and of Lewis 
established that the optic vesicle has an active role in lens induction.! In some 
species, however, it was found that lenses could be formed in the absence of the optic 
vesicle. While Spemann had found this to be the case in Rana esculenta,' others, 
notably Woerdeman,’ were unable to obtain independent lenses in esculenta. 

Studying determination of ciliary polarity in embryos of Amblystoma punctatum, 
Twitty® observed that low temperature slows the process of determination less than 
it retards the rate of morphological development. Following this lead, Ten Cate,‘ 
working in Woerdeman’s laboratory, has shown that the difference between Woerde- 
man’s and Spemann’s results with esculenta is attributable to the lower tem- 
perature at which Spemann kept his embryos prior to the extirpation of the retinal 
anlage. There still remained the problem of identifying the inductors involved in 
evoking a lens in the absence of the optic vesicle. 

Liedke® has demonstrated in A. punctatum that “lens competence is acquired 
during the open neural fold stages through activation by the anterior ento-meso- 
derm.”’ Considering together the effects of temperature and the role of the ‘an- 
terior ento-mesoderm,”’ Twitty® has suggested that 


. . . probably in all species lens formation is ordinarily elicited by the action of at least two 
supplementary or mutually reinforcing inductors, the head mesoderm and the retinal rudi- 
ment. Under average environmental circumstances and in most species the component 
events of development are geared with one another in such manner that the epidermis is 
insufficiently differentiated at the neurula stage to be fully receptive to the action of the 
mesoderm, and requires a subsequent and final impetus from the optic vesicle. At low tem- 
peratures, however, and especially in species like R. esculenta in which such treatment causes 
morphological development to lag substantially behind chemical differentiation, the lens 
epidermis has already become sufficiently mature and responsive at neurulation to require 
no further incitement beyond that provided by the head mesoderm. 


This suggestion has been substantiated in the present study by combining on a 
single species, 7'riturus torosus, the temperature and “head mesoderm”’ experi- 
ments and also by the results of isolating from the young neurula presumptive lens 
epidermis together with the subjacent tissues in vitro. 

As shown in Figure 1, the temperature at which torosus embryos are reared prior 
to excision of the retinal anlage in the open neural plate stage has a pronounced 
effect on the percentage of animals which develop lenses in the absence of the optic 
vesicle. An example of an embryo with an independent lens is shown in Figure 2. 

Liedke® had found in A. punctatum that gastrula epidermis transplanted to the 
lens site develops a lens only if it is in position during neurulation. Lenses de- 
veloped in the graft if the host was an early neurula but not if the host was a late 
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neurula. In the present experiments with torosus, this same distinction was evident 
but not so clear cut. Lenses developed in 71 per cent of the grafts to early neurula 
hosts and in 28 per cent of the grafts to late neurula hosts. 
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ait Fic. 1.—Percentage of torosus embryos which develop a lens in 
| the absence of the optic vesicle after rearing the groups of embryos 
H,) at the indicated temperatures prior to excision of the retinal rudiment 
y |. in the early neurula (stage 16). All groups were kept at 19° C. after 
Wg the operation. 
im 
iM In the early neurula the presumptive lens epidermis is underlaid directly by the 
H entodermal wall of the archenteron; the anterior edge of the lateral plate meso- 
} derm lies just posterior to the lens site. Presumptive lens epidermis from young 
an neurulae was isolated in combination with (1) the subjacent entodermal wall of the 


| archenteron; (2) the anterior portion of the lateral plate mesoderm; and (3) these 
two tissues combined (Fig. 3). Lenses were induced in a very substantial percent- 
\ age of cases in ail three series. No lenses have been found in presumptive lens 














Fic. 2.—Transverse section of a 12-mm. torosus larva in which a lens has developed in the 
absence of the retina. The embryo was reared at 13° C. until the early neurula stage, when the 
right retinal anlage was excised. 

Fig. 3.—Section of an explant of presumptive lens epidermis combined with the subjacent ento- 
dermal wall of the archenteron and the anterior portion of the lateral plate mesoderm; taken 
from a young torosus neurula (stage 16). A lens has been induced from the epidermis. 




















lO ——— — 


Vou. 41, 1955 ZOOLOGY: A. G. JACOBSON 525 


epidermis isolated alone. These explants were taken from embryos which had 
been reared at 9° or at 13°C. 

The results of these experiments indicate that lens induction is the result of the 
synergistic action of the optic vesicle and of the nonneural tissues underlying the 
presumptive lens epidermis during neurulation. In fact, these latter tissues may by 
themselves induce lenses, without the action of the optic vesicle, especially in 
embryos that have been reared at appropriate temperatures. The emergence of the 
lens under the influence of more than one inductor finds a close parallel in the in- 
duction of other head organs, such as the otocyst, the olfactory organ, and the hy- 
pophysis.’ 

'H. Spemann, Embryonic Development and Induction (New Haven: Yale University Press, 
1938), chap. ill. 
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